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1.0 INTRODOCTION 



1.1 ACID MINE DRAINAGE 

Acid mine drainage (AMD) , from coal mining in the East and 
both coal and hard-rock mining in the West, is one of the most 
persistent environmental pollution problems in the United States. 
Acid water draining from these sites lowers ground and surface 
water quality, impacts aquatic and terrestrial biota, and 
degrades domestic water supplies. An estimated 11,000 miles of 
rivers and streams in Appalachia have been negatively impacted by 
AMD (Brodie et al. 1986). More than 5,000 miles of waterways in 
the East and mid-West fail to meet federal water quality 
standards owing to acidic mine water (Kim et al. 1982). In 
Colorado alone, an estimated 10,000 abandoned and inactive metal 
mines are point sources of AMD, and at least 25 watersheds and 
450 stream miles have been rendered barren of aquatic life 
(Guertin et al. 1985). 

Artificial wetlands have been proposed as an approach to 
provide low-maintenance, broadly applicable solutions to AMD 
problems and numerous state agencies are currently investing 
considerable sums of money in the construction of wetlands to 
treat AMD. The efficacy of these projects is variable (Burr is 
1984) . Some wetlands improve water quality to meet drinking 
water standards. Others are less effective and provide only 
temporary metal uptake, selective metal uptake, seasonal 
effectiveness, or even complete abatement failure. Poor wetland 
performance may be attributed to inappropriate hydrology or 
vegetation. Flow channeling may cause polluted discharge to 
bypass remediation stages; transplanted aquatic plants may be 
unable to tolerate AMD loading or to take up sufficient metals to 
significantly abate AMD (McHerron 1985, Wieder and Lang 1986) . 
However, some wetlands with improved hydrologic conditions and 
well-established metal-tolerant plants still do not perform well, 
suggesting that other factors, such as microbial activity and 
sediment geochemistry, are of overriding importance at some sites 
(Kleinmann and Girts 1986, Lang and Wieder 1985, Peccia & 
Associates 1987) . 

The wetlands characterized in this study afforded a unique 
opportunity to observe the effects of more than 55 years of 
mining related effluent and AMD upon a natural wetland (Swamp 
Gulch wetland) . In addition, this wetland system could be 
compared to a non-polluted system, and to a recently constructed 
AMD wetland treatment system located near Sand Coulee, Montana. 
The findings will be valuable in the design of future artificial 
AMD wetland treatment systems by the Montana Abandoned Mine 
Reclamation Bureau. 



1.2 STUDY OBJECTIVES 

Specific objectives of this study were to: 

1) determine the areal and vertical gradients of metal 
loading in the natural wetland; 

2) assess contaminant transit time in the natural wetland; 

3) quantitatively describe plant associations and 
determine production by dominant species in the 
natural wetland; 

4) assess metal loading of dominant vegetation and potential 
impacts to wildlife in the natural wetland; 

5) estimate the effective longevity of the wetland to 
control the acid mine drainage in the natural wetland; 

6) correlate wetland mass and volume as a function of the 
quality and quantity of acid mine drainage in the natural 
wetland; and 

7) compare the relative biological activity of the natural 
wetland to the artificially constructed wetland. 



1.3 CONCLUSIONS 

1.3.1 Wetland Hydroaeochemistry 

1. The Swamp Gulch natural wetland is characterized by three 
basic stratigraphic units: 1) the surficial material 
(average thickness 0.6 m) composed of undecomposed organic 
debris, referred to in this report as the acrotelm; 2) the 
intermediate gray muck zone (average thickness 2.2 m) 
consisting of decomposed organic materials, silt and clay 
with occasional sandy gravel lens, generally referred to in 
this report as the catotelm; and 3) the underlying alluvial 
gravel zone. 

2. Trace metal concentrations were highest in the acrotelm near 
the Swamp Gulch discharge point and decreased both laterally 
and vertically away from this point. 

3. The only catotelm interval notably impacted by the acid mine 
drainage was near the discharge point, with trace metal 
levels in the remaining catotelm being approximately equal 
to measured background concentrations. 
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4. Trace element concentration distribution in the underlying 
alluvial gravel material showed only a nominal association 
with the Swamp Gulch acid mine discharge. 

5. The present iron treatment area of acid mine drainage is 
confined to approximately 3.9 ha (9.6 acres) and represents 
an acrotelm volume of about 23,000 m^ with an average 
thickness of 0.6 m. This volume treats a daily acid mine 
drainage inflow of 45 m-^ containing a volume weighted 
average of 21 mg/L iron, 3.9 mg/L manganese, 6.5 mg/L 
aluminum, 2.7 mg/L zinc and 0.71 mg/L copper at a pH of 4.0. 
The system has in the past, been effective in reducing iron 
levels in effluent waters. 

6. Approximately 550 metric tons of iron have been deposited in 
the wetland acrotelm from the Swamp Gulch source during the 
past 55 years. Under present (1987) loading rates (338 
kg/year) this represents a 1600 year accumulation, which 
strongly suggests past loading rates were much higher than 
those observed during the study period. 

7. Conditions conducive to the physio-chemical precipitation of 
iron may not be conducive to creating optimum conditions for 
microbial buffering of pH and sulfate reduction. This may 
explain why Swamp Gulch cores in the laboratory were 
ineffective in removing iron but, at the same time, the 
above mentioned iron was present in the acrotelm. 

8 . Lead is apparently effectively removed from the AMD within 
the wetland. The AMD attributed acrotelm lead mass is 1.3 
metric tons which represents 2,350 years of deposition at 
present loading rates. 

9 . The acrotelm copper mass due to AMD is approximately four 
metric tons which represents 325 years at the observed 
annual loading rate. This metal is apparently only 
partially ameliorated within the wetland. 

10. Zinc, manganese and cadmium acrotelm masses due to AMD are 
6.0, 1.1 and 0.004 metric tons, respectively. These 
represent deposition intervals of 137, 16.9 and 6.7 years, 
respectively. It is clear that the wetland efficiency in 
removal of these metals is limited. 

11. Using lead as a basis of comparison of the Swamp Gulch 
wetland efficiency in removing metals from AMD, the removal 
rates are: iron, 70 percent; copper, 14 percent; zinc, 5.8 
percent; manganese, 0.7 percent; and cadmium, 0.3 percent. 



12 . Bromide was a good choice for tracing water movement in the 
natural wetland studied. Employing a natural flow tracer 
technique and the bromide specific ion electrode method of 
analysis, sufficient bromide concentrations were detected. 

13. The most representative flow velocity of contaminated Swamp 
Gulch water in the wetland is that determined from a bromide 
tracer test. The tracer determined mean flow velocity is 
2.6 X 10"2 cm/s (73.7 ft/day) and ranged from 3.1 x 10"^ to 
6.2 X 10"^ cm/s. These estimated velocities are quite 
variable and represent a low flow regime (Autumn) . 

14. Field hydraulic conductivity (K) tests were performed on 
auger holes and wells in the area of expected metal loading 
and tracer movement. These hydraulic conductivity estima- 
tions ranged from 1.0 x 10"^ to 1.5 x 10"^ cm/s with a mean 
of 8.3 x 10~4 cm/s. The mean water flow velocity determined 
from the field hydraulic conductivity values is 7.1 x 10"^ 
cm/s. 

15. The large difference (three orders of magnitude) between 
mean hydraulic conductivity determined flow velocity and 
tracer determined flow velocity infers the measurement of 
two different flow systems. The tracer velocity represents 
preferential flow paths (channelling) within the shallow 
flow system (acrotelm) and the hydraulic conductivity 
determined velocity represents that flow found at greater 
depth (catotelm) . 

16. The distance the AMD water travels, before problem ameliora- 
tion, as estimated by the iron impacted area, is 200m. 
Using this 200m distance and the average tracer velocity, 
the average residence time of AMD water in contact with the 
wetland soils is estimated to be nine days. 

1.3.2 Wetland Veaetation/Wildlife 

17. The Swamp Gulch wetland was dominated by sedge ( Carex 
rostrata) . followed by birch ( Betula alandulosa var hallii ) , 
lodgepole pine ( Pinus contorta) , willow ( Salix boothii ) and 
itioss ( Sphagnum tenellum ) with percent canopy coverage of 
55.5, 11.3, 6.3, 5.5 and 4.7, respectively. Mesic sites 
were dominated by nearly pure stands of Carex rostrata while 
the more xeric sites were characterized by Douglas fir 
( Pseudotsuqa menzeisii ) , Englemann spruce ( Picea enael- 
mannii) and subalpine fir ( Abies lasiocarpa ) and other 
species typical of the surrounding upland forest. 

18. The composition of plant specieiB in the Swamp Gulch wetland 
was similar to that of the Hardscrabble Creek background 
wetland . 



19. Above-ground production at the Swamp Gulch wetland (3820 
kg/ha) was not significantly different than the Hardscrabble 
Creek background site (3750 kg/ha). 

20. Metal levels were generally elevated above background in the 
tissues of the dominant plant species. Levels were usually 
higher in below-ground rather than above-ground plant 
material . 

21. Carex rostrata accumulated elements in the following order: 
Cu>Pb>Al>Fe>Cd>Zn. This species apparently possesses an 
exclusionary mechanism which prevents the absorption of Mn 
and Ni . 

22. Of the shrubs species studied, Salix boothii accumulated 
elements to a greater concentration than did Betula 
glandulosa . Manganese and Ni were readily absorbed by these 
species as were Cd, Cu, Fe, Pb, Mn, Ni and Zn. 

23. The bryophyte Isopterygium pulchellum had tissue metal 
concentrations that exceeded background levels for Al, Cd, 
Cu, Fe, Pb, Mn, Ni and Zn. 

24. Some metal concentrations exceeded the maximum tolerable 
dietary intake levels for domestic animals, suggesting a 
potential problem for animals that consume these plants. 
The greatest risk was from high levels of Cd. 

1.3.3 Wetland Microbiology 

25. Removal of acidity in Swamp Gulch natural wetland sediment 
occured through the action of anaerobic microorganisms and 
was not solely a function of abiotic soil characteristics. 

26. Based on laboratory and greenhouse experiments, acidity 
remediation occured in the Swamp Gulch natural wetland but 
not in the Sand Coulee constructed wetland. 

27. Sul fate-reducing bacteria were present at all experimental 
sites, and appeared to be reducing sulfate in the Swamp 
Gulch sediments, but not the Sand Coulee wetland sediments. 

28. Adequate carbon and energy levels to support the growth of a 
variety of indigenous microorganisms including sul fate- 
reducing bacteria are evident in the Swamp Gulch sediments. 

29. In the laboratory experiments, iron removal from AMD was not 
evident in either the Swamp Gulch or Sand Coulee wetland 
cores . 



30. The total microbial biomass and densities of sul fate-reduc- 
ing bacteria were low in the Sand Coulee constructed wetland 
sediments, probably because of insufficient levels of carbon 
and other nutrients. 

31. The microbial populations at both Swamp Gulch and Sand 
Coulee wetlands were physiologically stressed. This may 
account for the inability of the Sand Coulee sediments to 
facilitate acid, iron or sulfate removal. Furthermore, it 
may be that the Swamp Gulch sediment microbiota are not able 
to achieve their maximum physiological capabilities for AMD 
amelioration. 

32. It should not be automatically assumed that poor performance 
in a constructed wetland is due solely to wetland size or 
hydrology. There is a significant microbiological contribu- 
tion in working wetlands that should be transferable to 
non-working wetlands. 

1.4 RECOMMENDATIONS 

1.4.1 Wetland Hydroaeochemistry 

1. These data suggest that any constructed wetland should be 
made as large as practically possible. It is doubtful that 
any constructed wetland could be built as large as the 
impacted area of this natural wetland (3.9 ha). 

2. The catotelm (muck zone) was generally ineffective in 
removing metals from AMD, but a limited thickness of this 
material (probably less than 7 cm) may be needed for a 
proper siobstrate for some types of wetland vegetation. 

3. The thickness of the natural wetland acrotelm averages 0.6 
m. This thickness may be desirable in constructed wetlands 
operating under similar circumstances. 

4. The specific ion electrode method of bromide tracer analysis 
could be improved by periodic double sampling and analysis 
with a different method (titration) , especially at low 
bromide concentrations. 

5. The acquisition of limited additional data from the Swamp 
Gulch wetland site would allow a determination of the 
system's present efficiency in removing metals from AMD. 
This will help determine the effective treatment longevity 
of the wetland. 



6. Precipitation of iron by physio-chemical processes may be a 
very important mechanism for removing iron from AMD (Wieder 
and Lang 1986) and likely played an important role in the 
Swamp Gulch wetland. Further evaluation of the wetland iron 
is recommended to determine if this material is organically 
derived or inorganic oxide precipitation. 

1.4.2 Wetland Vegetation/Wildlife 

7. Carex rostrata . Betula alandulosa and Salix boothii should 
be strongly considered for transplantation to constructed 
wetlands. These plants absorbed substantial quantities of 
Al, Cd, Cu, Fe, Pb, Mn, Ni, and Zn without ill effect. These 
plants are common to many parts of Montana and should be 
relatively easy to transplant. It is therefore recommended 
that a screening study be implemented to determine the 
ability of these species to tolerate transplantation and to 
proliferate at a constructed wetland. 

8. Potential impacts to wildlife from high plant tissue metal 
levels should be investigated. The concern has been 
expressed that in constructing wetlands for the amelioration 
of AMD one is increasing the risk of poisoning the wildlife 
resource. 

9 . The best indication of whether animals are at risk at the 
Swamp Gulch wetland is to analyze tissue samples from those 
animals that are most likely to be poisoned. It is therefore 
recommended that voles r Microtus spp.) be trapped and their 
kidneys be analyzed for metals, especially Cd. 

10. Until it can be demonstrated that animals will not be 
poisoned by consuming vegetation at constructed AMD 
wetlands, it is recommended that these systems be fenced to 
limit wildlife access. 

1.4.3 Wetland Microbiology 

11. Procedures need to be established to transfer microbio- 
logical activity of a working wetland to an ineffective 
wetland. 

12. The ester-linked phospholipid fatty-acid technique is an 
important tool for assessing the impacts of AMD on in- 
digenous microbiota and in evaluating the efficacy of 
constructed wetlands. 

13. The addition of carbon and other nutrient sources to 
wetlands would diversify the microbial community structure 
and be especially beneficial to the sul fate-reducing 
bacterial populations. 



14. Wetland constiruction practices should provide as much 
cation exchange capacity as possible to aid biological 
removal acidity and metals. Acidity control through the 
application of appropriate amendments should also be 
considered. 

15. Wetland construction practices should provide favorable 
conditions, such as anaerobiosis (e.g., by varying con- 
structed wetland depth) , for the proliferation of sulfate- 
reducing bacteria. If these bacteria are able to reduce 
sulfate to sulfide, metals should be immobilized through the 
formation of metal sulfide precipitates. 

16. Seasonal effectiveness of wetlands needs to be assessed in 
the context of the microbiota and site geochemistries. 

17. Emphasis on the exact nature of physiological stresses 
imposed on microbial populations in wetlands need to be 
elucidated . 

18. Vertical and lateral zones of microbial activity need to be 
mapped in both natural and constructed wetlands in order to 
develop engineering criteria for wetland construction. 



2.0 NATURAL SYSTEMS FOR AMD CONTROL 

2 . 1 WETLAND TREATMENT SYSTEMS 

Research has confirmed the theory that natural wetlands can 
improve the quality of acid mine water flowing through them 
(Wieder and Lang 1984, 1986). However, since natural wetlands 
are seldom located adjacent to acid mine effluent and because of 
legal barriers regarding the discharge of pollutants into 
wetlands, they are seldom utilized for AMD abatement. Construct- 
ing artificial wetlands is the obvious alternative to using 
natural wetlands. 

Artificial systems, both in the laboratory and in the field, 
have demonstrated the ability to improve water quality by 
removing certain metal ions (Kleinmann et al, 1983, Burris et al. 
1984, Gerber et al. 1985, Wieder et al. 1985 and Brodie 1986). 
These systems utilize bryophytes, particularly Sphagnum spp., 
which have generally been successful in removing Fe, Mg, Mn and 
other metals from solution. Simola (1977) and Chaney and 
Hundemann (1979) studied the ability of Sphagnum to absorb and 
tolerate Cd. Bryophyte species of the genera Rhvnchostegium . 
Rlytidiadelphus ^ Amblvstegium and Fontinalis have demonstrated 
metallophytic properties (Wehr and Whitten 1983, Brown and 
Beckett 1985) . Emergent vascular plant species (Typha. Carex, 
Eleocharis and Juncus ) and floating plants f Eichhornia . Salvenia 
and Lenuma) frequently occur in areas with acidic water and/or 
high metal concentrations. Haag and Covert (1987) reported Typha 
angustifolia and several grass species growing in AMD water in 
Missouri. Dinges (1982) reviews natural and artificial systems 
for improvement of polluted waters. Species of Typha are able to 
tolerate high metal levels, and therefore have been extensively 
studied (Taylor and Crowder 1983a, 1983b, 1984, McNaughton et al. 
1974, Mayer and Gorham 1951, Boyd 1970 and Bayly and O'Neill 
1972). In terms of biofiltration of AMD, vascular plants have 
not been studied as intensively as bryophytes. The ability of 
vascular plants to thrive in AMD water indicates that more 
research is needed in this area. 

Mechanisms thought to operate in pollution control include 
cation exchange of metal ions on the surface of organic matter, 
physical filtration, chelation or complexation by organic 
compounds, dilution, and biological sorption by bacteria, algae, 
fungi, animals and plants. Guertin et al. (1985) presents a 
review of each of these mechanisms. 

Field testing of artificial wetlands are underway at over 20 
sites in the eastern United States (Grits and Kleinmann 1986) , 
with the possibility of hundreds of wetlands being constructed in 



Appalachia in the next few years. Guertin et al. (1985) discuss 
a wetland system installed in Colorado and Heil and Kerins (1988) 
report results from an artificial wetland in central Montana. 

Primary design criteria for construction of artificial 
wetlands are based on estimates of system loading rates, since 
these directly effect maximum treatment efficiency and capacity 
of the wetland. Influent volxime and chemical composition in 
conjunction with treatment rate and required effluent concentra- 
tions provide the basis for determining the optimal AMD to 
wetland volume ratio, flow path length, and detention time. 
Baseline data on flow rate and chemical composition need to 
reflect yearly extremes. For a review of construction parameters 
and essential baseline information, consult Girts and Kleinmann 
(1986), Pesavento (1984) and Brooks (1984). Girts and Kleinmann 
(1986) present important criteria used for selecting, planting 
and transplanting various plant species. 



2.2 BIOSORPTION OF ELEMENTS BY PLANTS 

The idea that plant species absorb nutrients from the soil 
was first suggested by Glauber in 1656. Prior to this, it was 
generally accepted that plants accpaired all their sustenance from 
air and water. In 1804, de Saussure demonstrated that plants 
acquire a portion of their constituents from the medium in which 
they are rooted. Hewitt (1975) reported that higher plants 
absorb a portion of all the elements that are present in the 
rhizosphere solution, but exhibit absorption selectivity. 

Vascular plants characteristic of soils with high metal 
concentrations have been known for centuries. In 1588, Thalius 
noted that Minuartia verna was an indicator of metal presence in 
soil (Ernst 1965) . Metallophytes (plants that accumulate metal 
ions) often have tissue concentrations of ions several times 
higher than the medium in which they are growing. These types of 
plants have been used in geologic prospecting to indicate ore 
bodies (Cannon 1960) . 

Some bryophytes have also been identified as metallophytes, 
accumulating elements well in excess of their structural or 
metabolic requirements. The high ion exchange capacity, simple 
cellular morphology, and compact low growth form of many mosses 
can result in the rapid accumulation of metal ions, either by 
extracellular ion exchange or by entrapment of particulates. Lee 
et al. (1983) studied the brilliant lime green mosses that 
flourished in hillsides blanket bogs fed by spring water enriched 
with Zn, Pb and Cu. Bryophytes belonging to the Bryaceae and 
Bartramiaceae families were reported useful in identifying Cu and 
U deposits (Whitehead and Brooks 1969) . 
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2.3 METALLOPHYTES AND LAND/WATER RECLAMATION 

Metallophytes have proven useful in the reclamation of 
metalliferous waste sites (Gadgil 1968, Johnson et al. 1977) , and 
in the treatment of wastewater (Dinges 1982). The mechanisms of 
wastewater treatment by natural systems is discussed in detail by 
Dinges (1982) . The use of metallophytes to improve the quality 
of acid mine water is being studied extensively throughout the 
United States. 

2.4 CHEMICAL PROCESSES OF ACID MINE DRAINAGE 

2.4.1 Acid Water Formation 

Acid water occurs when pyrite (FeS2) is oxidized in air and 
water to form sulfuric acid (H2SO4) and iron hydroxide (Fe(0H)3). 
Jaynes et al. (1984) and Guertin et al. (1985) explain the 
stoichiometry of acid formation in detail. The overall chemical 
reaction was summarized by Schumate et al. (1971) as: 

FeS2 + 3.75 O2 + 3.5 H2O > H2SO4 + Fe(0H)3 

Acidic water is often indicated by the presence of a white to 
reddish iron hydroxide precipitate known as 'yellowboy'. 

2.4.2 Bacterial Mediation 

The kinetics of acid formation are dependent on the 
availability of oxygen, the surface area of the pyrite, the pH of 
the influent water, and the activity of iron-oxidizing bacteria. 
The principal bacteria involved in accelerating pyrite oxidation 
are the chemoautotrophic Thiobacillus thiooxidans (Kleinmann and 
Crerar 1979). Since these bacteria can increase the rate of this 
reaction by a factor of over 100 times, they can play a crucial 
role in acid production. Singer and Strum (1970) concluded that 
this bacterial ly mediated reaction was the rate determining step 
in acid formation. Walsh and Mitchell (1972) believe that the 
production of acid can be controlled by breaking the Fe-bacteria 
succession in the pH range of 4.0 to 4.5. Bacteria of the genus 
Metal loaenuiro have been found to catalyze the acid producing 
reaction (Parizik 1985) , while other microorganisms such as Tjl 
thiooxidans are important in bringing iron into solution from 
soil and organic matter (Oborn and Hem 1962). 

2.4.3 Influence of pH 

In general, metal ions increase in solubility several orders 
of magnitude with each unit decrease in pH. The solubility and 
stability of metals in solution is also a function of Eh, or 
oxidation-reduction potential . Hem and Cropper (1959) present 
data on the relationship between Fe ions, pH and Eh. In terms of 
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AMD treatment systems, Eh-pH diagrams are useful for obtaining 
rough estimates of the dominant ions in the effluent (Guertin et 
al. 1985) 

Traditional methods for treating acidic water are to 
oxygenate the water to remove volatile acids such as carbonic 
acid. Bases such as sodium hydroxide, limestone or sodium 
carbonate are often used as neutralizing agents. Injecting 
alkaline materials into acid producing spoil material to prevent 
AMD from occurring was reviewed by Ladwig et al. (1985). 

2.4.4 Pvrite Oxidation Control 

Methods for controlling pyrite oxidation include the use of 
bactericides (Barnes and Romberger 1968, Shearer et al. 1970), 
the use of bromo-substituted phenols to inhibit metabolism 
(Aleems 1972) , the introduction of heterotrophic bacteria to prey 
on the iron and sulfur oxidizing bacteria, removal of attachments 
for the stalked bacterium Metal loaenium . and increase the pH in 
the neighborhood of the bacteria (Walsh and Mitchell 1972). Some 
of these techniques have been tested in the laboratory while 
others have merely been suggested. 

Anionic surfactants, primarily sodium dodecyl sulfate (SDS) , 
alters the cytoplasmic membrane of the bacteria, allowing 
hydrogen ions access to the cell. Two organic compounds, sodium 
benzoate and potassium sorbate, have been demonstrated as 
effective as SDS (Erickson et al. 1985). Pelletized bactericide 
has raised the pH of spoil material, allowing it to be revege- 
tated (Zaburnov 1987) . 
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3.0 SITE DESCRIPTION 

3.1 LOCATION 

The Swamp Gulch Wetland study site is located 14 miles (22.4 
km) east of the town of Lincoln, Montana in Lewis and Clark 
County (section 2 0, TIN, R6W) (Figure 3-1). The site is bordered on 
the north by Montana highway 200 and on the south by the 
Blackfoot River. The entire wetland covers approximately 72.9 ha 
(180 acres) ; however, the study area proper covers about two 
thirds of that area of which 3.9 ha (9.6 acres) is apparently 
impacted by AMD. 

3.2 NATURAL ENVIRONMENT 

The wetland study site is in the headwaters of the Blackfoot 
River at an elevation of about 1573 m (5160 ft) . Elevations in 
the Swamp Gulch drainage which feeds a portion of the study area 
range from about 1574 m (5165 ft) to 1908 m (6260 ft) . The Swamp 
Gulch watershed contains 72.6 ha (179 acres) and is composed of 
generally well vegetated southeast to southwest aspects with a 
total length of about 1.6 km (1 mi). The site is on the west 
side of the continental divide and is therefore influenced more 
by Pacific weather patterns than by the continental climate in 
the region east of the divide. However, moisture from the Gulf 
of Mexico occasionally spills over the divide during the summer 
producing intense thunderstorms. Mean annual precipitation is 
approximately 73.7 cm (29 in) (SCS 1981) much of which occurs as 
snow (Coffin and Wilke 1971) . Precipitation during the 1987 
water year and annual year averaged 60 percent of normal at the 
Lincoln Ranger Station. Record (1949-1987) maximum and minimum 
annual precipitation at the Lincoln Ranger Station are 91.49 cm 
(36.02 in) and 31.14 cm (12.26 in), which occurred in 1975 and 
1973, respectively. Average January and July temperatures are - 
7.2°C (19.1°F) and 16.4°C (61.6°F), respectively (NOAA 1987). 

The Blackfoot River drainage was significantly influenced by 
glaciation. Sediments deposited by retreating valley glaciers 
created a series a moraines and outwash plains in the Lincoln 
area, allowing for the formation of wetlands in drainage impeded 
areas on the upslope sides of the moraines (Alt and Hyndman 
1972). The lithology of the study area is dominated by the 
Spokane Shale Formation, part of the Belt series super group. 
This formation consists of red or red-purple shale with numerous 
green beds locally with quartzite (Ross et al. 1955). Parent 
material of the valley floor apparently consists of a thin layer 
of glacial till with fluvial gravel deposits. 
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Swamp Gulch and many other drainages of the region have been 
impacted by mining. The Carbonate Mine is within the Heddleston 
mining district and in 1933 included an adit about 850 feet long 
and a shaft approximately 300 feet deep (Pardee and Schrader 
1933) . A 75 ton mill was erected in the summer of 1944 and it is 
reported that 2500 tons of lead-zinc ore was taken out in 1945 
(Shea 1947) . 

The geology of the Carbonate Mine is similar to that of the 
Mike Horse Mine which is the largest mine in the district (Shea 
1947). The mine workings are in a diorite sill that intruded 
argillite of the Spokane Formation, Belt Precambrian sedimentary 
rocks (McClernan 1983) . Intruding into the sill and Belt rocks 
is a decomposed rhyolite dike (Shea 1947) . Pardee and Schrader 
(1933) describe the lode as consisting of a sheared and altered 
diorite partly replaced by sulphides. The mine lode contains 
abundant pyrite in several forms including fine grained, coarse 
grained and radial. Ore minerals are galena, sphalerite, silver 
and gold. Gangue minerals occur primarily as sericite, pyrite 
and quartz and also include "a carbonate containing calcium, iron 
and a little manganese" (Pardee and Schrader 1933) . 

The dominant soil in the study area is a Hydric Borofibrist 
(Soil Survey Staff 1975) . The gravels most often encountered at 
depth are fluvial in origin as evidenced by many subrounded 
pebbles. These were deposited by the meandering river following 
glaciation and were likely derived from glacial outwash mate- 
rials. 

Vegetation within the Swamp Gulch wetland is dominated by 
sedge (Carex rostrata) . This species is sometimes associated 
with shrubs ( Betula qlandulosa and Salix boothii ) or lodgepole 
pine ( Pinus contorta ) , but is most often encountered in nearly 
pure stands. Where the water table is relatively low, plant 
species more indicative of the surrounding forest are found. 

Upland forest vegetation is predominantly subalpine fir 
( Abies lasiocarpa ) , Douglas-fir ( Pseudotsuqa menzeisii) and 
ponderosa pine ( Pinus ponderosa ) climax forest (Ross and Hunter 
1976) . Understory species include pinegrass ( Calamaarost is 
rubescens ) , huckleberry f Vaccinium s pp. ) , bluebunch wheatgrass 
(Agropyrgn spicatum ) , mountain brome ( Bromus marainatus ) and 
Columbia needlegrass ( Stipa columbiana ) . 

3.3 NATURAL WETLANDS 

3.3.1 Definition and Classification 

Marshes, swamps and bogs have been well known terms for 
centuries, but only recently have attempts been made to group 
these ecological systems under the single term "wetland" . This 
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change has developed out of a need to understand and describe the 
characteristics and values of all types of land so that they can 
be wisely and effectively managed (Cowardin et al. 1979). 

In general, wetlands are lands where the soil or substrate 
is at least periodically saturated with or covered with water. 
Saturation with water is the dominant factor in soil development 
and in determining the types of plants and animals that are in 
the soil and on its surface. 

The most widely used wetland classification system in the 
United States is that of Martin et al. (1953), which was 
republished in the U.S. Fish and Wildlife Service Circular 39 
(Shaw and Frederick 1956) . Wetland types were based on criteria 
such as water depth, water permanence, water chemistry, vegeta- 
tion life form and dominant plant species. The system defined by 
Golet and Larson (1974) refined freshwater wetland types of 
Circular 39 for the glaciated Northeast, but did not recognize 
the coastal (tidal) freshwater wetlands. Stewart and Kantrud 
(1971) devised a system for classifying wetlands in the glaciated 
prairies. This system also relied on water and vegetation 
characteristics . 

In 1979, a wetland classification scheme was developed as a 
tool for managers to inventory wetlands and deep water habitats 
of the United States (Cowardin et al. 1979). This system defines 
wetlands by plants (hydrophytes) , soils (hydric soils) and 
frequency of flooding. The structure of the system is hierar- 
chical, progressing from system to class, subclass and dominance 
types. Modifiers for water regime, water chemistry and soils are 
applied at the class, subclass and dominance levels. Special 
modifiers are available for man or beaver modified wetlands. 
This scheme was used to classify the Swamp Gulch study area. 

3.3.2 Classification of the Swamp Gulch Wetland Study Site 

Based on the Cowardin et al. (1979) scheme, the Swamp Gulch 
site can be classified as a "Palustrine" system. The Palustrine 
system includes all nontidal wetlands dominated by shrubs, trees, 
perennial emergents, or emergent mosses or lichens. This 
category was developed to group the vegetated wetlands tradition- 
ally called marshes, swamps, bogs, fens, prairies or ponds. 
Palustrine systems may be situated almost anywhere on the 
landscape. 

In terms of the subclass, the majority of the Swamp Gulch 
site was classified as a "Persistent Emergent Wetland" because of 
the dominance by erect, rooted and persistent herbaceous 
hydrophytes. This subclass is typified by such plant types as 
rush ( Juncus ) , cattails (Typha) , sedge ( Carex ) , bulrush ( Scirpus ) 
and certain species of true grasses. The portion of the Swamp 
Gulch site dominated by shrubs was classified as a "Palustrine 
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Broad-leaved Deciduous Scrub - Shrub Wetland". This subclass is 
typified by alder f Alnus ) . willow (Salix) , and birch (Betula) • 
The western one-third of the site was classified as "Palustrine 
Needle-leaved Evergreen Forested Wetland" . This subclass is 
reserved for wetlands having trees taller than 6 meters. 

To completely describe the Swamp Gulch site (and any other 
wetland) certain modifiers at the class level need to be applied. 
These modifiers are for water regime, water chemistry (Eh and 
pH) , soil factors, and special modifiers for disturbed wetlands. 
For simplicity, we will limit our discussion of special modifiers 
to the dominant subclass present at the Swamp Gulch site, the 
'Palustrine Persistent Emergent Wetland'. The water regime for 
this subclass was 'nontidal, semipermanently or seasonably 
flooded ' . These descriptors cover the range from saturated 
surface soil to areas consisting of open water during the growing 
season. The salinity modifier was 'fresh' (i.e. less than 800 
umhos/cm) , the pH modifier was 'acid' (i.e. less than 5.5), and 
the soil modifier was 'organic'. The special modifier was 
'impounded', due to the activities of beaver. 

The full wetland classification for the majority of the 
Swamp Gulch site according to Cowardin et al. (1979) is, 
"Palustrine, Persistent, Emergent Wetland, nontidal, semiper- 
manent or seasonably flooded, fresh, acid, organic, impounded". 
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4.0 SURFACE HYDROLOGY OF A NATURAL WETLAND 



The Swamp Gulch surface wetland water hydrology was 
evaluated: 1) to determine the amount of AMD inflow and hence the 
calculated metal loading input to the wetland (see Section 7.0); 
2) to evaluate the general surface flow pattern to determine the 
area of impact and if channeling was evident; and 3) to evaluate 
the effect of AMD dilution by other unpolluted water sources 
flowing into the wetland. The following sections present these 
aspects of the study. 

4 . 1 METHODS 

The measurement of the AMD into the wetland was determined 
using a 22.9 cm (0.75 ft) throat Parshall type flume equipped 
with a Stevens Type F recorder. This unit was installed June 15, 
1987 in the upstream toe of the tailings dam on Swamp Gulch. The 
area beneath the flume was saturated and a mixture of tailings 
and bentonite was used to seal the area. The flume discharged 
directly into a wooden box culvert through the dam. The box 
culvert had sufficient cross-sectional area and slope so that no 
submersion of the flume outlet was observed. Complete records 
were obtained from June 15 through November 17, 1987, at which 
time the unit froze. The rate of flow was calculated using the 
equation: 

Q = 4WHa 1«522 yi 0.026 (Parshall 1950) 

where: Q = discharge, cubic feet per second (cfs) 
W = throat width (feet) 
Ha = upper gage head (feet) 

The rate calculations were integrated over time intervals varying 
from 0.25 to 24 hours as determined from the recorder charts for 
total daily flow. 

Precipitation was measured with a Belfort #5-780-6 universal 
recording rain gage installed near piezometer site E-1 on June 2- 
3, 1987 (Figure 4-1). The gage was mounted on 3.05 m (10 ft) 
long by 1.9 cm (0.75 in) inside diameter flanged pipes driven 
into the wetland approximately 2 m (7 ft) . The installation is 
equipped with an Alter (1937) type wind screen mounted at the 
orifice level, which is approximately 2.1m (7 ft) above the 
wetland surface. The gage was calibrated in 2.54 cm (1 in) 
increments from to 15.24 cm (0-6 in) following installation. 
Records are complete to date (Appendix B-2) , but some snowfall 
may have been missed from approximately December 15, 1987 through 
January 11, 1988, due to snow bridging. 
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Figure 4-1. Swamp Gulch Wetland Study site (with water sampling 
stations) . 



4.2 RESULTS AND DISCUSSION 

Inflow Of surface water to the wetland site is from two 
general sources: Swamp Gulch and the Blackfoot River. The 
Blackfoot River water is dispersed by several beaver dam systems 
below the confluence with Pass Creek (Figure 4-1) . Direct 
drainage from these dams flows in moderately well-defined 
channels southwestward past site AAA-3 into the study area. The 
Blackfoot River water drains into the beaver pond between sites 
C-3 and c-4 and then, via a well defined channel, flows past site 
C-3 to the pond directly east of site D-2. The flow at site C-3 
may approach 1 cfs at times and may average better than 0.25 cfs. 
Outflow from the pond east of site D-2 drains westward to a small 
pond adjacent to site E-2 and then through two ponds north of 
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site F-2. The pond immediately northeast of F-2 receives 
additional Blackf oot River water diverted via a channel near 
site E-3. Drainage from the pond northwest of F-2 flows 
approximately 37 m (120 ft) in a well defined channel into the 
Blackf oot River. This system likely determines the current 
southern extent of the Swamp Gulch acid drainage impact area. 
The pond northeast of site F-2 is also the likely point of 
accumulation of wetland drainage resulting from Swamp Gulch. The 
area topography generally confines this flow to the area 
immediately south of Highway 200 (Figure 4-2) . 

It should be noted that fish. Eastern Brook Trout, have been 
occasionally observed in this system from site AAA-3 to the pond 
near site E-2. These were observed frequently in the channel 
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from AAA-3 to the pond between C-3 and C-4, and only once at the 
pond near E-2 (during the April, 1987 sampling trip). 

The area surrounding the Swamp Gulch AMD discharge point is 
relatively uniform with no distinct channeling. The discharge is 
spread in a wide arc of roughly 180 degrees. 

The age of the present wetland drainage system is important 
m relation to mining and post-mining periods in determining the 
actual extent of the area potentially impacted from this source 
and if there have been any changes in the impacted area. A 
comparison of aerial photography from the years 1939, 1966, 1972 
and 1987 indicates that the present beaver dam system was in 
place before July 15, 1939, and that apparently mature Pinus 
contorta were established at that time in the same beaver dam 
areas where they are presently found. Several dead specimens 
were obtained and tree ring counts indicated ages from 22 to 70 
years (Table 4-1) . The oldest specimen exhibited intact bark and 
some needles suggesting it has been dead less than five years. 
Other samples had deeply checked surfaces with no bark with tree 
ring counts ranging from 22-63, which suggests these are some of 
the dead trees evident in the 1972 photographs. Based on these 
data, the minimum age for dam systems would be 78 years, having 
been established sometime prior to 1909. It is likely that the 
lodgepole pine could not have become well established on these 
dam areas until unsaturated surface soil conditions developed, 
possibly following abandonment. It is, therefore, probable that 
the system has remained essentially unchanged since &t least 1900 

Table 4-1. Beaver dam tree ring counts. 



Sample Ring 

Number Counts Condition Location 

1 70 Bark intact, some needles, 30 ft SW C-3 

12 inch diameter. 

2 22 No bark, deeply checked, 25 ft SW C-3 

3.5 inch diameter. 

3 47 No bark, deeply checked, 30 ft S D-2 

5.5 inch diameter. 

4 27 No bark, deeply checked, 25 ft W D-2 

7.5 inch diameter. 

5 63 No bark, deeply checked, 15 ft N C-4 

8.5 inch diameter. 
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and the mining impact area has not changed during this interval. 
The 1939 aerial photographs also indicate that no tailings dam 
was present adjacent to the highway at that time, which supports 
the apparent 1944 construction date for this feature. 

Flume records indicate flow rates varied from a minimum of 
0.1 L/s (0.004 cfs) on November 15, 1987, to a maximum of 14 L/s 
(0.5 cfs). The peak flow occurred July 17-18, 1987, in response 
to 6.6 cm (2.60 in) of precipitation received at the site during 
the two day period. Daily minimum and maximum flows were 12. 3m-^ 
(432 ft-*) and 734m3 (25,840 ft^), respectively. The average 
daily flow was 45m3 (1,588 ft^) . A comparison of flume hydro- 
graphs and precipitation gage charts suggests the time of 
concentration (the time interval from initial storm precipitation 
to the storm hydrograph peak) at the fltime is about ten hours 
under conditions of low antecedent precipitation and moderately 
large storms. The maximum intensity precipitation rate was 
recorded on July 22, 1987, when 1.4 cm (0.55 in) fell within one 
hour. The maximum daily total (24 hours) was 5.72 cm (2.25 in) 
which occurred July 17, 1987 (Appendix B-1) . This event was 
close to the reported ten year, 24 hour storm event of 5.6 to 6.1 
cm (2.2 to 2.4 in) (Miller et al. 1973). It is apparent that 
even the moderately intense storms do not produce peak runoff 
sufficient to cause any channeling in the wetland. Extreme 
events would be attenuated by the capacity of the highway culvert 
and highway overtopping would likely spread the flow to reduce 
wetland channeling. There is no evidence that this type of event 
has occurred in the last few years. Normal peak flow would be 
expected during spring snow melt and no data has yet been 
recorded for this time interval. 

Precipitation records from the site (Appendix B-2) suggest 
warm season precipitation is, in general, similar to that 
observed at the Lincoln Ranger Station with the exception of 
isolated storms such as the July 1987 event. The limited record 
is insufficient to determine if these isolated storms have a 
significant effect on the annual long term precipitation. The 
November 1987 through February 1988 records indicate consistently 
more snowfall occurs at the wetland than occurs at the Lincoln 
Ranger Station. This would be expected due to the higher 
elevation of the wetland. The insufficient record length at the 
present time prevents an accurate estimate of average annual 
precipitation, but it is apparent that the site does receive 
amounts in excess of the 48.1 cm (18.93 in) received at the 
Lincoln Ranger Station. A linear regression of monthly data from 
the two sites for the June 1987 through February 1988 records 
gives the following relationship. 

WP (cm) = (1.45) LRSP - .47 cm 

where: WP = wetland precipitation 

LRSP = Lincoln Ranger Station precipitation 

r2 =0.69 (r=0.83) 
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Using the regression equation with the average annual precipi- 
tation at the Lincoln Ranger Station gives a 69.3 cm (27.27 in) 
estimate of the annual precipitation at the wetland. This is 
somewhat less than the 73.66 (29 in) estimated from SCS (1981) 
isolines. A rough free water surface evaporation estimate as 
determined from Kohler et al. (1959) is approximately 66 cm (26 
in) . 

In summary, the Swamp Gulch AMD daily average flow is 
estimated to be 45m-» with a range of 12.3 to 734m3. The Swamp 
Gulch discharge into the wetland is generally well distributed 
without any notable localized channels, although there are 
preferred flow paths (see Section 5) . Blackfoot River and Pass 
Creek surface water is introduced into the wetland via beaver 
dams and two well-defined channels which have apparently been in 
place during the complete history of the mine. The channels have 
likely limited the extent of the AMD impact area and provide a 
sufficient unpolluted water volume to dilute the measured AMD 
volume to low levels. The Swamp Gulch wetland receives precipi- 
tation in excess of that known for the Lincoln Ranger Station. 
The precipitation is likely also slightly in excess of the 
estimated free water surface evaporation at the site and hence, 
no concentration of dissolved solids in wetland water likely 
occurs from evaporation. 
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5.0 HYDR06E0IOGY OF A NATURAL WETLAND 



The purpose of this part of the Swamp Gulch wetland study 
was to determine if residence time of mine-contaminated water in 
a natural wetland could be determined using a bromide tracer. 
Specific objectives of the study were as follows: 

1) Investigate the ability of bromide to serve as a tracer 
in the wetland. 

2) Identify water flow velocity in the wetland. 

3) Relate water flow velocity to wetland physical 
properties (hydraulic conductivity) . 

4) Determine residence time in the wetland area affected 
by AMD. 



5.1 LITERATURE REVIEW 

5.1.1 Hydraulic Principles of Tracer Use 

In hydrology, a tracer is either matter or energy carried by 
water which will give information concerning the direction and/or 
velocity of the water. When sufficient data are collected, 
tracers can assist in the determination of hydraulic conduc- 
tivity, porosity, dispersivity and other hydrogeologic parameters 
(Davis et al. 1980). Introducing a tracer at one point in the 
flow field and observing its arrival at other points, is the most 
direct method of determining groundwater velocity (Freeze and 
Cherry, 1979) . 

One factor which controls velocity is hydraulic conductivity 
which represents the ease with which water moves through the soil 
or aquifer. Hydraulic conductivity, "K", is often used to 
characterize an aquifer because it includes the properties of the 
fluid and the field of gravity as well as the properties of the 
porous medium such as permeability. Hydraulic conductivity has 
the dimensions of length/time (L/T) or velocity (Fetter 1980). 
Hydraulic conductivity, porosity and permeability values can vary 
widely in space and time. The relationship between velocity and 
hydraulic conductivity is derived from Darcy's Law and is shown 
in the following equation: 

V = (-K/ne) (dh/dl) (Eq. 5-1) 

Where: V = average linear velocity 
K = hydraulic conductivity 
ne = average effective porosity 
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dh/dl = average hydraulic gradient 
dh = change in head 
dl = change in length 

Interpretation of the results of tracer tests involves 
plotting the concentration of a tracer as a function of time or 
volume of water passing through the aquifer. In the resulting 
"breakthrough curve" the concentration is commonly given as a 
ratio of the measured concentration at the observation well, "C", 
to the initial tracer concentration injected, "Co". The average 
travel time of a non-reactive (conservative) tracer can be 
determined from a breakthrough curve for transport from the 
injection point to the observation point. The first arrival time 
of a tracer as it moves through the system represents the maximum 
velocity of the groundwater. The peak concentration of the 
tracer represents the average transit time of groundwater through 
the system, if a conservative tracer is used. Retardation of the 
transit time of a tracer is related to the breath of the 
breakthrough curve. Figure 5-1 gives hypothetical examples of 
breakthrough curves for a mixture of tracers injected as a single 
slug into an aquifer. Notice the change in the curve shape when 
the tracer is not conservative and interacts with the aquifer 
system. 




Breakthrough curves of four different tracer 
types: (a) is conservative, (b) some effect of 
sorption, (c) large effect of sorption, and (d) 
precipitated or destroyed (From Davis et al, 
1985) . 
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5.1.2 Tracer Flow Retardation and Tracer Dilution 

Tracers are not perfectly conservative and the concentration 
distribution of a water soluble substance which is transported in 
a porous medium by groundwater is affected by sorption, molecular 
diffusion and hydrodynamic dispersion (Gustafsson and Klockars 
1981) . Sorption includes adsorption and absorption processes. 
The many chemical processes which contribute to sorption result 
in retardation of tracer movement. Thus, velocity of the tracer 
is slower than that of 

the groundwater. Therefore in order to design a meaningful 
tracer experiment, the sorptive characteristics of the tracer 
must be known (Davis et al. 1985). Davis et al (1985) and 
Gustafsson and Klockars (1981) give equations to describe the 
effect of tracer sorption in relation to groundwater flow. 

Hydrodynamic dispersipn and molecular diffusion have the 
effect of diluting the concentrations of artificially induced 
tracers. Hydrodynamic dispersion generally affects short term 
tracer tests and molecular diffusion affects the concentration of 
slow moving tracers in heterogeneous materials (Davis et al. 
1985). Freeze and Cherry (1979) also state that molecular 
diffusion is important only at low flow velocities. Figure 5-2 
depicts movement by molecular diffusion. Note that no water 
movement is required for the dye to spread out (diffuse) in a 
direction tending to equalize concentrations in the blotter. 




spot of dye 

soaked blotter- 

(no water movement) 

initial conditions 




three hours 




one hour 



one day 



Figure 5-2. 



One dimensional example of movement by molecular 
diffusion (from Davis et al. 1985). 
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Hydrodynamic dispersion is the spreading which occurs both 
perpendicular to and in the direction of groundwater flow, of a 
water-soluble substance that is transported with the groundwater. 
The transit velocity of specific water-soluble substances may be 
higher or lower than the average groundwater velocity. Hydrody- 
namic dispersion is dependent on the velocity distribution in the 
medium and on molecular diffusion. Bear (1972) gives a quite 
detailed evaluation of hydrodynamic dispersion. Dispersion of a 
solute requires groundwater flow in a medium with a system of 
pores or channels. An example of hydrodynamic dispersion which 
is caused by unequal velocities of the ground water is shown in 
Figure 5-3. Tracer particles released at the same time and 
carried by the groundwater have different flow paths. This 
results in a more widespread distribution of tracer particles 
with time. 



General direction of water motion 
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Figure 5-3. 



One dimensional example of hydrodynamic dispersion 
for tracer particles A - G (from Davis et al. 
1985) . 
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5.1.3 Tracer Methods 

The initial step in conducting a tracer study is to collect 
as much hydrologic information as possible about the study site 
(Davis et al. 1985). This infonnation should include homogeneity 
of the aquifer, layers present, fracture patterns, porosity, flow 
system boundaries, hydraulic gradient and hydraulic conductivity 
(Davis et al. 1985). This hydrologic information is used to 
assess the groundwater flow direction and velocity. Direction 
and velocity are usually estimated by the use of Darcy's law or 
by performing a preliminary tracer test (Davis et al. 1985). The 
second step is to determine the best tracer to use for the 
conditions and objectives at the site. The third step in 
conducting a tracer study is to determine the correct amount of 
tracer to be used. This amount is based on dilution expected, 
natural background concentration and detection limit possible 
for the tracer. The next step is to determine the correct tracer 
method. 

One method involves the use of an environmental tracer which 
is a substance that exists in the soil before the investigation 
begins. This tracer can be artificial (man induced), semi- 
artificial (tritium) or natural (natural radioisotope) . A good 
environmental tracer must be free from chemical reactions, such 
as ion exchange and precipitation and must not react with the 
medium (Fried 1975) . Initial and boundary conditions must be 
known, chiefly by knowing the amount of tracer added and its 
history. The wetland system under investigation has the 
characteristic of a high organic matter content (large exchange 
capacity) . Trace metals (environmental tracers) commonly react 
with the organic material, making them unacceptable in this 
situation. 

Freeze and Cherry (1979) describe four main types of field 
dispersivity (tracer) tests. These are (1) single-well 
withdrawal -inject ion tests, (2) natural-gradient tracer tests, 
(3) two-well recirculating withdrawal -inject ion tests and (4) 
two-well pulse tests. In a natural gradient tracer study, the 
direction and velocity of the groundwater flow are very important 
(Davis et al. 1985). In the natural -gradient test, the tracer is 
introduced into the system and its migration is then monitored at 
one or more sampling points. Dispersivity values are obtained by 
fitting an analytical or numerical model to the experimental 
data. Davis et al. (1985) notes that: 

"It is not at all uncommon to inject a tracer in a well and 
not be able to intercept that tracer in a well just a few 
meters away, particularly if the tracer flows under the 
natural hydraulic gradient which is not disturbed by 
pumping." 
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A tracer may be injected as a slug or as a continuous source 
input. Mixing of the tracer during injection is important in 
most types of tests and can be as simple as pouring it in the 
water to be studied (Davis et al. 1985). For shallow wells a 
plunger can be surged up and down in the hole or the tracer can 
be released through a pipe with many perforations . The ideal 
condition is to inject the tracer into the water instantaneously 
as a slug. 

Taking field measurements of electrical conductance within, 
ahead of, and behind a tracer slug can minimize the number of 
field samples kept for laboratory analysis (Lee et al. 1980). 
Electrical conductance is the simplest and most inexpensive 
detection and analysis technique for ionic tracers and can be 
used as a break-through indicator (Slichter 1902) . In the 
wetland system being studied, the naturally high ionic content 
may make electrical conductance data difficult to interpret. 

5.1.4 Tracer Problems 

Many water tracing attempts are unsuccessful. Aley and 
Fletcher (1978) have found that the major causes of failure are: 
(a) insufficient hydrological field work before the tracer is 
injected, (b) tracing attempts during low flow conditions, and 
(c) failure to allocate sufficient time for tracing effort. 
Davis et al. (1980) found that tracer test failures are most 
commonly a result of incorrect choice of tracers, insufficient 
concentrations of tracers and a lack of understanding of the 
hydrogeologic system being tested. 

For a given head drop, expected travel time is a function of 
the distance squared, and therefore increases very rapidly with 
the distance. Davis et al. (1985) note that this relationship 
causes one of the most common errors in tracer tests, which is to 
conduct tests between points which are separated by too great a 
distance. Freeze and Cherry (1979) describe four main disad- 
vantages to the determination of groundwater velocity by the 
direct tracer method. 

1) Undesirably long periods of time, because of the fact 
that groundwater velocities are usually low, are 
normally required for tracers to move significant 
distances through flow systems. 

2) Numerous observation points are usually required to 
adequately monitor passage of the tracer through the 
study area because geological materials are typically 
quite heterogeneous. 

3) Because of (1), only a small and possibly nonrepresent- 
ative sample of the flow field is tested. 
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4) Because of (2), the flow field may be significantly 
distorted by the measuring devices. 

The concentration of ion to be injected should be well above 
the natural background concentration level found at the test site 
and high enough to ensure detectable levels in observation wells 
(Davis et al. 1985). The ion concentration injected must also be 
kept low enough so that density effects do not effect flow of the 
tracer. An insufficient quantity of tracer will result in an 
unsuccessful trace: too much tracer wastes materials and can 
degrade water quality. Dilutions of a tracer in transit from 
injection to sampling wells are almost always at least tenfold 
for "slug" injections and dilutions of ten thousandfold are 
common (Davis et al. 1980) . 

Lenda and Zuber (1970) describe a method to estimate the 
adequate amount of tracer and carrier to be injected. The 
effects of adsorption and the uncertainty from whether the 
observation well is exactly at the center of the tracer path 
should be accounted for by the use of a safety factor on the 
order of 10. Skibitzke and Robinson (1963) used tracers to show 
that solid particles (sand grains) retard diffusion in a porous 
medium. Biggar and Nielson (1962) concluded that the mere 
presence of a tracer downstream from the point of injection is a 
poor indicator of the velocity of the fluid. Pore geometry, 
water content changes and the magnitude of the interaction 
between the tracer and the porous medium are important in 
determining an accurate estimate of the fluid velocity. 

5.1.5 Bromide as a Tracer 

There is no such thing as the perfect tracer but Davis et 
al. (1980) note that the ideal ground-water tracer is nontoxic, 
inexpensive, moves with the water, is easy to detect in trace 
amounts, does not alter the natural direction of the flow of 
water, is chemically stable for a desired length of time, is not 
present in large amounts in the water being studied and is 
neither filtered nor sorbed by the solid medium through which the 
water moves. Different types of tracers include: water tempera- 
ture, solid particles, ionized substances (Br-) , stable isotopes, 
radioactive tracers, organic dyes, gases, and f luorocarbons . 
Davis et al. (1980) report that some of the most useful general 
tracers are bromide, chloride, rhodamine WT, and various f luoro- 
carbons. Most tracers have relatively limited or specialized 
uses. 

Tracer selection should be based on purpose of the study, 
type of aquifer system, aquifer characteristics, natural 
background concentration of the ions in the groundwater, and 
analytical techniques available (Davis et al. 1985) . In most 
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cases, anions are not affected by the aquifer medium (Davis et al 
1985) but the characteristics of some aquifers will cause 
retention or exclusion of anions moving through the system. 

Anionic tracers such as bromide (Br"") and chloride (Cl~) are 
particularly useful because of their low susceptibility to 
adsorption or ion exchange processes of natural aquifer material- 
s. Bromide does not appear to be lost by precipitation, 
adsorption or absorption and is biologically stable (Schmotzer et 
al. 1973) and therefore can be considered a conservative tracer. 
Bromide offers one of the best possibilities as a general tracer 
for groundwater studies (Davis et al. 1980). Most bromide 
compounds also have relatively low toxicities. Davis et al. 
(1980) note that bromide samples, being nonvolatile can be stored 
indefinitely without concern for tracer loss to the atmosphere 
and sampling can be done using inexpensive air-lift pumps. 

Davis et al. (1985) relate that bromide is perhaps the most 
commonly used ion tracer. These authors list advantageous 
characteristics of the use of bromide as a tracer as inexpensive, 
stable, low limit of detection, low background concentrations, 
low toxicity and no sorption. Bromide as a tracer is commonly 
injected as NaBr, CaBr2/ or KBr. 

The concentration of bromide in natural ground waters is 
roughly 1/300 that of chloride and usually <1 mg/L (Davis et al. 
1980; Vinogradov 1959). Detection of bromide is relatively 
simple with a specific ion electrode which has a lower limit of 
detection of about 0.4 mg/L. If natural water has 30 mg/L of 
chloride (suggesting the natural presence of 0.1 mg/L of Br~) and 
if the bromide tracer is introduced with a concentration of 1000 
mg/L, then a dilution factor of 10^ is possible before it is 
masked by the natural background (Davis et al. 1980). Concentra- 
tions in our introduced chloride tracer should not exceed about 
3000 mg/L because of increased density of the solution (Davis et 
al 1980) . Since the halides bromide and chloride behave 
similarly (Bohn et al. 1985) it is thought that bromide should 
not exceed the 3000 mg/L concentration. 

Davis et al. (1985) state that an advantage of anions, such 
as bromide, used as tracers is that they do not decompose and are 
not lost from the system. This statement is generally true but 
under certain circumstances anions such as bromide may be 
affected by anion exclusion and/ or anion exchange. As anions 
move through soil they do not come in contact with all of the 
soil water. This is termed anion exclusion and occurs in 
response to the fluid flow rate and the fact that water near 
negatively charged soil surfaces is relatively immobile. The 
result of this exclusion is that anions can move through the soil 
faster than one would predict on the basis of uniform association 
with all the soil water (Smith and Davis 1974) . The association 
between anion exclusion and cation exchange capacity is strong. 
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Anion exclusion is a manifestation of the unequal ion distribu- 
tion in the diffuse double layer surrounding charged colloid 
surfaces (Bohn et al. 1985) . Factors affecting anion repulsion 
(exclusion) include: 1) anion charge and concentration, 2) 
species of exchangeable cation, 3) pH, 4) presence of other 
anions, and 5) nature and charge of the colloid surface (Bohn et 
al. 1985). Thomas and Swaboda (1970) suggest that in soils with 
high cation exchange capacities anion exclusion causes anions to 
move much faster than they would if no interaction with clays 
were to occur. This suggestion agrees with the theory that anion 
exclusion is a function of negative charge. 

Grim (1968) mentions two types of anion exchange in clay 
minerals. One is the replacement of hydroxide ions and the other 
factor is related to geometry of the anion. The geometry of the 
bromide anion does not fit that of silica tetrahedral sheets, 
thus it cannot be so absorbed. Anion exchange would take place 
around the edges of the clay minerals, not on the basal plane 
surface. A factor which complicates anion exchange studies is 
that any free or exchangeable iron, aluminum or alkaline earth 
elements present in the clay may form insoluble salts with the 
anions (Grim 1968) . Mattson (1929) has shown that the adsorption 
of anions was found to be negative in neutral and alkaline 
solutions and that as pH decreases the capacity of clay minerals 
for holding anions increases. The interaction of cations with 
clays is much more frequent than with anions. Berg and Thomas 
(1959) found that sulfate and chloride anions are adsorbed in 
soils high in kaolin clays and aluminum and iron oxides. 
Chloride will desorb readily at pH values found in most field 
conditions, but at low pH values, chloride ions were not easily 
desorbed. They also found that sulfate is held much more 
tightly to these soil types than is chloride. 

The study site is in the headwaters area of the Blackfoot 
River and bromide toxicity was a concern. Alexander et al. 
(1981) studied the effect of sodium bromide on Fathead minnows. 
Sodium bromide has a low toxicity to Fathead minnows with the 
average LC50 (lethal concentration for 50% of the test popula- 
tion) for a 96 hour period being 16,479 mg/L. Tests by Barnes et 
al. (1981) show that the mean number of organisms (thirty-five 
species of algae and zooplankton were studied) increased with 
time due to the addition of 1000 ppb potassium bromide and 
nutrients to the system. Species diversity decreased slightly in 
this situation. Schmotzer et al. (1973) reported that bromide 
has a very low toxicity in humans at 50-100 mg of bromide/100 ml 
of blood. This translates to a human having to drink 12 liters 
of 200 mg/L bromide to be toxic. The oral LD50 (lethal dose) 
in rats is 3.5 g/kg (Merck Index 1983). 
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Martin et al. (1956) concluded that plant tolerance to 
bromide is quite variable. Carrot tops contained 2.5% (25,000 
ppm) bromide with no reduction in growth but citrus trees 
containing 0.2% (2000 ppm) demonstrated reduced plant growth. 

The desirable characteristics of using bromide as a soil 
water tracer were demonstrated by Onken et al. (1977) to be easy 
detection, unlikely contamination of the environment and lack of 
reaction with soil and soil constituents. Schmotzer et al. 
(1973) conducted a fairly extensive study of using bromide as a 
groundwater tracer and found bromide to successfully fulfill the 
requirements devised by Schmotzer to be as close as possible to 
being the ideal tracer. The favorable characteristics of bromide 
include low toxicity, high sensitivity of detection, little loss 
through precipitation, adsorption and absorption, high stability, 
low background concentrations, small sample size requirements 
(post sampling activation analysis) , low cost, government 
approval is relatively easily attained and bromide is biologi- 
cally stable. 

Vinogradov (1959) found that bromide correlated well with 
iodide content and both bromide and iodide content was propor- 
tional to the amount of organic material in the soils studied. 
The concentration of bromide is greater in humic soils and there 
IS practically no dependence of chloride content on soil organic 
matter content. The studies by Vinogradov (1959) suggest that 
bromide is sorbed by peats and that the amount sorbed decreases 
somewhat with aging of the peats. Vinogradov (1959) found that 
as soil organic carbon content increased the bromide and iodide 
content increased, but no effect was observed on the chloride 
content . 

Smith and Davis (1974) found that bromide is a good tracer 
for mimicking the movement of nitrate (NO3) through subsoils. 
Merrill et al. (1985) used bromide as KBr to trace NO3-N movement 
and indicate water flux in a study to develop an understanding of 
plant growth response to soil thickness over sodic minespoils. 
Onken et al. (1977) used sodium bromide and sodium nitrate to 
show that nitrate and bromide move together in the soil profile. 
They note that both nitrate and bromide are readily absorbed by 
plants but the rates of removal from soil are different. 

Tennyson and Settergren (1980) used a sodium bromide tracer 
to evaluate percolating water and ion movement in an irrigation 
saturated surface soil. Background levels of bromide in soil 
water, groundwater, and precipitation were measured and bromide 
movement was quantified by soil water sampling and post-sampling 
neutron activation analysis. They suggested that laboratory 
measured hydraulic conductivity was not adequate in evaluation of 
the site because the bromide tracer moved through the soil much 
more rapidly than the hydraulic conductivity suggested. Peak 
concentrations of the bromide tracer moved through 0.9 m of these 
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soils within 3.75 hours after field application. Tennyson and 
Settergren (1980) found indication of bromide retention occurring 
because bromide concentrations above background levels were 
present in the soils studied three weeks after application of the 
tracer. 

5.1.6 Wetland Hvdroloav 

Wetland systems can be considered quite heterogeneous. 
Figure 5-4 is shown to exemplify this point. Notice that not 
only is there seven orders of magnitude of hydraulic conductivity 
variation reported in the literature shown but there is often 
quite a large range of conductivity within a study site. Freeze 
and Cherry (1979) and Davis et al. (1985) note that a problem 
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Ranges of published field data on K (hydraulic 
conductivity) of peat: 1. Baden and Eggelsman 
(1961, 1963, 1964); 2. Eggelsman and Makela 
(1964); 3. Boelter (1965); 4. Ingram (1967); 5. 
Galvin and Hanranahan (1968); 6. Romanov (1968) 
7. Sturger (1968); 8. Dowling (1969) ; 9. Irwin 
(1970); 10. Yamamoto (1970); 11. Knight et al. 
(1971); 14. Paivanen (1973); 15. Galvin (1976); 
16. Dasberg and Neuman (1977) ; 17. Chason and 
Siegel (1986) . Adapted from Chason and Siegel 
(1986). 
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arises in determining the direction of water movement and travel 
time when local differences in hydraulic conductivity amount to 
several orders of magnitude. The following paragraphs describe 
some of the findings of specific authors in regard to wetland 
characteristics. 

Hydrogeologic factors (Winter and Carr 1980) that affect 
wetland ground water flow include: 

1) geometry of the geologic framework through which ground 
water flows, including flow boundaries; 

2) hydraulic conductivity of the geologic materials, 
including anisotropy — the ratio of vertical to 
horizontal hydraulic conductivity; 

3) recharge and discharge of the ground water system. 

In particular. Winter and Carr (1980) found three major inter- 
relationships between wetlands and ground water: 

1) some wetlands appear to recharge ground water; 

2) some wetlands are flow-through types where ground water 
enters one side and surface water seeps into the ground 
on the other side; 

3) some wetlands are discharge areas for ground water. 

These interrelationships are further complicated by changing 
throughout the year. 

One of the more important physical properties that affect 
the hydrologic features of bog areas is the hydrologic conduc- 
tivity of the peat soil horizons (Boelter 1965) . Dai and 
Sparling (1973) suggested using the piezometer method to 
determine the velocity of flow from the hydraulic conductivity. 
Ingram et al. (1974) used field experiments to find that 
estimates of hydraulic conductivity in humified peat increased 
with head, which is not in accordance with Darcy's law. It is 
suggested that Darcy's law may only apply in peat of low 
humidification. Boelter (1965) used field measurements of 
hydraulic conductivity (K) and found they covered a wide range of 
values. Chason and Siegel (1986) found that horizontal K was 
significantly greater than vertical K. 

Sturger (1968) studied the hydrologic properties of a 
mountain bog in Wyoming. Undecomposed surface peat had the 
lowest bulk density, 0.160 g/cc, and much of its volume consisted 
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of large voids which emptied at suctions of less than 0.10 bar. 
Greater suctions were required to drain the smaller pores of the 
higher bulk density, decomposed material found at depth. 

In 1969 Boelter found that the classification of peat 
materials based on degree of decomposition as measured by fiber 
content (> 0.1 mm) and bulk density can relate significant 
infoirmation about the hydraulic conductivity, water retention, 
and water yield coefficient. Regression equations were employed 
to determine a range in hydraulic conductivity and other physical 
properties for fibric, hemic and sapric peat materials (Table 5- 
1) . In studying a Minnesota spring fen-raised bog complex Chason 
and Siegel (1986) found that hydraulic conductivity, bulk density 
and humicity all vary widely through the peat profile and in most 
cases were not mutually dependent. There are two major hydro- 
logic zones within a peat soil column (Ingram 1983): 

1) the acrotelm, a thin upper aerated zone of fluctuating 
water conditions composed of undecomposed dead and live 
vegetation, 

2) The catotelm, an underlying anaerobic zone in a constant 
waterlogged state composed of more humified peat. 

Some authors (Dasberg and Neuman 1977; Ingram 1967; Sturger 
1968) note erosion channels or fractures at depths within wetland 
systems. The results of field and laboratory studies show that 
the properties of peat change drastically when it becomes 
partially desaturated (Dasberg and Neuman, 1977). Two layers 
result from this change, a permanently saturated layer below the 
zone of water table fluctuation, and an overlying partially 
unsaturated layer with properties that vary with depth. The 
saturated peat layer studied had a high porosity (90%) ; organic 
matter content (60% by weight) ; specific yield (20 to 30%) ; and 
a very high specific storage (0.7 x lO^^, to 1.7 x 10"^). The 
unsaturated peat layer was found to have a high bulk density, up 
to 0.8 g/cc near the surface. The hydraulic conductivity of the 
unsaturated layer is generally higher than 1.0 mm/h (3.0 x 10"^ 
cm/s) , but much more variable than the saturated peat layer. In 
summer, shallow cracks often develop which may increase the 
hydraulic conductivity by several orders of magnitude in the 
unsaturated surface layer. 

Ingram (1967) described a "water track" as a mire surface 
feature having a higher rate of water movement due to more 
steeply inclined water tables and/ or peat with higher hydraulic 
conductivity. These water tracks are associated with the 
presence of more eutrophic plant communities and are sites of 
greater ion supply. Erosion channels occur both on the surface 
and at depth in peat deposits with the associated vegetation 
often being unlike the surrounding vegetation. Subterranean 
erosion tunnels are frequently encountered in blanket bogs but 
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Table 5-1. Range of important physical characteristics of 
fibric, hemic, and sapric peat materials from 
northern Minnesota bogs (From Boelter 1969) . 



Bulk Total 0.1 bar Hydraulic 

Organic density porosity H2O content conductivity 

Material fa/cc) C%J f%^ fi o'^cm/sec) 

Fibric <.075 >90 <48 >180 

Hemic .075-. 195 85-90 48-70 2.1-180 

Sapric >.195 <85 >70 <2.1 



little is known of these concealed free drainage systems (Ingram 
1967) . Wetland hydraulic conductivity measured in situ by the 
piezometer method (Sturger 1968) at depths of 46 and 91 cm was 
23.9 X 10"3 and 16.1 X 10"^ cm per day respectively. Two of five 
piezometers at the 91 cm depth showed much higher hydraulic 
conductivity, indicating the piezometers terminated in or near a 
fissure. The fissures running through the deeper peat are filled 
with water under positive pressure and there are no surface 
physical features to indicate fissure location. Overall wetland 
hydrologic characteristics are found to be complicated by the 
fact that they vary greatly between wetlands, within wetlands and 
with time. 

Tracers have seldom been applied to wetland hydrology, 
"primarily because of the difficulty in choosing an appropriate 
tracer and the disturbance of the experimental area which results 
from frequent sampling (Girts 1986)." Knight et al. (1971) 
successfully used a tritixim tracer to reveal variable patterns of 
movement in a wetland under different drainage conditions. The 
flow rate variation was found to be attributed to factors within 
the peat rather than to analytical procedure. Girts (1986) 
successfully used bromide and chloride tracers to yield flow 
rates within the range of published values (see Figure 5-4) . 

Bowmer (1987) used bromide and dye tracers to study a man 
made wetland used to treat sewage. In this study theoretical 
detention time of effluent was predicted using the pore volume of 
the system. Preferential flow paths and "dead zones" were found 
using the tracers and a substantial proportion of the effluent 
was found to travel through the system faster than predicted by 
the theoretical retention time. The hydrology of the system was 
further complicated by occasionally high evapotranspiration rates 
and diurnal changes in water consumption. 
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5.2 MATERIALS AND METHODS 
5.2.1 Site Instrumentation 

Thirty-eight shallow small diameter wells (piezometers) were 
installed throughout the study area during the spring and summer 
of 1987. Twenty-eight of these wells are paired, with one well 
shallow and one deep resulting in 24 total well sites (Figure 4- 
1) . The shallow wells were completed to a total depth between 
1.2 to 1.5 m (4 to 5 ft) and deep wells to a total depth between 
1.8 and 4.9 m (6 and 16 ft). Wells were constructed of 2.5 cm (1 
m) diameter PVC pipe and Timco 0.025 cm (0.01 in) slotted, 
capped PVC screen. All wells were completed with a 20-30 grain 
sand pack and bentonite pellets. Figure 5-5 depicts a typical 
paired well site construction. All of these wells were developed 
by surging and overpumping. 

Well holes were constructed using 5.1 cm (2 in) diameter 
Giddmgs soil core barrels hand driven into the ground. A 7.6 cm 
(3 m) silage auger (Arts 1987) modified with a core cutter/hol- 
der was used to obtain soil samples in the to 1.2 m (0 to 4 ft) 
depth at each well site. An additional 32 auger holes were 
excavated near the highway culvert outflow to a depth ranging 
between .5 and .9 m (1.6 and 3 ft). These additional auger holes 
were positioned (Figure 5-6) so as to facilitate tracer sampling. 
A Parshall flume equipped with a Stevens Type F recorder was 
installed June 15, 1987 on Swamp Gulch below the acid mine seep 
to quantify flow rates into the wetland. A Belfort weighing 
precipitation gauge was also installed at the site in June 1987. 
The position of both the flume and precipitation gauge are noted 
m Figure 4-1. 

5.2.2 Acfuifer Characteristics 

Porosity, bulk density and particle density were determined 
on soil samples obtained during excavation of auger holes. Bulk 
density was determined using the Core Method described by Blake 
and Hartge (1986a). The Pycnometer Method was used to determine 
particle density (Blake and Hartge 1986b) . Total porosity, n, 
was calculated with the following equation (Danielson and 
Sutherland 1986) : 

n = (1 - pb/pp) , (Eq. 5-2) 

where pb = bulk density 

pp = particle density 

5.2.3 Hvdraulic Conductivity 

Saturated hydraulic conductivity was measured in the field 
using both auger hole and seepage tube (piezometer) methods. A 
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variable head method modified from the hydrostatic time-lag 
method of Hvorslev (1951), as reviewed by Cedergren (1977), was 
used to calculate saturated field hydraulic conductivity. This 
method is based on the rate at which water rises or falls in a 
hole after a known volume (slug) has been removed or added. In 
the case of the auger holes a bailer was used to remove a volume 
of water and the rate of water rise was measured using equipment 
designed similar to that described by Beers (1983) . For the 
piezometers a volume of a solid steel rod was added and the rate 
of water fall was measured using an electric tape. 
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Figure 5-5. Example of study site instrumentation and lithology. 
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Hydraulic conductivity was calculated for auger holes and 
piezometers, respectively, using the following equations: 



and 



K = R X fho - h i) (Eq. 5-3) 

16DS (t2 - ti) 

K = r^ In L lnjChi/h2l (Eq- 5-4) 

2L ( R ) t2 - ti 



Where: R = radius of cavity 

D = depth of cavity below static water table 

S = shape factor coefficient 

h^ = head at time one 

h2 = head at time two 

ti = time one 

t2 = time two 

r = radius of the well 

L = length of screen 

Hydraulic conductivity was converted to velocity (V) to 
obtain a preliminary estimate on which to base the tracer test 
sample timing. These velocities determined from hydraulic 
conductivity also facilitate comparison with velocities deter- 
mined by the tracer test. The formula used to make the K-V 
conversion is shown in equation 5-1. 

5.2.4 Tracer Experiment 

A preliminary dye tracer test was run using fluorescein dye. 
This experiment showed preferential flow direction, helped 
determine the tracer input point and resulted in an even better 
idea of flow rates near the tracer input point. In this 
experiment approximately 7 grams of fluorescein powder was added 
to 18.9 liters (5 gal) of water, mixed and introduced at 
different points along the mine drainage flow path. Tracer 
movement was monitored visually. 

Bromide was chosen as the tracer in the experiment for use 
in determining water flow velocities. Prior to adding bromide to 
the wetland, background (natural) bromide concentrations were 
measured. These background concentrations were measured using 
the bromide specific ion electrode as described below and by the 
EPA (1979) titrimetric Method 320.1. In the titration method the 
EPA (1979) found bromide recoveries from 83% to 99% and standard 
deviations from +.13 to +.44 for sample concentrations ranging 
from 0.3 to 20.3 mg/1. The titration method as used in this 
study showed 100% accuracy based on two laboratory spikes. 

The natural gradient tracer experiment utilized 208 liters 
(55 gal) of wetland water at a concentration of 1000 mg/L 
bromide. This water was added at the point of highway culvert 
outflow at 10:30 AM on September 9, 1987 (time zero). Fluor- 
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escein dye (18.9 liters - 5 gal) was added at ^his time as a 
visual aid in bromide sampling. The tracer was added to the 
system as an extended slug input at a rate of 0.39 1/s {0.bz:> 
gpm). 

Flow of the tracer plume was monitored by sampling the auger 
holes and wells and analyzing for bromide. Water sampling was 
done using portable peristaltic pumps. Care was taken in 
sampling so that a representative water sample was acquired, but 
not so much water that the flow field was greatly affected. This 
care usually amounted to pulling one well volume to ^inse the 
sample bottle three times and taking the 100 ml sample with the 
second well volume. 

Water samples were then analyzed for bromide in the field 
using a bromide selective ion electrode in conDunction with a 
pH/millivolt (mv) meter. Calibration curves (mv vs mg/L) were 
generated using standards measured under the same conditions as 
the samples (Orion Research 1982). Using the bromide electrode 
the concentration range is from 0.4 to 79,900 mg/L with a ±2% 
reproducability . 

The sampling intervals at the start of the experiment were 
nearly continuous near the input point. Sampling had decreased 
to once every other day, by September 14 (124 hours) . From 
September 26 (412 hours), to the end of the field experiment, 
November 21, 1987 (1755 hours), sample sets were taken once each 
week, water sampling was terminated after the tracer Pl^^^ had 
passed a particular well or auger hole. November water sampling 
was hampered by freezing conditions. The passing of a tracer 
plume was determined from the bromide breakthrough curves. 

Bromide breakthrough curves are plotted as bromide con- 
centration in milligrams per liter (mg/L) vs time m hours, 
concentrations were generated from the sample analysis data from 
a particular site and calibration curves from the time of sample 
analysis. Times represent the time a particular sample was taken 
from the wetland. The average velocity of water flow was taken 
to be the distance from the tracer input point to a particular 
well or auger hole divided by the time of peak concentration for 
that site breakthrough curve. 

5.3 RESULTS AND DISCUSSION 
5.3.1 Aquifer Character istics 

Field evaluation of auger hole and well cores showed a 
fairly distinct boundary between the upper, undecomposed (fibric) 
oraanic matter (sedge peat) and the underlying moderately 
decomposed (hemic) organic matter. Little distinction was made 
between moderately and well decomposed zones which are grey ro 
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black organic matter muck with occasional sand and/or gravel (<5 
cm) stringers and rotted wood pieces (<10 cm) . The horizoniza- 
tion scheme used distinguishes between two soil types 1) upper, 
undecomposed organic matter (acrotelm) , and 2) underlying muck 
(catotelm) (Figure 5-5) . The upper undecomposed zone is aerated 
and subject to fluctuating water conditions and water is 
transmitted comparatively faster in the upper zone compared to 
the underlying muck zone which is anaerobic and constantly 
waterlogged. In the wetland being studied, the acrotelm 
thickness ranges from 0.24 m (0.8 ft) to 1.04 m (3.4 ft) with an 
average of 0.58 m (1.9 ft) (Figure 5-7) and catotelm thickness 
ranges from 0.58 m (1.9 ft) to 4.30 m (14.1 ft) with an average 
of 2.16 m (7.1 ft) (Figure 5-8). This information is summarized 
in Table 5-2. Note from figures 5-7 and 5-8 that the greatest 
thickness for both zones are near the highway culvert outflow. 
These greater thicknesses may be due to sediment deposition from 
Swamp Gulch or metal precipitation originating from AMD water. 
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Figure 5-7. Acrotelm isopach map with 0.2 m (0.7 ft) contours. 

43 



CARBONATE MINE/> // 




A WELL SITE 

O SURFACE WATER SAMPLE 

SITE 
D SURVEY CONTROL STATION 
H FLUME 
O PRECIPITATION GAGE 



SEC20, TI5N,R6W 
LEWIS AND CLARK CO .MONTANA 

FEET 
200 



50 
METERS 



Figure 5-8. Catotelm isopach map with 0.5 m (1.6 ft) contours. 



Both sediment deposition and metal precipitation are 
indicated in the wetland soils data (Appendix C) . Sediment 
deposition m the acrotelm near the highway culvert outflow site, 
C-1 (0 to 91 cm), is indicated by a higher sand content (37%) and 
lower organic carbon content (12.4%) than the other three wetland 
acrotelm sites analyzed for these parameters. Acrotelm metal 
loading at this site is indicated by the highest observed iron 
content m the study site at 415,000 ppm (41%). The sediment 
deposition and metal loading at site C-1 has had an interesting 
effect on cation exchange capacity (CEC) and anion exchange 
capacity (AEC) . The low organic carbon content and high sand 
content of the acrotelm at site C-1 h^ve likely contributed to 
^tl^ °^f®^®^^i°^ ^EC of 4.9 meq/lOOg. : The resulting increase in 
CEC with depth at this site is the opposite of that expected and 
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observed at the other wetland sites measured for CEC. The AEC 
exhibits an inverse relationship with CEC and decreases with 
depth at site C-1. The high AEC (90 meq/lOOg) at site C-1 may 
have the effect of retarding the bromide anion flow velocity by 
sorption of bromide on the exchange sites. 

Values for bulk density, particle density and porosity were 
measured on three surface soil samples (acrotelm, to 30 cm 
depth) . Undisturbed samples were not possible below the 30 cm 
depth due to the sampling methods used and the nature of the 
wetland soils. The average bulk density (0.26 g/cc) is higher 
than that found by other studies at this depth in a wetland 
(Sturger, 1968 - 0.184 g/cc; Chason and Siegel, 1986 - 0.06 to 
0.165 g/cc and; Boelter, 1969 - <.075 g/cc). The distinctness 
between the measured bulk density and that from other studies may 
be due to differing vegetation types and/or soil compaction 
during the sampling procedure. 

Particle density values range between 2.077 and 2.942 g/cc 
with an average of 2.595 g/cc. This average density is less than 
the 2.65 g/cc used as an approximation for many mineral soils 
(Danielson and Sutherland 1986) , but higher than the usual case 
for humus which is <1.5 g/cc (Blake and Hartge 1986b). The 
relatively intermediate, average particle density is presumed to 
be the effect of the high organic matter content and high heavy 
metal loading in the wetland. 

Bulk density and particle density were used to find the 
total porosity (equation 5-2). Total porosity ranged from 87.35% 
to 94.27% with an average of 90.30%. This average porosity is 
similar to that found by Boelter (1969) which was >90% (Table 5- 
1), but higher than that calculated by Girts (1986), 66.1% to 
76.1%. The values found by Girts (1986) were for partially and 
greatly decomposed peat and assumed no water was held within 
plant tissues. Boelter (1969) noted that total porosity 
decreased with degree of decomposition or depth in a wetland. 

Effective porosity is less than total porosity because in 
small pores the retention forces are greater than the weight of 
water. Effective porosity in practice may be considered equal to 
the specific yield of an unconfined aquifer (Kruseman and Ridder 
1983). Generally, specific yield is the water yielded by gravity 
drainage from water bearing material (Lohman 1979) . Boelter 
(1965) found specific yields in undecomposed peat ranged between 
0.79 and 0.33 cc/cc and that moderately to well decomposed peat 
ranged from 0.22 to 0.10 cc/cc. The averages of the above 
ranges, 0.54 and 0.186 cc/cc respectively, were used as the 
effective porosity for converting field hydraulic conductivities 
to velocities. The decrease in effective porosity with decom- 
position is appropriate because the more decomposed peat 
materials have more small pores which are not easily drained. 
The screens of the deeper wells were completed within or near the 
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gravel zone lying below the wetland. The effective porosity of 
deep wells completed in gravel, as used in conductivity to 
velocity conversions, was 0.25 cc/cc (Driscoll 1986). 

Results of water level measurements in both wells and auger 
holes (Appendix B-4) indicate both upward and downward vertical 
hydraulic gradients in parts of the wetland. The artesian 
gradient (upward vertical) can be explained by two different 
arguments: 1) a confined aquifer system at depth, or 2) the 
wetland serves as a groundwater discharge area. Hydraulic 
conductivities which will be discussed later, suggest that the 
wetland system is capable of transmitting water throughout the 
soil profile (not confined). Therefore, the wetland is con- 
sidered both a groundwater discharge and recharge area. Field 
observations, such as at the culvert outflow, suggested water 
flowing into the wetland also seeped into the soil profile and 
became groundwater. Surface flow channels which would disappear 
into the ground and apparently to reappeared down gradient were 
also noticed in the wetland. Situations where wetlands function 
as groundwater discharge, recharge and flow-through areas are not 
uncommon (Winter and Carr 1980) . 

Deep and shallow piezometric maps and a potentiometric map 
are shown in Figures 5-9, 5-10, and 5-11, respectively. A 
general comparison of these maps reveals that the hydraulic 
gradient does not change between water surfaces measured, and 
there is a hydrologic and topographic (Figure 4-1) high near the 
culvert outlet. The average hydraulic gradient across the 
wetland where the hydraulic conductivity tests were performed is 
0.025 m/m. The change in head within a shallow and deep well and 
auger hole nest is generally less than 0.30 m (1 ft). This 
relationship can be seen in Figure 5-5. Overall piezometric 
surface levels stayed constant throughout the period measured. 
The potentiometric surface fluctuated approximately 0.30 m for 
short periods of time (days) during the tracer test. These 
fluctuations may be attributed to increased water use by 
phreatophytic plants (Mayboom 1967) during period of high air 
temperatures and the lack of evapotranspiration during the 
winter. This is further evidenced by slight diurnal fluctuations 
in water flow rate at the Swamp Gulch flume. The rate of water 
flow into the wetland from the culvert exhibited little change 
during the tracer study as indicated by the flume recorder charts 
(Appendix B-1) . Precipitation during the tracer study was very 
low (Appendix B-2) . 

5.3.2 Hydraulic Conductivity 

Hydraulic conductivity tests were run on auger holes and 
wells in the area of expected metal loading and expected tracer 
movement. This area was determined from water chemistry data and 
hydraulic gradients in the wetland. Field determined hydraulic 
conductivity values and the calculated water velocities for each 
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Figure 5-9 



Deep piezometric map on 7/27/87 with 0.5 m 
contours . 



auger hole and well tested are given in Table 5-3. 

When field hydraulic conductivity values were grouped, some 
differences were found between auger holes and wells (Table 5-4) . 
Lowest conductivity values were found in the shallow wells with 
their screened intervals within the catotelm (well decomposed 
organic matter, muck) . Highest conductivity values were found in 
the deep wells with their screened intervals within both the 
catotelm and the underlying gravel. Auger holes which have their 
open intervals within the acrotelm (undecomposed organic matter) 
had conductivities intermediate in comparison to the wells (Table 
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Figure 5-10. 



Shallow piezometric map on 7/27/87 with 0.5 m 
contours . 



5-4) . These values fall within the range of conductivities found 
in other wetlands studied (Figure 5-4). The range of hydraulic 
conductivities from the wetland, 1 x 10"^ to 1.5 x 10"^ cm/s, 
compares to the range commonly found for silty sand (Freeze and 
Cherry 1979) . The decrease in conductivity with depth, 
excluding the deep wells completed in gravel, is a common 
occurrence when moving from the acrotelm to the catotelm (Chason 
and Siegel 1986) . 

An argument to account for the small conductivity difference 
between auger holes and wells is that the tests measured part of 
the same zones. In viewing the depths of auger holes tested 
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Table 5-3. 


Field 


hydraulic 


conductivit 


:y (K), effec 


;tive porosity (ne) , 




average hydraulic gradient 


(dh/dl) and 


calculated average 




velocity (V) for each site. 










Site 


K 1 


'cm/s) n« 


rcc/cc) dh/dl (m/m) 


V( 


cm/s) 


B-2-AH 


1.06 


X 


10-3 


.540 


.025 


4.91 


X 


10-5 


BC-1-2-AH 


5.60 


X 


10-4 


.540 


.025 


2.59 


X 


10-5 

r- 


C-l-AH 


7.02 


X 


10-4 


.540 


.025 


3.25 


X 


10"5 


C-1.5-AH 


1.06 


X 


10-3 


.540 


.025 


4.91 


X 


10-5 


C-2-AH 


9.53 


X 


10-4 


.540 


.025 


4.41 


X 


10-5 


CD-1-2-AH 


1.18 


X 


10-3 


.540 


.025 


5.46 


X 


10-5 


CD-l-AH 


4.54 


X 


10-4 


.540 


.025 


2.10 


X 


10-5 


D-l-AH 


4.18 


X 


10-4 


.540 


.025 


1.94 


X 


10-5 


E-l-AH 


3.40 


X 


10-4 


.540 


.025 


1.57 


X 


10"5 
10-5 
10-5 
10-3 

r- 


E-2-AH 


1.08 


X 


10"3 


.540 


.025 


5.00 


X 


F-l-AH 


1.11 


X 


10-3 


.540 


.025 


5.14 


X 


B-2-D 


1.01 


X 


10"2 


.250 


.025 


1.01 


X 


B-2-S 


1.71 


X 


10-4 


.186 


.025 


2.30 


X 


10-5 
10"^ 
10-5 


BC-1-2 


3.89 


X 


10-5 


.186 


.025 


5.20 


X 


C-l-D 


2.81 


X 


10-4 


.186 


.025 


3.78 


X 


C-l-S 


6.02 


X 


10-4 


.186 


.025 


8.09 


X 


10-5 
10-5 
10-5 
10-5 

rr 


C-1.5 


3.00 


X 


10-4 


.186 


.025 


4.03 


X 


C-2-D 


1.22 


X 


10-4 


.250 


.025 


1.22 


X 


C-2-S 


1.44 


X 


10-4 


.186 


.025 


1.94 


X 


CD-1-2 


2.22 


X 


10-4 


.186 


.025 


2.98 


X 


10-5 


CD-I 


1.40 


X 


10-4 


.186 


.025 


1.88 


X 


10-5 
10-4 


D-l-D 


1.69 


X 


10-3 


.250 


.025 


1.69 


X 


D-l-S 


6.95 


X 


10-5 


.186 


.025 


9.30 


X 


o in in L 
o o o 

H H 1-1 


D-2-D 


6.62 


X 


10-4 


.186 


.025 


8.90 


X 


D-2-S 


1.92 


X 


10-4 


.186 


.025 


2.58 


X 


E-1 


9.16 


X 


10-5 


.186 


.025 


1.23 


X 


10-= 


E-2 


3.21 


X 


10-5 


.186 


.025 


4.30 


X 


10"° 
10-5 


F-l-D 


5.05 


X 


10-4 


.186 


.025 


6.79 


X 


F-l-S 


1.53 


X 


10-5 


.186 


.025 


2.10 


X 


10-^ 



that the auger holes also measured the decomposed organic matter 
zone (catotelm) layer thereby possibly reducing conductivity. A 
similar argument is true of the shallow wells some of which had 
the top of the screened interval and sand pack very near or 
within the acrotelm. The effect of this situation is that well 
tests measured a portion of the acrotelm in addition to the 
catotelm and may have greatly overestimated the actual conduc- 
tivity of the catotelm. The overall result is an average 
conductivity from the water table to shallow well total depths, 
but not the two distinct conductivity differences that might be 
expected. The average conductivity of both the auger holes and 
shallow wells is 4.9 x 10*4 cm/s and the velocity is 3.0 x 10"^ 
cm/s. 
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Table 5-4. Average hydraulic conductivity (K) from field 
results and calculated average velocity (V) . 



Site Type K fcm/s) V (cm/s) 

Auger Holes 8.11 x 10-4 3.75 x lO'S 

Shallow wells 1.68 X lO'^ 2.26 X 10 ^ 

Deep wells 2.22 x 10*3 2.31 x 10 ^ 

overall Mean 8.38 x lO'^ 7.14 x 10 => 



The hydraulic conductivity determined in this study is the 
effect of both vertical and horizontal conductivity. No 
distinction between vertical and horizontal conductivity was 
possible with the methods used in the field. It is presumed from 
other studies (Girts 1986; Chason and Siegel 1986; and Dai and 
Sparling 1973) that horizontal conductivity would be greater than 
vertical conductivity in the sedge peat wetland being studied. 

Average velocities calculated using field hydraulic 
conductivities and equation 5-1 are shown in Tables 5-3 and 5-4. 
Velocities determined were lower than the associated field 
hydraulic conductivity. This order of magnitude change (Table 5- 
4) is explained, in part, by the effective porosity chosen for 
the velocity calculation. From equation 5-1 it can be seen that 
as effective porosity decreases the water flow velocity will 
increase if the hydraulic conductivity and gradient remain 
constant. Another explanation for the lower velocities deter- 
mined for the field sites measured is the small average hydraulic 
gradient in the wetland (0.025). The hydraulic gradient of the 
potent iometric surface in the area of tracer input, culvert 
outflow to C-1 is 0.058. Incorporating this increased gradient 
into velocity calculations increases the average velocity from 
auger hole sites to 8.7 x lO'^ cm/s and velocity from shallow 
well sites to 5.2 x lO'^ cm/s. 

5.3.3 Tracer Experiment 

Water flow velocity was also determined in a bromide tracer 
experiment. Bromide was chosen as the primary tracer in this 
study for a number of reasons. Many authors have successfully 
used bromide as a tracer in a variety of environmental conditions 
(Onken et al. 1977; Smith and Davis 1974; Merrill et al. 1985; 
and Tennyson and Settergren 1980; to name a few). In addition. 
Girts (1986) and Bowmer (1987) successfully used bromide to trace 
wetland water systems. Davis et al. (1985) note that bromide is 
perhaps the most commonly used ion tracer. Schmotzer et al. 
(1973) consider bromide to be as close as possible to the ideal 
tracer. Bromide is a good tracer choice, relative to other 
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tracers, in the wetland system studied because it is inexpensive, 
stable, has a low limit of field detection with a specific ion 
electrode (0.4 mg/L) , has low background concentrations, exhibits 
low toxicity and has low sorption. The relatively few disad- 
vantages of bromide as a tracer are anion exclusion, anion 
exchange and the presence of other ions which interfere with 
operation of the electrode. Data from this study indicate 
bromide is a good choice in the system studied. 

Bromide background concentrations and the concentrations of 
ions which interfere with the electrode were analyzed before the 
tracer was introduced to the wetland. Background concentrations 
were analyzed by both the specific ion electrode and by titration 
(Table 5-5) . No auger hole background concentrations were above 
the 0.4 mg/L limit of detection using the bromide specific 
electrode. Well water which had background values above the 
electrode lower limit of detection, 0.4 mg/L, had that value 
subtracted before plotting the breakthrough curve. The 
subtraction of background values determined from the bromide 
electrode was necessary only for five wells. The average bromide 
value subtracted was 0.47 mg/L with a range from 0.43 to 0.55 
mg/L. This subtraction was found to change the shape of the 
breakthrough curves, but not the placement of the concentration 



Table 5-5. 



Bromide background concentrations in the tracer 
study area of the Swamp Gulch wetland. 



Site 



Bromide Concentr ation rma/L) 

Determined by Ion Determined by 

Electrode 9/8 - 9/87 Titration 7/28 - 31/87 



AH-6 

AH-14 

AH-17 

AH-22 

AH-29 

AH-32 

BC-I-2I 

BC-l-S 

C-l-S^ 

C-1.5^ 

C-2-S 

CD-l^ 

CD-I-2I 

D-l-S 

D-2-S 

BLANK 

BC-1-SURF 



0.40 
0.40 
0.40 
0.40 
0.40 
0.40 
0.43 
0.40 
0.44 
0.50 
0.40 
0.45 
0.55 
0.82 
0.40 
0.40 



< 





.10 


< 





.10 


< 





.10 




0. 


10 


< 





.10 




0. 


20 


< 





.10 




0. 


20 



^Backgrounds values subtracted from breakthrough curve. 
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peak. Bromide background levels measured by titration were never 
larger than 0.2 mg/L and 80% of all the wetland samples measured 
were <0.1 mg/L (detection limit). Review of the two methods of 
analysis suggests that the titrimetric method may be more 
accurate than the field electrode method, especially at low 
bromide concentrations. The advantages of using the field 
millivolt meter and bromide electrode are quick, inexpensive 
results at the site. This allows for rapid on-site sampling 
decisions and no need to save a large number of samples for 
expensive laboratory analysis. 

Ions which interfere with the bromide specific electrode, in 
order of increasing interference, include S""^, CN", I", NH3, Cl~ 
and OH" (Orion Research 1982). Sr''"^ (strontium) must be absent 
(Cole-Parmer 1987). Water quality analysis (Appendix B) of the 
acidic wetland system indicate that, of the interfering ions 
above, only chloride and ^trontium may cause inaccurate electrode 
readings. Williams (1979) relates that a chloride concentration 
of 20 times the bromide concentration will cause error in the 
bromide determination. At the lower limit of detection with the 
bromide specific electrode (0.4 mg/L) a chloride concentration of 
8 mg/L can cause error. Water analysis at well D-l-S and surface 
water near BC-1 showed chloride concentrations of 16 mg/L and 8 
mg/L respectively on 7/29/87 (Appendix B) . The most likely 
source of this chloride is highway road salt runoff. All wetland 
sites sampling strontium showed less than the lower limit of 
detection (0.1 mg/L) except for well F-l-S which showed 0.1 mg/L. 
Chloride and strontium may thus have created limited variation in 
the readings for bromide during the tracer test. 

The time involved in a uniform (natural) flow tracer method 
can be quite long and the method is usually limited to short or 
intermediate distances not more than 100 m (Davis et al. 1985). 
The velocities determined from the field hydraulic conductivities 
and preliminary tracer test suggested that the piezometer network 
might not be close enough to the tracer input point. In addition 
to variable velocities, conductivity was noticed to decrease with 
depth and this suggested that flow would be predominantly in the 
acrotelm. Therefore, a grid of shallow sampling sites was 
installed nearer the tracer input point. Inexpensive auger holes 
were used as the sample sites. The preliminary dye tracer test 
exposed preferential water flow paths and auger holes 27 to 32 
were placed along these paths (see Figure 5-6). 

Bromide standards ranging from 0.1 to 1000 mg/L were run 
during each sample session to determine the corresponding mv 
value for that time period. An example of the resulting 
calibration curve is shown in Figure 5-12. Regression equations 
were determined for each calibration curve. The coefficient of 
determination, (r^) , for each regression was maintained at 0.99. 
Values which were below the 0.4 mg/L limit of detection were 
excluded from the regression. Field mv values for each sample 
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Figure 5-12 



Typical calibration curve generated for the tracer 
study . 



site were converted to mg/L values using the regression equation 
for that corresponding time period. Distilled water blanks were 
run during each sample session and all were less than the 0.4 
mg/L limit of detection. 

Bromide breakthrough curves were developed using actual 
bromide concentrations on the y axis instead of the commonly used 
concentration/concentration injected (c/co) . This was done to 
facilitate calculations and has no effect on the shape of the 
breakthrough curve. Changing the y axis to concentration/ concen- 
tration injected is commonly done in tracer studies to facilitate 
comparison of different types of tracers introduced at the same 
time. An example of a bromide breakthrough curve generated is 
presented for site AH-14 in Figure 5-13. Based on the work by 
Davis et al. (1985) (Figure 5-1), the shape of this curve (steep 
slopes) indicates little sorption of tracer in the porous medium 
between the input point and sampling site. The curve shape also 
exemplifies the expected result that bromide can be considered a 
fairly conservative tracer in the wetland studied. 

There are many possible explanations for shapes of the 
breakthrough curves. Some of these explanations have been 
discussed earlier such as retardation by sorption and dilution by 
dispersion and diffusion. Site AH-25 generated data which 
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Figure 5-13. Tracer study breakthrough curve for site AH-14, 6.1 
m from tracer input showing little effect of 
sorption. 



resulted in only one point above the lower limit of detection. 
All the neighboring sites, 1.5 to 3.7 m (5 to 12 ft) away 
generated many bromide breakthrough curve points up to 4.8 mg/L. 
This suggests zones of very low flow which are described as "dead 
zones" by Bowmer (1987). At some sample sites the sampling 
interval was not frequent enough to determine the exact break- 
through peak. In the case of site AH-9 (Figure 5-14) the 
resultant curve is flat topped between hours 169 and 220. Both 
of these times were used to calculate velocity in addition to the 
peaks at 30 and 54 hours. Subtraction of the background values 
resulted in no detectable tracer peaks on the breakthrough curve 
for some of the well sampling sites such as C-l-D, D-l-S, and D- 
2-S. 

Placement of the breakthrough curve peak is important in 
that the associated time represents the passing of the tracer 
plume. This time is also used for determination of average 
velocity. Widely spaced sampling intervals may cause the exact 
breakthrough curve peak to not be found and the velocities 
determined may be too large or too small. Bromide values below 
the 0.4 mg/L limit of detection were not considered as possible 
peaks on the breakthrough curves. A common observation of the 
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Figure 5-14. Example of infrequent sampling effect on break- 
through curve for site AH-9, 6.1 m from tracer 
input . 

breakthrough curves shows more than one peak of tracer concentra- 
tion. An example of such a curve is given in Figure 5-15 for 
site AH-19. This curve seems to indicate two bromide tracer 
plumes moving past the site with the second plume being more 
effected by sorption. 

How breakthrough curves acquire more than one concentration 
peak is not specifically discussed in the literature reviewed. 
Possible explanations of multiple breakthrough peaks could 
include: inadequate tracer mixing during input, a heterogenous 
medium with layers of different flow velocity, periodic flushing 
and/ or cross contamination of samples. The inadequate mixing 
theory is negated by the breakthrough curve shape of the first 
sampling site encountered (Figure 5-16). This site, AH-3, which 
showed no multiple peaks and exhibited a high concentration and 
short time interval to the breakthrough curve peak for this short 
distance from tracer input. A heterogenous medium is an entirely 
plausible explanation as layering was indicated during soil 
sampling (Table 5-2) and preferential flow paths were found in 
the wetland. These situations are also found in other wetlands 
(Ingram 1967; Dasberg and Neuman 1977; Sturges 1968). The 
possibility of samples being drawn from more than one zone of 
flow is good because many sampling sites have their screened 
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Figure 5-15. Breakthrough curve representing two tracer plumes at 
site AH-19, 9.1 m from tracer input. 
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Figure 5-16. Steep breakthrough curve indicating adequate tracer 
mixing at site AH-3, 3 m from tracer input. 
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interval in contact with more than one soil horizon (Figure 5-5 
and Table 5-2). It is theorized that pulse flow could cause 
multiple breakthrough peaks by flushing higher concentration 
water which was previously held in pore spaces during different 
flow rates. Periodic flushing is indicated by mildly fluctuating 
water levels during the study period (Appendix B-4) . The diurnal 
variation noticed in flume charts could also contribute to this 
type of flow. 

The possibility of cross contamination of samples, analysis 
error and/or some other means of acquiring faulty values is one 
which must be considered in any experiment. Blank samples 
(distilled water) were run through the same analysis procedures 
as the wetland samples. The bromide concentration in blanks 
never was greater than the 0.4 mg/1 limit of detection and cross 
contamination was not a problem. Another explanation for 
multiple breakthrough curve peaks is error during analysis for 
bromide in the field. Analysis methods were constant throughout 
the study and should not have created the noted variation. The 
final explanation for multiple peaks is a periodic contamination 
by the interfering ions, chloride and strontium which seems 
unlikely. 

Given these possible explanations for multiple tracer peaks 
it seems that they are true bromide values and have been 
processed accordingly. The previously shown breakthrough curve 
for site AH-14 (Figure 5-13) yields a velocity of 1.4 x 10"-* cm/s 
(6.1 m/121.7 hr) . Table 5-6 shows the velocities associated with 
the bromide breakthrough curve peak/peaks chosen for each site. 
The velocities shown represent the rate at which a tracer plume 
passed that particular site. Site locations are shown in Figure 
5-6. The overall mean velocity represented by all of the sites 
and peaks is 2.6 x IC^ cm/s with a range of 3.1 x lO"-*- to 6.2 x 
10-4. 

The velocities were compared in three different ways 1) 
between auger holes and wells, 2) between the highest peak of the 
breakthrough curve and the remaining peaks, and 3) between sample 
sites within the grid area, <18.3 m, and sites outside the grid 
area, 27 to 38 m (Figure 5-6) . These results are shown in Table 
5-7. Several observations were apparent: 1) there is little 
difference between tracer velocities determined at auger holes 
and wells, 2) the velocities determined from other than the 
highest peak of the breakthrough curve are notably slower than 
from the remaining peaks, and 3) there is little difference (same 
order of magnitude) in velocities determined near the tracer 
input and those velocities from farther away. 

The apparently similar velocities observed at auger holes 
and wells can be explained by the fact that portions of the 
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Table 5-7. Tracer velocity comparison. 



Comparison Number Mean 

Type of Cases Velocity rcm/s) 

Auger Holes 83 2.9 x 10"2 

Wells 18 1.4 X 10~2 

High Peaks 45 5.0 x IC^ 

Remaining Peaks 56 6.4 x 10"-^ 

Grid Area 70 3.1 x 10"2 

Outside Grid 31 1.4 x lO"'^ 



-2 



Overall 101 2.6 x 10 



screened interval of most of the wells are within the same layers 
of conducting material (Table 5-2 and Figure 5-6) as the auger 
holes. The total depth of the shallow wells has little effect in 
that they are able to pick up tracer water from upper layers 
(acrotelm) . This "pick up" can occur through the well sand pack 
which is usually 30 cm (1 ft) from ground surface, below the 
bentonite seal. The auger holes are subject to the acrotelm 
flowing water slightly sooner than the shallow wells as indicated 
by the slightly faster velocity. The comparison of the high peak 
and remaining peaks supports the theory of residual tracer water 
being flushed through after the main tracer plume. These 
remaining tracer peaks could represent tracer water flowing 
through the less conductive (catatolem) peat layers hence 
reducing flow velocity. The higher velocities determined in the 
auger hole grid area, although within the same order of magni- 
tude, may reflect the steeper gradient in the area near the 
culvert outflow. Considering all the velocities from the 
different comparisons it is concluded that the overall mean of 
2.6 X 10"2 cm/s is the most representative of the wetland system 
measured . 

The presence of preferential flow paths is demonstrated by 
results of the bromide tracer experiment. Breakthrough curves 
peak at just 1.4 hours from tracer input for sites AH-28 and AH- 
31, both 15.2 m (50 ft) from tracer input (Table 5-6). This 
results in a tracer plume of the shape represented in Figure 5- 
17 . These preferential flow paths indicate the tracer is moving 
in both a south-southwesterly and easterly direction. The 
bromide concentration at site AH-31 is more than twice as high 
(52 mg/L) as site AH-28 (20.5 mg/L) at the same time (1.4 hrs) 
and indicates that more tracer is flowing in the south-south- 
westerly direction. This southwesterly skewed radial pattern is 
expected as the culvert outflow is a topographic high point and 
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the regional hydraulic gradient and topographic slope is to the 
southwest. Preferential flow paths are found in man made 
wetlands (Bowmer 1987) and in natural systems are sometimes 
called water tracks (Ingram 1967) . Such paths may cause a 
substantial proportion of the water (tracer and/or AMD polluted) 
to travel through the system faster than predicted by the field 
hydraulic conductivities. Brooks (1984) recommends man made 
wetlands be designed to avoid channelization (preferential flow 
paths) and encourage low velocity sheet flow by use of flow 
obstructions. Burris et al. (1984) shows that as flow rate 
decreases effective iron removal from AMD increases in a wetland. 

Tracer flow velocities were faster than flow velocities 
estimated from field hydraulic conductivity tests on auger holes 
and wells. The relationship between tracer determined velocity 
and field hydraulic conductivity determined velocity is compli- 
cated by a number of factors. The measurement of the conduc- 
tivities and the tracer test were done at different times. In 
effect, the measures are of different systems, but in practice, 
are assumed the same in view of the minimal change noticed in the 
system over an extended period. The field conductivity deter- 
mined velocities are more representative of the wetland at depth. 
The tracer velocities, which include preferential flow paths, are 
more representative of the true movement of AMD waters. The 
amount of area and time represented by the tracer test is large 
compared to the small area and time that the field conductivities 
represent. This infers that the tracer velocities are closer to 
the true value than the conductivity determined values. 



5 . 4 SUMMARY 

Bromide was a good choice for tracing water movement in the 
wetland studied. Using the natural flow tracer technique, 
sufficient bromide concentrations were detected up to a distance 
of approximately 38 m from tracer input. The bromide specific 
ion electrode method of analysis was found to be adequate. This 
method could be improved by periodic double sampling and analysis 
with another method (titration) especially at low bromide 
concentrations . 

The wide range of velocities determined from this study are 
the effect of a very heterogeneous wetland flow system and the 
methodology used in determining the velocities. Some of the 
heterogeneity involved a change in hydraulic gradient, changing 
flow rates and periodic flushing over time. Preferential flow 
paths and "dead zones", the result of variable hydraulic 
conductivity, also add to the heterogeneity. The most represen- 
tative velocity of the AMD water in, the wetland was that 

determined by the tracer experiment. These velocities were found 
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to be a reflection of the heterogeneous system and ranged between 
3.1 X lO"-"- to 6.2 X 10"4 cm/s with a mean of 2.6 x 10"^ cm/s 
(73.7 ft/day) . 

Comparison of field determined hydraulic conductivity, from 
tests of auger holes and wells, to tracer determined velocity was 
accomplished by converting conductivity to velocity utilizing the 
estimated values of effective porosity and hydraulic gradient. 
These field conductivity determined velocities are reasonable in 
that they estimate water flow velocity in the small area where 
conductivity was determined. When comparing the conductivity 
determined velocity to tracer velocity even the average velocity 
of auger hole sites (fastest of the conductivity determined 
velocities) was three orders of magnitude lower than average 
tracer velocity. This order of magnitude relationship does not 
change even if effective porosity is changed to 0.1 (low) and the 
hydraulic gradient is changed to the steepest observed in the 
study area (0.058 m/m) . This relationship infers the comparison 
is that of two different flow systems. The tracer velocity 
represents preferential flow paths within the shallow (acrotelm) 
system and the conductivity represents that flow found at greater 
depth (catotelm) . 

The determination of water residence times in the wetland is 
contingent upon the velocity of the water system to be considered 
and the distance the water travels. The concern in this study 
involves the AMD water which emanates from the culvert outflow 
into the wetland. The velocity representing this source is that 
of the tracer which was input at the same point that the AMD 
water is introduced to the wetland. The distance the AMD water 
travels, before complete problem amelioration, was not deter- 
minable from the tracer methods employed, but is estimated to be 
approximately 200m as observed in acrotelm soil iron concentra- 
tions (Appendix C and Section 8) . The shallow groundwater iron 
concentrations show a trend similar to the soils by decreasing in 
the down-gradient direction and approaching background iron 
concentrations near site F-1. With the 200m distance chosen and 
using the mean tracer velocity, the average residence time of AMD 
water in contact with the wetland acrotelm soils is estimated to 
be 215 hours (9 days) . This estimate is considered to be quite 
variable as noted earlier and represents a low flow regime 
(Autumn) . This residence time estimate should be viewed with 
caution as it pertains only to iron amelioration and the velocity 
used represents only that measured at the time of the tracer 
study for the very limited area of the wetland near the AMD 
discharge point. Assuming the wetland studied is a viable system 
for the amelioration of AMD, the challenge in the design of man 
made wetland systems will be to determine how to emplace similar 
hydrologic characteristics, and limit preferential flow paths. 
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6.0 WETLAND INFLUENCE ON MINE DRAINAGE WATER QUALITY 

The water quality sampling program for the Swamp Gulch 
study was initiated to determine: 1) the concentrations of 
parameters in the AMD; 2) the background water quality in Swamp 
Gulch above the Carbonate Mine and on the Blackfoot River above 
the wetland; 3) the changes in AMD water quality as it pro- 
gressed through the wetland and; 4) changes in the water quality 
of the Blackfoot River below the wetland discharge point. 



6.1 SAMPLE COLLECTION AND PROCESSING PROCEDURES 

Water samples were collected by various means, each ap- 
propriate to the sample source. All samples were initially 
collected in 4 liter (1 gal) linear polyethylene (LPE) (Nalgene) 
decontaminated bottles (acid/deionized water cleaned) . These 
bottles were rinsed three times with sample water prior to fill- 
ing. The filled 4 liter bottles were placed in coolers with ice 
until divided into appropriate aliquots for preservation, 
usually within minutes after collection. 

Water level depths were recorded at all wells prior to 
pumping. These measurements were referred to the top of the 
well casing (as surveyed) and recorded to the nearest 0.01 foot 
using a Soil Test Model 444-000 or Johnson UOP water level 
indicator. 

Well samples were acquired with a peristaltic pump (Master- 
flex) using a length of silicone tube. The apparatus was rinsed 
with distilled or deionized water after each sample, and if 
organic deposition on the hose was evident, new hose was in- 
stalled. A minimum of three well volumes (including sand pack 
volume) of water was pumped before sampling at the higher yield 
sites. At low yield sites, wells were pumped dry and the re- 
charge was sampled. 

Surface water samples were obtained by submersing the four 
liter bottles below the water surface where sufficient water 
depth was present or by filling the four liter bottles with a 
peristaltic pump where small volumes of water required this 
method. Flow measurements or estimates were made at the time of 
sample collection for all flowing surface water sites. 

Various aliquot portions were taken from the four liter 
bottles for the required analyses. The samples were preserved 
in three separate bottles: 1) 125 or 250 ml of 0.45 micron 
filtered sample with 1 ml of concentrated HNO3 preservative 
(metal analyses), 2) 125 or 250 ml of 0.45 micron filtered 
sample with one ml of concentrated H2SO4 preservative (NO3 
analysis) and 3) 250 or 500 ml of nonfiltered sample with no 
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added preservative (pH, electrical conductivity (EC) , major 
cations, HCO3, F, P, CI, SO4 analyses). Special bottles were 
prepared for the limited number of dissolved organic carbon and 
cyanide analyses. 

Samples were filtered utilizing a Horizon Ecology peristal- 
tic pump (silicone hose) and a Geotech 102 mm filter apparatus. 
The pump and filter assembly were rinsed thoroughly with deio- 
nized water between samples. Filters and pref liters were 
changed between all samples except field replicates. All sample 
bottles were rinsed three times with sample water (filtered or 
nonfiltered according to aliquot) prior to filling. A minimum 
of one field blank was included with each sample collection and 
inserted blind within the sample train. The blank was processed 
in a manner identical to that of the normal samples. Blind 
field replicates were processed from the same four liter initial 
sample. Blind EPA field standards were processed in the same 
manner as were field blanks. All processed samples were immedi- 
ately placed in coolers with ice or refrigerated. 

Temperature, conductivity (EC) and pH measurements were 
made (in-situ where possible) at the time of sample collection 
or, under adverse conditions (ambient temperatures below freez- 
ing) , they were determined within one to three hours in a 
sheltered environment. The temperature was determined with a 
mercury thermometer, dial thermometer or a Cole-Parmer SL-5985- 
80 pH meter with a digital temperature readout. These instru- 
ments were evaluated with an ice bath and a laboratory mercury 
thermometer for accuracy and, in all cases, were accurate within 
1°C. Sample pH was determined using either a Cole-Parmer Model 
5986-00 or a Cole Parmer Model SL-5985-80 pH meter. The latter 
instrument was rated to maintain calibration up to one year. It 
was calibrated at the beginning of each sample trip and checked 
daily. The instrument was recalibrated if the pH varied more 
than 0.2 pH units from standards. The Model 5986-00 pH meter 
was standardized for each sample. 

The EC of samples was measured with a Lamotte Chemical 
Model DA-1 Conductivity Bridge or a Cole-Parmer Model R-1481-50 
Conductivity Meter. The Lamotte unit was standardized for each 
analysis. The Cole-Parmer unit was standardized at the begin- 
ning of each sample trip and checked twice daily thereafter. 
This unit is rated to maintain calibration up to 90 days. 

All data were recorded in the project field book or on site 
matrix sheets. 

At the completion of sample collection and processing, 
samples were stored under security until shipped by bus to 
Energy Laboratories in Billings, Montana. Although strict EPA 
custody procedures were not followed due to the expense, at no 
time were samples left out of sight, unattended or unsecured. 
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The water sampling program for the Swamp Gulch wetland 
study began with initial sampling during the April 29 - May 2, 
1987 period, following completion of the project installation. 
The largest surface flows through the wetland were observed at 
this time, but the flume on Swamp Gulch had not yet been instal- 
led so only a one point estimate of Swamp Gulch flow was noted 
(estimated at 20 to 30 gpm) . The second set of samples was 
collected from July 28 to July 31, 1987, and included samples 
from piezometers F-l-S and F-l-D which were installed on July 
28, 1987. The third set of samples, collected during the 
January 11 to January 13, 1988 period, consisted of a partial 
set confined largely to the area of estimated maximum acid 
drainage impact with several appropriate background sites. Low 
flow conditions prevailed during both the second and third 
sampling trips. 

6.2 ANALYSES METHODS 

All analyses were made utilizing standard EPA procedures 
(Appendix A) . These are detailed in Methods for Chemical 
Analysis of Water and Wastes (EPA 1983). 

6.3 QUALITY CONTROL/ QUALITY ASSURANCE 

The water quality program was designed to evaluate the 
validity of these data. Field sampling included blind field 
replicates, blind field water blanks and blind field standards. 
Laboratory procedures included duplicates and spikes. These 
data have been utilized to generate accuracy and precision 
statements following standard procedures for evaluating data 
from the EPA Contract Laboratories Program (Appendix A) . These 
data are summarized in the following sections. 

6.3.1 Accuracv 

An evaluation of five Environmental Protection Agency (EPA) 
standards submitted with samples indicate laboratory analyses 
somewhat overestimated Ca, Mg, SO4 , HCO3, TDS and NO3 (Table 6- 
1) . The laboratory did not recover any phosphate as ortho- 
phosphate from standards containing 0.50 and 0.27 mg/L 
phosphate. This problem is the result of running phosphate 
analyses on the unpreserved aliquot which did not contain the 
PO4 standard. The laboratory spikes for PO4 suggest an adequate 
recovery (96 percent) for this parameter. Three EPA standards 
were submitted for trace element analyses. These data generally 
exhibit low standard deviation values and little bias (Table 6- 
2). The variation between the number of samples was due to the 
changing of the suite of analyses to be performed as the project 
continued. 
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6.3.2 Precision 

Precision calculations were based on blind field replicates 
and laboratory duplicates. The maximum range of uncertainty for 
laboratory duplicates (at the 90% confidence level) was 7.69% 
for HCO3 (Table 6-3) . Precision calculations based on blind 
field replicates generally indicated good precision (at the 90% 
confidence level) except for sodium, fluoride, manganese and 
total suspended solids (TSS) (Table 6-4) . The high imprecision 
for sodium is the result of values at or near the instrument 
detection limit where a 1 to 2 ppm difference in analyses 
produced large relative standard deviations. High uncertainty 
values for fluoride and manganese were due to the small number 
of valid samples (2 each) and the resulting high t values. The 
precision of total suspended solids (TSS) analyses based on 
blind field replicates was strongly influence by one sample 
(from a total of 2) . The TSS sample and replicate sampled on 
April 29, 1987, exhibited values of 16 and 8 ppm, respectively, 
for this parameter which resulted in a 67 percent relative 
standard deviation and a precision value of +298 percent. 

6.3.3 Cross Contamination 

Evaluation of the field water blanks indicated minimal 
cross contamination (Table 6-5) . No chemical parameters were 
above the instrument detection limit except bicarbonate for the 
April 1987 and January 1988 sampling trips. Some bicarbonate is 
usually present in all blanks due to carbon dioxide present in 
the atmosphere, a portion of which dissolves in water forming 
bicarbonate. One out of the two blanks for the July 1987 trip 
exhibited a sulfate value of 40 mg/L, likely as the result of 
sulfuric acid (H2SO4) used for preservation. The bottles may 
have been mislabeled or the laboratory may have run sulfate on 
the wrong aliquot. It is not felt that an actual contamination 
by sulfate occurred as all other parameters associated with this 
blank, with the exception of TDS, were below the detection 
limits. Sulfate values associated with normal water samples for 
the July 1987 sample trip do not show any associated jump in 
concentration . 



6.4 SURFACE WATER QUALITY 

6.4.1 Swamp Gulch Above the Carbonate Mine 

The water quality of Swamp Gulch above the mine is good 
and meets most criteria for drinking water (USPHS 1962; EPA 
1978, 1979), Livestock (NAS 1974), aquatic life (EPA 1983; NAS 
1974) and irrigation (NAS 1974) . The only exception was an 
enormously high (0.21 mg/L) manganese value obtained during the 
January 1988 sample collection. The April and July, 1987 values 
for manganese were both below the 0.02 mg/L instrument detection 
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Table 6-4. 



Precision, based on field replicates, at the 90% 
confidence level. 



Parameter 


Precision 


Completeness 


Comments 


Ca 


±0.0% 


100% 




Mg 


± 4.4% 


100% 


— 


Na 


±46.0% 


100% 


— 


K 


+ 0.0% 


50% 


— 


SO4 


± 1.2% 


75% 


— 


HCO3 


±6.0% 


100% 


- 


CI 


±0.0% 


50% 


- 


F 


±81.0% 


50% 


— 


P 


NC 


0% 


All values <IDL 


NO3 


±13.6% 


75% 


- 


Al 


NC 


0% 


All values <IDL 


Cd 


NC 


0% 


All values <IDL 


Fe 


5.6% 


75% 


— 


Pb 


NC 


0% 


All values <IDL 


Mn 


27% 


50% 


- 


PH 


2.2% 


100% 


.. 


EC 


2.5% 


100% 


— 


TDS 


8.0% 


100% 


— 


TSS 


298.0% 


67% 





NC 



Not calculated 



limit. The high sample value may have been due to the amount of 
sediment disturbed during sampling (ice removal) . Swamp Gulch 
water above the mine has exhibited little variation during the 
period of study and is characterized by a calcium bicarbonate 
water type. Calcium, magnesium, bicarbonate and sulfate percen- 
tages in eguivalents/L are 54.6, 37.9, 80.4 and 19.6, respec- 
tively. Slight increases were noted for both sodium (2 mg/L 
increase) and sulfate (7 mg/L increase) during the January 1988 
sample period and are of little significance. 

6.4.2 Swamp Gulch Below the Carbonate Mine 

Swamp Gulch water below the mine consistently exceeded 
drinking water criteria (USPHS 1962, EPA 1978, EPA 1983) for 
aluminum, cadmium, iron and manganese. The water was charac- 
terized by an iron sulfate water type in which iron and sulfate 
comprised 31.1 and 99+ percent of cations and anions, 
respectively, assuming there was an even balance between ferrous 
and ferric iron forms. Other major cations were calcium (27.6 
percent), magnesium (18.7 percent), aluminum (16.9 percent) , 
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manganese (3.3 percent) and sodium (2.3 percent). This water 
also exhibited an increase in sulfate and iron in the January 
1988 sampling. The increase in iron was especially notable 
(from 18.8 to 48.4 mg/L) and likely represents the combined 
result of several factors: 1) low flow which has less dilution; 
2) the presence of an ice cover, which may have restricted 
oxygen availability (and subsequent precipitation of iron as 
iron-hydroxide) and, 3) somewhat lower pH which would allow more 
iron in solution. This observation suggests physio-chemical 
precipitation of iron may be important in this system and that 
this process starts before the AMD is discharged into the 
wetland. 

6.4.3 Blackfoot River 

The general water type of the Blackfoot River at Pass 
Creek, located upstream of the wetland study area, was consis- 
tent throughout all three sample periods (Figure 3-1) . The 
dominant water type, expressed in milliequivalent per liter, was 
a calcium bicarbonate in which calcium and bicarbonate composed 
53 percent of both cations and anions respectively. Magnesium 
and sulfate represented 44 and 47 percent of cations and anions. 
Little change occurred in these percentages during high or low 
flow. The water type for the Blackfoot River at the Meadow 
Creek road bridge, located downstream of the wetland study area, 
exhibited a marked change during low flow. The April and July 
1987 water samples exhibited a similar cation-anion balance to 
that found at the Pass Creek site (54, 43, 52 and 48 percent for 
calcium, magnesium, bicarbonate and sulfate, respectively) . 
However, during the January 1988 sample period, a slight 
decrease in calcium (to 51 percent), an increase in sodium (from 
2.1 to 4.6 percent) and a gross increase in sulfate (69.7 
percent) at the expense of bicarbonate (30.3 percent) was 
evident. The ramifications of this will be discussed later. 

6.4.4 Background Statistical Comparisons 

Two general comparisons have been evaluated: 1) the Swamp 
Gulch water quality above and below the Carbonate Mine and 2) 
the Blackfoot River above and below the wetland discharge area. 
The Swamp Gulch drainage exhibited significant to very highly 
significant increases in the concentrations of SO4, Al, Mn, Zn, 
Cd, Fe, Pb and Cu when compared using a paired "t" test. No 
significant increases were observed for As, Cu, Co, Sb, Se and 
Ag at the given detection limits (Appendix B-3) . The "t" test 
comparison of the Blackfoot River at Pass Creek and at the 
Meadow Creek Road Bridge indicated no significant differences. 
However, the limited number of samples (3) for each site did 
suggest the Meadow Creek Road Bridge sample site was slightly 
higher in iron, manganese and sulfate. The slight increases in 
parameters in Blackfoot River samples below the Swamp Gulch 
wetland discharge point as opposed to those above the wetland 
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are great enough to account for approximately one half of the 
iron and all of the manganese contributed by the Swamp Gulch 
acid mine drainage. Zinc is the only other trace metal present 
in the Blackfoot River at notable levels and was apparently 
diluted in the river reach adjacent to the wetland. Both the 
April and July sample sets show the highest zinc levels in the 
upper two river sample sites with each downstream sample site 
showing progressively lower zinc concentrations. This effect 
could not be determined during January 1988 due to the reduced 
number of samples taken at that time. These data strongly 
suggest the source of the zinc is either Paymaster Gulch, the 
main fork of the Blackfoot River immediately above the study 
area or some discharge directly into the river or wetland . The 
latter may be some flowing abandoned drill holes, one of which 
was measured for EC and pH (352 umbos/ cm and 3.97, respectively) 
on April 29, 1987. 

6.4.5 Wetland Surface Water 

In addition to the Blackfoot River and Swamp Gulch samples, 
surface water samples were obtained from six surface sites 
within the wetland during the July sample trip. These were 
obtained at sites AAA-3, BC-1, C-1.5, C-3, CD-I and F-1. Water 
quality parameters at site AAA-3 were very similar to the 
Blackfoot River at Pass Creek samples and exhibited a calcium 
bicarbonate water type and low levels of all metals (iron at 
0.07 mg/L) . This would be expected from the observed drainage 
system. Parameters at site C-3 were also similar to the river 
samples but iron and manganese levels were elevated to 6.94 and 
0.24 mg/L, respectively. This suggests some wetland drainage is 
mixing with the diverted river water at this point. The pond 
sample at BC-1 was intermediate in composition between values 
observed at C-3 and the Swamp Gulch drainage which suggests this 
pond is being partially recharged from the upgradient, non- 
polluted wetland area. The water quality at sites C-1.5, CD-I 
and F-1 was similar to the Swamp Gulch AMD discharge. The water 
type at C-1.5 was an iron sulfate and at sites CD-I and F-1 it 
was a calcium sulfate. Sulfate was the only major anion in the 
latter samples. High levels of sulfate, iron and manganese 
would be expected at sites C-1.5 and CD-I, located 37.5m (123 
ft) and 32.9m (108 ft), respectively, from the Swamp Gulch 
discharge point. However, the high levels observed at site F-1 
raise serious questions as to the actual effectiveness of this 
wetland in ameliorating AMD water. The process may be one of 
partial dilution which cannot be accurately accessed with the 
present instrumentation. It should also be noted that these 
conclusions are based on one set of surface samples which could 
be subject to analytical or sampling errors which are not 
apparent . 
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6.5 GROUND WATER QUALITY 

The ground water of the wetland research area has been 
characterized by samples from two zones, shallow and deep. The 
former is generally from 0.6 to 1.2 m (2 to 4 ft) and the latter 
varied up to a maximum of 4.9 m (16 ft). The shallow zone 
screened interval generally spanned portions of both the 
surf icial organic layer (acrotelm) and the underlying muck layer 
(catotelm) . The deep wells bottomed in the gravel underlying 
the muck layer. However, extreme difficulty was encountered 
trying to keep the deep bore holes open long enough to place the 
screen within the gravel zone and as a consequence, most deep 
wells have a significant portion of their screened interval 
within the muck zone. All wells were sandpacked and, on the 
deep wells, the sand pack extended into the partially caved 
gravel portion of the bores. The water quality data strongly 
suggest most water withdrawn from the deep wells is from .the 
gravel aquifer. The following sections describe water quality 
parameters of the shallow and deep zones. An evaluation of the 
two wells in the acid drainage seep have been evaluated 
separately in section 6.5.3. 

6.5.1 Shallow Ground-water System 

The shallow ground-water system (0-1. 2m, 0-4ft) is 
dominated by iron-calcium sulfate and calcium-magnesium 
bicarbonate water types. The former is present throughout the 
main Swamp Gulch AMD impact area, the latter is generally 
confined to the southwestern portions of the study area (Figure 
6""1) • A calcium-magnesium sulfate water type is conspicuously 
present at two sites near the Blackfoot River (C-5 and F-3) . 
The only other site at which the calcium-magnesium sulfate water 
type was dominant is the site AA-2. It is noteworthy that 
background site AAA-3 exhibits the same water type found in the 
AMD impact area although overall concentrations are lower at 
site AAA-3. The zone dominated by the bicarbonate anion is 
confined to the area which apparently receives water from the 
Blackfoot River via channels flowing past sites c-3 and E-3. 
The single sample from site F-1 was the only analyses in which 
bicarbonate was the dominant anion and which was not affected by 
Blackfoot River water. 

A plot of mean dissolved iron values for all well sites 
clearly indicates that the shallow ground-water system is 
enriched in iron as the result of AMD (Figure 6-2). Background 
means for iron ranged from <0.03 to 27 mg/L with the very low 
values (<0.03 and 0.25 mg/L) associated with sites E-3 and F-3, 
respectively. Sites E-3 and F-3 are apparently recharged 
directly from the Blackfoot River. Samples from shallow well 
sites near AMD discharge point are 2 to 3 times above both the 
background shallow ground-water concentrations and the inflow 
AMD concentrations. The highest mean (73 mg/L) was found for 
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Figure 6-1. Shallow ground-water types for the wetland system. 

site C-1, the closest site to the AMD discharge i oint. 
Dissolved iron concentrations fall abruptly e.. ard towards 
site B-2, at least partially the result of lowei ground-water 
gradients in this direction. Concentrations fall more gently 
towards the southwest, which suggests more AMD Hows in this 
direction which is consistent with topography and piezometric 
maps. Iron concentrations are apparently reduced 0.65 to 0.98 
percent per meter (0.2 to 0.3 percent per foot) away from the 
maximum concentration observed at site C~l in this ^^etland. It 
is apparent that much of the iron in the AMD hp.s been deposited 
in the 3.9 ha (9.6 ac) area surrounding the Ai€D discharge point. 

Isolines of dissolved manganese concentrations in the 
shallow ground-water system show the highest concentrations are 
associated with sites CD-I and D-1 (Figure 6-3). The zone of 
elevated manganese concentrations is generally confined to the 
area between the F-2 to C-2 line and Highway 200 or roughly 1 ha 
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Figure 6-2. Iron concentrations for the wetland shallow ground- 
water system. 



(2.5 ac) . Concentrations in this zone are 3 to 11 times 
background concentrations. High manganese concentrations have 
also been observed at sites C-5 and E-3 but it is doubtful that 
these sites are impacted from the Swamp Gulch AMD. 

Manganese is apparently attenuated to a much lesser degree 
than is iron in this wetland. The manganese concentration at 
site F-1, near the wetland outflow, are still about 4 times 
above background levels. These results are consistent with 
other studies which have shown Sphagnum wetlands are more 
effective in removing iron than they are for removing manganese 
from AMD (Wieder and Lang 1986, Burris et al. 1984). 
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Figure 6-3. 



Manganese concentrations for the ^';pt and shallow 
ground-water system. 



Isolines of sulfate concentrations follow a pattern very 
similar to that found for the bromide tracer (Figure 6-4) . From 
a high concentration zone in the D-1 to C-1 area, levels 
decrease gently towards sites E-1 and B-2, apparently along the 
preferential flow paths. Sharp concentration reductions over 
short distances occur at the margins of the flow paths. 
Concentrations are typically 3 times above background levels in 
the 0.7 ha (1.8 ac) AMD impacted zone. The sulfate isolines 
suggest a low level of AMD impact does occur in the area defined 
by site B-2, A-3, C-3 and C-2. This is supported by the 
increase in iron and manganese concentrations observed in the 



79 



CARBONATE MINEX- 




A WELL SITES 

O SURFACE WATER SAMPLE 

SITE 
D SURVEY CONTROL STATION 
H FLUME 

O PRECIPITATION GAGE 
— ISOLINES OF SO4 CONCENTRATIONS 
ftvg/L) 



SEC20, TI5N,R6W 
LEWIS AND CLARK CO. .MONTANA 

FEET 
200 



Figure 6-4. 



Sulfate concentrations for the wetland shallow 
ground-water system. 



surface water sample taken at site C-3 (section 6.4.5). This 
wetland is apparently effective in attenuating shallow ground- 
water sulfate concentrations to levels approaching background 
values. 

Considerable variation was apparent in TDS concentrations 
within the wetland. Mean values ranged from a low of 119 mg/L 
at site D-2 to a high of 403 mg/L at site D-1. Higher TDS con- 
centrations were generally confined to a 0.6 ha (1.6 ac) area in 
the same general vicinity where elevated manganese 
concentrations were found (Figure 6-5) . Total dissolved solids 
concentrations are reduced from values about 2 times above 
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Figure 6-5. 



Total dissolved solids concentratio. 
wetland shallow ground-water systr" 



for the 



background values to levels similar to background concentrations 
within the wetland area. 

Concentrations of aluminum, cadmium, lead and zinc were 
generally low throughout the wetland. The maximum aluminum 
value was 2.8 mg/L in well BC-l-S for April, 1987 saiaple trip. 
Aluminum was detected in 18 other samples from 7 other wells 
which all but one sample exhibited values lass than 1 mg/L. 
Only two samples contained cadmium in excess of 0.01 mg/L, the 
recoxamended criteria for potable and irrigation water supplies. 
The two samples were obtained from sites BC-1-2 (0.032 mg/L) and 
BC-l-S (0.044 mg/L) both of which are in close proximity to the 
Swamp Gulch AMD discharge point. These two wells were also the 
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source of the only two samples in which detectable levels of 
lead were found (both 0.03 mg/L) . All zinc concentrations were 
well below applicable criteria for all common uses. The maximum 
zinc concentration in the impacted area was 0.49 mg/L from well 
D-l-S for the April, 1987 sample trip. The maximum zinc value 
observed in the wetland was 1.18 for site C-5-S, but this site 
is not likely affected by the Swamp Gulch AMD. 

Numerous other trace metals were analyzed in the initial 
April 1987 samples. These included arsenic, chromium, cobalt, 
copper, nickel, antimony, tin, selenium and silver. No shallow 
wells had values above detection limits for arsenic, chromium, 
copper, nickel, selenium and silver. Only site C-5-S had 
detectable cobalt at 0.09 mg/L. Antimony concentrations at or 
above the 0.01 detection limit were found at 14 sites. The 
highest antimony levels found were 0.03 mg/L at BC-1-2, D-l-S 
and C-5-S. These data suggest there is some enrichment of 
antimony due to the Swamp Gulch AMD, but the effect is minor. 

In summary, the levels of aluminum, cadmium, lead, zinc, 
sulfate and TDS observed in the wetland shallow ground water 
system should have little overall adverse impact on vegetation 
or wildlife. However, the high iron and manganese 
concentrations observed exceed all applicable water use criteria 
and can be expected to adversely impact all species not adapted 
to high levels of these metals. The elevated cadmium 
concentrations observed near the AMD discharge point may also be 
of concern and may be detrimental to wildlife through 
bioaccumulation (Section 9) . 

6.5.2 Deep Ground-water Svstem 

The overall water type for the deep ground-water system 
exhibited markedly less variation than the shallow ground-water 
system. A calcium-magnesiiom sulfate water type was dominant in 
all sites except AAA-3, B-2 and C-2, which were characterized by 
an iron-calcium sulfate water type, and sites D-1 and F-1 in 
which calcium-iron bicarbonate and calcium-magnesium bicarbonate 
dominated, respectively. Site C-3-D exhibited a sodium-iron 
sulfate water type which was likely the result of problems with 
bentonite contamination during construction. 

Isoline map evaluations for concentrations of iron, man- 
ganese, sulfate and TDS suggest the AMD has a much more limited 
impact on this aquifer system. Only iron and manganese show 
increases in concentrations associated with AMD and the effect 
is generally subdued as compared to the shallow aquifer system 
(Figures 6-6 and 6-7). The areas in which iron and manganese 
levels are clearly elevated above background concentrations 
roughly coincide with the impact area observed for elevated TDS 
levels in the shallow ground-water system and cover approxi- 
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Figure 6-6. 



Iron concentrations for the wetland deep ground- 
water system. 



mately 1.5 ha (3.6 ac) and 1.7 ha (4.3 ac) , respectively. 
Concentrations and attenuation gradients in the deep ground- 
water system for both iron and manganese are considerably less 
than in the shallow ground-water system which suggests a very 
limited AMD input to the deep system. This conclusion is 
supported by sulfate and TDS concentrations which were clearly 
elevated above background levels only at wf^ll C-l-D. 

No deep wells exhibited lead values above the 0.01 mg/L 
detection limit. The highest wetland aluminum values from deep 
wells were found at sites AA-2 and B-2 with maximums of 1.7 and 



83 



CARBONATE MINE/^ 




A WELL SITE 

O SURFACE WATER SAMPLE 

SITE 
D SURVEY CONTROL STATION 

•m FLUME 

PRECIPITATION GAGE 
iSOLINES OF MAN6ANESE(mg/L) 



SEC 20, T15N.R6W 
LEWIS AND CLARK CO .MONTANA 

FEET 
200 



Figure 6-7. 



Manganese concentrations for the wetland deep 
ground-water system. 



3.0 mg/L, respectively. No other sites were found with aluminum 
values greater than 1.0 mg/L and many were below the 0.1 mg/L 
detection limit. Detectable cadmium concentrations were found 
at sites AAA-3, AA-2, B-2, BC-1, C-1 and C-2 . Maximum values 
were consistently found for site AA-2 with a mean and maximum of 
0.007 and 0.012 mg/L, respectively. This was the only site 
exhibiting levels above 0.01 mg/L, a recommended criteria level 
for potable and irrigation water. As with shallow ground -water 
samples, arsenic, chromium, cobalt, copper, nickel, antimony, 
selenium and silver were determined on the April 1987 sample 
set. Arsenic was detected at only two sites, AAA-3 at 0.012 
mg/L and AA-2 at 0.005 mg/L. Cobalt was detected only at site 
C-5 which contained 0.02 mg/L. Copper was present in detectable 
quantities only at sites AA-2 (0.02 mg/L) and B-2 (0.07 mg/L). 
The maximum antimony concentration found (0.03 mg/L) was for 
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site C-1. Sites C-5 and D-1 both exhibited antiaony values of 
0.02 mg/L and all other sites were at or below the 0.01 mg/L 
detection limit for antimony. No detectPble amounts were found 
for chromium, nickel, selenium and silver in the deep wells. 
With the exception of antimony, the secondary trace metals have 
no apparent relationship with the AMD. 

It is possible that some of the observed effect of elevated 
concentrations of iron and manganese in deep wells is due to 
subsurface seepage from old mine workings that migrates below 
the highway in fractured bedrock rather than vertical perco- 
lation through the wetland muck layer. There is insufficient 
data at present to evaluate this factor. 

6.5.3 Acid Mine Drainage Seepage Wells 

The water quality of the shallow and deep "seep" wells is 
important as it represents the actual water quality that is 
likely emanating from old mine workings prior to dilution by 
surface water (Figure 5-9 and 5-10) . The general water type for 
this site is a calcium sulfate with only small to moderate 
differences noted for most major parameters between deep and 
shallow samples. 

This water was characterized by low pH (3.0) and moderate 
to high levels of most trace metals. Mean values for TDS, EC 
and sulfate are 1838 mg/L, 1924 umhos/cm and 1228 mg/L, 
respectively. Parameters which did exhibit notable differences 
between shallow and deep samples were aluminum, iron, manganese, 
zinc and several minor metals. Aluminum and zinc were 
consistently higher in shallow samples, 70.4 mg/L versus 3.6 
mg/L and 31.9 mg/L versus 8.5 mg/L, respectively. This trend 
was also noted for arsenic, cadmium, copper, nickel and 
antimony. Both iron and manganese were higher in the deep well 
than in the shallow well (57.5 versus 5.0 mg/L a i 65.1 versus 
40 mg/L, respectively) . Silver was also apparently present at 
slightly higher concentrations in the deeper V'^" " at this site. 
It is readily apparent that many of the trace el-ment 
concentrations present in water at these two wells are quickly 
diluted by high quality upper Swamp Gulch water before entering 
the wetland. This process likely eliminates most of the 
toxicity problems that might otherwise occur. The high levels 
of many metals, sulfate and fluoride would make this water 
unsuitable for any common use including potable supplies, 
livestock, irrigation and aquatic life usage. 

6.6 SUMMARY 

Of the parameters analyzed (Appendix B) , only aluminum, 
cadmium, copper, iron, lead, manganese, zinc and sulfate were 
significantly elevated in the Swamp Gulch watershed drainage due 
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to AMD. Concentrations of these parameters in the wetland 
shallow ground-water system are elevated to levels above the 
concentrations found in the AMD near the AMD inflow point and 
decrease both laterally and vertically from this point. The 
only catotelm area notably affected by AMD is in the immediate 
vicinity (30 m - 100 ft) of the AMD discharge point and, in 
general, this stratigraphic unit is only minimally impacted by 
the AMD. 

Minor changes in the water quality of the Blackfoot River 
above and below the wetland discharge point suggest that 
manganese and possibly other several metals may not be 
effectively ameliorated within the wetland. However, the Swamp 
Gulch wetland area is at least partially effective in removing 
iron from the AMD. These data suggest about half of the 21 mg/L 
weighted mean iron concentration in the mean 45m3 daily AMD 
inflow may be precipitated within the 3.9 ha (9.6 acre) iron 
impacted acrotelm area (Section 8). These data suggest some 
means of physio-chemical precipitation of iron should be 
incorporated in artificial wetland construction. This system 
should provide aeration of the AMD to enhance iron 
precipitation, possibly through the use of waterfalls. Since 
optimum treatment of AMD pH and sulfate by microbial action 
(Section 10) requires reduced conditions, the aeration stage 
should be constructed near the discharge end of the artificial 
wetland, under the assumption that it will be much easier to add 
oxygen to the AMD than it would be to remove it. 
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7.0 DISSOLVED METAL LOADING OF A NATORAT. ..iTLAND 



An estimation of the AMD metal load supplied to the Swamp 
Gulch wetland through the April 1987 - March 1988 year has been 
made based on water sample concentrations and estimated flow 
volumes. A comparison of this loading with observed metal 
masses in the acrotelm (Section 8) permit an estimation of past 
wetland efficiency in amelioration of AMD. 



7 . 1 METHODS 

The calculation of metal loading into the wetland was 
hampered by the short duration of runoff records necessitating 
estimations be made for a number of factors. The assumptions 
that have been made are: 

o the Swamp Gulch input represents all metal loading to 
the study area; 

o the three sample dates, April 29, July 29/ 1987 and 
January 11, 1988, are representative of all flows for 
the periods April 1 through June 30, 1987, July 1 
through November 16, 1987 and November 17, 1987 through 
March 31, 1988, respectively; 

o the observed flow on April 29, 1987 is representative of 
the entire month; and 

o that the mean daily flow rate for the period May 1 
through June 15, 1987, was equal to the value for the 
period June 15 through June 30, 1987. 

The latter assumption is based on the Lincoln Ra: jer Station 
precipitation which was similar for the two montl.s at 4.1 cm and 
3.4 cm (1.62 and 1.33 inches) for May and June, t aspectively. 
The Lincoln Ranger Station records also indicate 3.99 cm (0.39 
in) of precipitation fell during April, 1987, and that 
precipitation for December, 1986, January and February 1987 were 
far below normal. Casual field observations during April and 
May indicated little snow melt runoff occurred and thus the 
assumption was made that the April 1987 flow estimate was 
reasonable for the entire month. The time intervals were 
determined on the general flow regime and the time of sampling. 
These are, in general, spring flow consisting of some snow melt 
and precipitation; summer flow consisting of base flow and 
precipitation runoff; and, winter low flow (base flow) when no 
runoff and little precipitation input occur and during which 
freezing conditions persist. 
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The wetland metal loading figures are thus basou on a 
limited record for a period in which abnormally low precipi- 
tation conditions prevailed. Extrapolation of these data to 
cover the annual period of April 1987 through March 1988 may 
involve appreciable error and further extrapolation to other 
years is not generally recommended. 



7.2 RESULTS AND DISCUSSION 

The approximate weights (mass) of dissolved metals dis- 
charged via Swamp Gulch into the wetland from April 1987 through 
March 1988 varied from 338 kg (744 lbs) for iron to values less 
than 100 grams (3.5 oz) for arsenic and silver (Table 7-1). 
High values for aluminum (107 kg) , manganese (63 kg) , zinc (44 
kg) and copper (12 kg) are also evident. All other metal 
loading values were less than 1 kg and of these, only cadmium 
with its ability to be readily bioaccumulated, may be of 
additional concern. The implications and fate of these metals 
is discussed further in section 8. 



Table 7-1. 



Total dissolved metal loading (kg) on the Swamp Gulch 
wetland, April 1987 through March 1988. 



Apr. 1- June 30 
Flow Volume fm-^) 7.940 



PERIOD ■ 

July l-Nov.16 
7.360 



Nov.l7-March Annual 
1.110 



Parameter 

Al 

As 

Cd 

Cr 

Co 

Cu 

Fe 

Pb 

Mn 

Ni 

Sb 

Se 

Ag 

Zn 



48 


52 


6.1 


107 


<0.04 


<0.04 


<0.01 


<0.09 


0.302 


0.213 


0.027 


0.542 


<0.159 





- — 





<0.079 


<0.074 


<0.011 


<0.164 


6.274 


4.931 


410 


11.6 


46 


138 


5' 


338 


.278* 


0.221 


0.044 


.543* 


30.8 


28.8 


3.75 


63.4 


0.32 


0.29 


0.03 


0.64 


0.08 


<.07 


<0.06 


0.08 


<0.04 





. 





<0.04 


<0.04 


<.01 


<0.09 


22.6 


17.9 


2.9 est. 


43.6 



*Estimate 



Annual Total 565 kg 
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8.0 WETLAND GEOCHEMISTRY AND THE CONTROL OF ACID HINE DRAINAGE 

The Swamp Gulch wetland area core samples have been 
analyzed for numerous chemical and physical parameters (Appendix 
C) . These analyses were made to determine 1) the extent of past 
metal loading; 2) the distribution of various parameters with 
depth and stratigraphic changes; and 3) the basis for comparison 
of elemental concentrations observed in vegetation growing on 
this site. 

8.1 METHODS 

Samples for geochemical evaluation were obtained during 
construction of well and auger hole installations. Core samples 
were obtained by methods outlined in section 5.2.1. These 
samples were placed in plastic sample bags and shipped to Energy 
Laboratories in Billings, Montana. Samples for metal analyses 
were digested using hydrogen peroxide, EPA Method 3050 (EPA 
1982) . The very high iron content of most wetland soil samples 
necessitated special corrections to eliminate interference in 
samples analyzed by inductively coupled argon plasma spectro- 
scopy (ICP) . These corrections were calculated at the 25 
percent iron content and, as a result, were insufficient for the 
very high iron content found at site C-1. The apparent 
concentrations reported for several metals are depressed at site 
C-1 for this reason but the laboratory reported values are still 
included in Appendix C. Another result of the iron interference 
was the high tin and antimony values (positive interference) 
reported in the initial screening samples upon which the 
selection of parameters to be analyzed was based. As a 
consequence, tin and antimony were run on samples from several 
sites and, following interference corrections, all values were 
below the instrument detection limit. Two Natior\l Bureau of 
Standards (NBS 1982) samples of river sediment, t ;andard 1645, 
were run with wetland samples. Recoveries were vithin 10 
percent for all parameters except aluminum, arst- t Ic and iron. 
The recoveries were: aluminum, 28%; arsenic, 79%; cadmium, 90%; 
chromium, 97%; copper, 92%; iron, 79%; mercury, 100%; manganese, 
92%; nickel, 90%; lead, 95%; and zinc 94%. The low recoveries 
for aluminum make interpretation of these data difficult and the 
findings should, therefore, be considered only relative. No 
adjustment of these data has been made since, with the exception 
of aluminum, most elements were within or very near the reported 
range of estimated uncertainties for the standard. Chloride, 
fluoride, bicarbonate and sulfate were analyzed from a saturated 
paste extract. The cation exchange capacity was determined by 
Method 8-4, American Society of Agronomy Monograph 9 (1965) for 
acid soils. Organic carbon content was calculated from loss on 
ignition at 550°C for 2 hours. The data base consists of single 
samples for each selected interval at each site, generally based 
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on field observed stratigraphic changes. Fourteen laboratory 
duplicates and the two standards were analyzed for most 
parameters. 

Selected metal deposition in the wetland due to the Swamp 
Gulch discharge has been estimated utilizing 1) isopachs of the 
acrotelm (Figure 5-7), 2) isoline maps of chemical concentra- 
tions of specific metals, and 3) summation of volumes multiplied 
times concentrations multiplied times bulk density. Background 
concentrations were subtracted before isoline maps were con- 
structed for this purpose. 

8.2 RESULTS AND DISCUSSION 

In general , metal and other chemical parameters are more 
varied than similar data for ground-water which is likely due to 
a more heterogeneous nature of the soils. Occasional 
anomalously high or low values for several elements have been 
noted and may be due to sampling or laboratory ICP interference. 
Samples obviously in error have not been used in the following 
interpretation. Total masses of metals deposited in the 
acrotelm due to AMD have been calculated based on the 3.9 ha 
(9.6 ac) area found to be impacted by iron deposition. The 
total volume of the impacted acrotelm area is approximately 
23,000 cubic meters (30,000 yd^) . 

Most trace metal concentrations in the acrotelm show 
similar areal distribution. Aluminum, cadmium, copper, lead, 
manganese and zinc all exhibit the highest concentrations in the 
area immediately south of Montana Highway 200 and generally west 
of site C-1. Iron distribution more closely follows the 
preferential ground-water flow paths (Section 5) . The following 
paragraphs describe the distributions of specific metals in the 
wetland soil materials. 

8.2.1 Areal Distribution of Selected Elements 
8.2.1.1 Iron 

Iron concentrations in the wetland acrotelm vary from a 
maximum of 415,000 ppm at site C-1 to 36,400 ppm at site AA-2 
(Figure 8-1). The distribution pattern is similar to the 
ground-water flow paths found during the tracer study which 
strongly suggests the Swamp Gulch drainage is the source of most 
of the iron found in the acrotelm. Concentrations of iron fall 
from the maximum (about 8 times background levels) to values 
about 2 times above background levels at site F-1. This 
indicates the wetland has been effective in past removal of iron 
from the AMD. 
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Figure 8-1. Iron concentrations in the wetland acrotelm. 



An estimation of acrotelm iron mass due to the Swamp Gulch 
AMD discharge indicated approximately 550 metric tons (600 U.S. 
tons) of iron were deposited from this source. This amount 
suggests the present annual estimated iron loading [(340 kg/year 
(750 lbs/year - Section 7) has been greatly exceeded in the 
past. It is likely that iron discharges to the wetland were 
notably higher during the active mining/milling period and that 
the annual estimated amount may have been underestimated due to 
abnormal climatic conditions during the study period. 
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8.2.1.2 Manganese 

The distribution of manganese in shallow wetland soils is 
marked by high concentrations in the site D-1 to site F-1 area 
(Figure 8-2) . Values at site D-1 are approximately an order of 
magnitude above background sites AA-2 and AAA-3. Concentrations 
of soil manganese decrease with increasing distance from site D- 
1 suggesting some manganese is being removed from the ground- 
water and sequestered in the acrotelm. However, the manganese 
concentration at site F-1 is still three times greater than 
background levels which indicates a notable amount of manganese 
may likely be discharged from the wetland. Approximately one 
metric ton (2,400 lbs) of manganese has been deposited in the 
wetland acrotelm from the Swamp Gulch AMD. The boundary for the 
elevated manganese impacted area was less well defined than that 
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Figure 8-2. Manganese concentrations in the wetland acrotelm. 
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for iron and, therefore, the one ton estimate may be conserv- 
ative. However, it is clear that the efficiency of this wetland 
in the removal of manganese is at least an order of magnitude 
less than for iron. 

8.2.1.3 Aluminum 

The relationship between the Swamp Gulch discharge and 
elevated aluminum levels in the acrotelm are the least distinct 
of the metals evaluated, possibly due to ICP iron interference 
as suggested by the low recoveries found in the NBS standards. 
Relative concentrations at several sites (BC-1, C-2, D-1, E-1, 
E-2 and F-1) are apparently two to three times higher than 
background levels at AAA-3, AA-2 and A-3 (Figure 8-3). However, 
the reported aluminum values at sites CD-I, CD-1-2, and C-1.5, 
all close to the Swamp Gulch discharge point, were similar to 
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Figure 8-3. Aluminum concentrations in the wetland acrotelm. 
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background concnetrations . The relative aluminiim levels in the 
acrotelm are consistently higher to the west of the C-1 to C-5 
line than they are upgradient to the east. It is possible that 
iron interference may have produced the low aluminum values 
reported at sites CD-I, CD-1-2 and C-1. 5 and may have produced 
lesser uncorrected interference with all samples. Due to the 
apparent problems evident with aluminum analyses, no estimation 
of the aluminum mass present in the acrotelm has been attempted. 

8.2.1.4 Copper 

Copper levels in the wetland acrotelm exhibit a marked 
increase in concentration near the Swamp Gulch discharge point. 
Background concentrations as represented by sites AAA-3, AA-2 
and A-3 range from 210 mg/kg to 44 mg/kg. Sites B-2 and B-3 
also exhibit low copper concentrations at 19 and 40 mg/kg, 
respectively. Acrotelm copper values rise abruptly west of a 
line running due south from site BC-1, which roughly corresponds 
to the 100 mg/kg isoline (Figure 8-4) . The highest concen- 
trations are found in the C-1. 5 to F-1 zone where copper levels 
are generally 20 to 30 times the background mean. The remaining 
wetland area to the south of this zone is characterized by 
values 6 to 8 times the mean background concentration. The 
distribution pattern suggests some amelioration of copper has 
occurred in this system, but the high levels still evident at 
site F-1 suggest much copper may pass through the system and 
discharge to the Blackfoot River. The acrotelm copper mass is 
about four metric tons (4.4 U.S. tons) which represents 325 
years of deposition at the 1987-88 annual rate of 11.6 kg (25.5 
lbs). While this period would be well in excess of the period 
from mine initiation to the present, past loading rates were 
likely higher (as observed with iron and lead deposition) and 
hence, removal of copper from the AMD is likely considerably 
less than 100 percent. The 40 to 50 mg/kg range, typical of the 
wetland background areas, is common in numerous agricultural 
soils (EPA 1987b) . Other studies have found soil total copper 
levels of 500 to 700 ppm to produce severe yield reductions in 
rye grass and bush beans (Wallace et al. 1977, McGrath et al. 
1982) . Numerous other crop species are significantly impacted 
at soil total copper levels as low as 100 ppm (EPA 1987b) . 

8.2.1.5 Lead 

Lead concentrations in the wetland acrotelm vary from less 
than 1 mg/kg to values greater than 1000 mg/kg. Background 
values upgradient from the Swamp Gulch discharge point are 
generally less than 100 mg/kg but are higher than the 11.6 to 
32.2 mg/kg range that Chattopadhyay and Jervis (1974) found for 
the Holland Marsh muck soil. The anomalously high level (450 
rog/kg) found at background site AAA-3 is likely the result of 
other influences and not the result of the Swamp Gulch AMD 
discharge. The wetland system is apparently quite efficient in 
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Figure 8-4. Copper concentrations in the wetland acrotelm. 



the removal of lead, likely due in part to the low mobility of 
lead in soils. The Swamp Gulch acrotelm contains a AMD derived 
lead mass of 1.3 metric tons (1.4 U.S. tons) which represents 
2,350 years of deposition at the calculated 0.54 kg (1.2 lbs) 
annual loading rate. It is clear that past lead loading rates 
were higher than those observed during this study. While some 
of the observed lead may have been derived from natural 
mineralization in the area, it is doubtful that much lead from 
this source would have been deposited in the acrotelm due to the 
semi-isolation of the acrotelm provided by the more limited 
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hydraulic conductivity of the underlying catotelm and the 
limited mobility of lead. Elevated levels found in the zone 
roughly described by sites C-1, C-1.5, D-3 and E-1 are an order 
of magnitude above background sites. The concentrations at 
sites E-2 and F-l are approaching background levels. Of 
particular interest are the high lead concentrations found in 
near surface materials at sites F-2 and F-3 (Figure 8-5). These 
are likely the result of mining activities on the upper 
Blackfoot River above this site. A literature review on the 
phytotoxicity of selected metals (EPA 1987a) suggests lead 
concentrations observed in the elevated zone would be phytotoxic 
to most crops and may inhibit maximum vegetation production in 
the wetland. 
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Figure 8-5. Lead concentrations in the wetland acrotelm. 
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8.2.1.6 Zinc 

Zinc concentrations in the wetland acrotelm exhibit several 
anomalies, including high levels at background sites AAA-3 (1780 
mg/kg) and A-3 (2092 mg/kg) and a low concentration at site C-1. 
The cause of these variations is not apparent but may be due in 
part to analytical interference. It is evident that zinc 
concentrations in the acrotelm are consistently elevated south 
of Montana Highway 200 from sites C-1. 5 to F-1 (Figure 8-6). 
Zinc concentrations in this area are generally an order of 
magnitude higher than most background levels. It is also 
evident that zinc concentrations are not very effectively 
reduced down-gradient from the Swamp Gulch discharge point. 
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Figure 8-6. Zinc concentrations in the wetland acrotelm. 
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Concentrations of zinc in the elevated zone are consistently 
above 1500 mg/kg, a level that would be expected to result in 
nearly a 100 percent yield reduction for most agricultural crops 
(EPA 1987a) . 

The Swamp Gulch wetland acrotelm zinc mass due to AMD is 
approximately 6.0 metric tons (6.6 U.S. tons) which would 
represent 137 years of deposition at the observed 43.6 (95.9 
lbs) annual loading rate. Although this time interval is in 
excess of historic mining and AMD discharges, past discharges of 
zinc were undoubtedly higher and hence, the actual time interval 
represented is likely considerably less. 

8.2.1.7 Arsenic 

Arsenic concentrations at the four wetland sites analyzed 
for this element (AA-2, C-1, C-4 and F-1) exhibit similar 
distribution patterns with depth. All sites are characterized 
by highest concentrations located at the top of the muck zone. 
The maximum arsenic value is 86 mg/kg for the 0.46 to 0.91 m 
interval at site C-4. Values at AA-2 and C-1 are similar and 
suggest the wetland acrotelm has little effect on AMD arsenic 
content at this study site. Total arsenic levels as low as 10 
ppm have been shown to be detrimental to production for some 
crop species (EPA 1987a) . Surf icial acrotelm values range from 
18.5 mg/kg at site AA-2 to 29.5 mg/kg at site C-4, both of which 
sites are not likely impacted by the Swamp Gulch AMD. Some 
impact to wetland vegetation and wildlife may occur from the 
observed arsenic concentrations. 

8.2.1.8 Barium and Bervllium 

Barium and beryllium distributions are generally similar to 
arsenic. These two elements are more generally distributed 
within the muck zone with the highest concentrations tending 
more towards the central area of this stratigraphic zone. A 
notable exception to this observation occurred at site C-1 where 
the maximum beryllium concentration (240 mg/kg) occurred in the 
acrotelm, strongly suggesting this element is enriched due to 
the Swamp Gulch discharge. The maximum barium concentration 
(2140 mg/kg) was found in the 1.77 to 1.89 m depth interval at 
site F-1. 

8.2.1.9 Cadmium 

Cadmium concentrations in the wetland soils range from <0.5 
to 45 mg/kg. This metal is apparently enriched 10 to 20 times 
in wetland soils due to the Swamp Gulch AMD. The highest 
concentrations are associated with the acrotelm layer at site F- 
1 (44 and 45 mg/kg) . At site F-1, the muck zone maximum is 4 
mg/kg. This is in contrast to site C-1 where there is 
apparently little difference between most muck values and the 
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acrotelm (1.2 to 21 mg/kg and 22 mg/kg, respectively) . The AMD 
derived acrotelm cadmium content is only 4 kg (8.8 lbs) which 
represents 6.7 years of deposition at present loading rates. 
These data suggest some Swamp Gulch discharge is percolating 
through the muck zone at site C-1 and that the wetland is 
generally less effective in ameliorating cadmium in AMD. At 
background site AA-2, the cadmium concentration in the acrotelm 
is 2 mg/kg and a maximum concentration of 8.6 mg/kg is 
associated with the gravels beneath the muck zone. All muck 
zone values at this site are less than 2 mg/kg. The 22 to 45 
mg/kg values found in the acrotelm at sites C-1 and F-1 can be 
expected to severely impact production of many crops (EPA 1987a) 
and can be expected to add significant cadmium to the food 
chain, possibly impacting wetland wildlife (Section 9) . 

8.2.1.10 Cation exchange capacity 

Cation exchange capacity (CEC) of the acrotelm determined 
on sites AA-2, C-4 and F-1 (Appendix C) , had a weighted mean of 
35.6 raeq/lOOg. The range of CEC in the wetland, including site 
C-1, was 4.9 to 45.0 meq/lOOg. These CEC findings are lower 
than those typically found in wetlands. Mathur and Farnham 
(1985) found wetland peat materials with a mean CEC of 148.5 
meq/lOOg. These low CEC findings in the Swamp Gulch wetland may 
be the effect of the method of CEC analysis. Presuming the CEC 
found is true, the low CEC may be a reason for the limited level 
of wetland effectiveness in ameliorating the AMD metals. 

All samples from sites AA-2, C-1, C-4 and F-1 that have 
been analyzed for mercury, antimony, and thallium are below the 
1 mg/kg instrument detection limit for these elements. 

8.2.2 Distribution of Elements with Depth 

An evaluation of selected trace metals distributions with 

depth indicates that, in general, the anaerobic muck zone below 

the acrotelm contains lower concentrations of most elements 
(Table 8-1) . 

Aluminum is a notable exception, and is often found at 
higher concentrations in the muck zone. A comparison between 
mean background levels and concentrations within the AMD 
impacted area suggests the muck zone is impacted by AMD (Table 
8-1) . All metals evaluated in the impacted muck zone were 
higher than corresponding background values. However, due to 
the lower hydraulic conductivity of the muck zone, it likely has 
limited effective remediation potential for AMD amelioration. A 
comparison of muck zone concentrations with the underlying 
gravel materials reveals a less clear situation in which, 
compared to the catotelm, some metals are higher and some lower 
in the gravels, zinc levels were often higher in the gravels. 



99 



Table 8-1. Mean selected elemental concentration (mg/kg) 

changes with depth for background and AMD impacted 



areas, 



Site 



Al 



Cu 



Fe 



Mn Pb 



Zn 



Zone - Acrotelm 

Background 

(AA-2,A-3) 9,930 46 64,100 225 125 1380 

Impacted Zone 

(C-1,D-1,F-1) 27,150 1287 236,000 898 689 2580 



Zone - Catotelm 

Background 

{AA-2,A-3) 22,300 92 42,300 123 117 188 

Impacted Zone 

(C-1,D-1,F-1) 37,800 597 63,800 468 291 704 



Zone - Gravels 

Background 

(AA-2,A-3) 23,200 285 26,250 95 77 1190 

Impacted Zone 

(C-1,D-1,F-1) 17,830 563 32,870 413 170 693 



especially near the top of this unit (Appendix C) . Relative 
sulfate levels, as determined from the saturated paste extracts, 
suggest the distribution of sulfate is similar to zinc. The 
highest sulfate values for background sites AAA-3, AA-2 and A-3 
are found in the underlying gravels. This distribution is also 
evident at sites D-1 and D-2. Sites B-2, C-5 and D-3 have the 
highest sulfate concentrations in the acrotelm. The only site 
where high sulfate values are associated with the catotelm is C- 
1. Low sulfate values in the catotelm would be expected due to 
reducing conditions which would tend to reduce sulfate to 
insoluble sulfides which would not be recovered in a saturated 
paste extract. It is possible that the total sulfur 
concentration in this zone would exceed that in other strati- 
graphic units. This observation of sulfate removal is 
substantiated by the studies on microbiological sulfate removal 
presented in Section 10, Phase II. 
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8 . 3 SUMMARY 

Annual deposition of metals into the acrotelm have 
probably decreased from levels present during the period of 
active mining and milling. Acrotelm masses of lead, iron, 
copper and zinc would, at present depositional rates, represent 
accumulation periods ranging from 2,350 to 137 years (Table 8- 
2) . If the assumption is made that the apparent depositional 
interval (2,350 years) found for lead represents 100 percent 
removal of this metal, and that this interval can serve as a 
relative basis of comparison with other elements, the wetland 
efficiency in removal of metals is: iron, 70%; copper, 14%; 
zinc, 5.8%; manganese, 0.7%; and cadmium, 0.3% (table 8.2). 
Assuming past loading rates were higher than those observed at 
present, it is apparent that cadmium, manganese and likely zinc 
are all passed through the wetland without much amelioration. 
Copper may also be passed through the system, but it would 
appear that it is at least partially ameliorated. Both lead and 
iron may be effectively removed within the wetland area. Much 
of the apparent decrease in wetland water sample concentrations 
of cadmium, manganese, zinc and possibly copper with increasing 
distance from the AMD discharge point, is probably due to 
dilution processes rather than immobilization. The higher 
levels of metals observed in the gravels may be derived from two 
sources: 1) from the original gravel materials which contain 
pebbles and grains of mineralized materials derived from 
upstream outcrops; and 2) from mineralized ground-water flowing 
laterally in these materials. It is also 

possible, especially at site D-1, that there could be an upward 
ground-water gradient beneath the gravels in bedrock materials. 
This ground -water could be derived in part from seepage 
emanating from old underground mine workings. 

Table 8-2. Acrotelm metal masses and the relative efficiency of 
wetland metal removal. 







Annual 




Relative 




Acrotelm 


Input from 


Years represented efficiency 




metal mass 


AMD 


by deposition 


compared 




f metric tons) 


(Kq) 


rate 


to leadf%) 


Cadmium 


0.004 


0.54 


6.7 


0.3 


Copper 


3.8 


11.6 


325. 


14. 


Iron 


553. 


338. 


1,640. 


70. 


Lead 


1.3 


0.54 


2,350. 


100. 


Manganese 


1.1 


63.4 


16.9 


0.7 


Zinc 


6.0 


43.6 


137. 


5.8 
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It is clear that most observed sequestering of AMD metals 
occurs in the acrotelm, with the catotelm playing a minor role. 
Efforts should be concentrated on constructing effective acrotelm 
zones in artificial wetlands, and over time, a natural catotelm may 
form from material derived from the acrotelm. 
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9.0 WETLAND VEGETATION AND THE CONTROL OF ACID MINE DRAINAGE 

The objectives of this portion of the Swamp Gulch wetland 
study were to: 

1) quantitatively describe plant associations and determine 
production by dominant species, and 

2) assess metal loading of dominant vegetation and 
potential impacts to wildlife in the natural wetland. 

9.1 SAMPLE COLLECTION, PREPARATION AND ANALYSIS 

Vegetation observations were recorded during three trips to 
the Swamp Gulch wetland from April through August, 1987. Plant 
reference specimens and vegetation samples were collected from 
the Swamp Gulch study area and from the Hardscrabble Creek 
background wetland during August 1987. Specimens have been 
deposited in the Reclamation Research Herbarium. 

The composition of the vegetation at the Swamp Gulch site 
was determined by estimating plant canopy coverage. Eight 
transects were located across the study area, along which 
quadrats (20 X 50 cm) were placed at regular, preselected 
intervals. A total of 104 quadrats were used. Species specific 
canopy coverage was estimated at each sampling station to the 
nearest percent (Sensu Daubenmire 1959) . Depth to the water 
table was also estimated at each sampling station. At the 
background site on Hardscrabble Creek, canopy coverage values 
for the dominant species were visually estimated for the site as 
a whole. 

Aboveground production for the dominant herbaceous species 
(Carex rostrata ) was estimated at the Swamp Gulch and the 
Hardscrabble background site using 0.25 m^ (0.5 x 0.5 m) 
quadrats. Plants were clipped at ground level from seven 
quadrats at the Swamp Gulch site and in two quadrats at the 
background site. These samples were placed in paper bags, oven- 
dried at 50° C to a constant weight and then weighed to the 
nearest 0.1 gram. Values were converted to kg/ha and 
comparisons were made between the Swamp Gulch site and the 
Hardscrabble site at the 0.05 level using the Students t-test. 

Vegetation samples from throughout the Swamp Gulch area and 
background sites were collected for chemical analysis. To 
expedite this process, the Swamp Gulch study area was divided 
into seven collection zones. From each zone, whole plants were 
collected for Carex rostrata and the mosses, Isopterygium 
pulchellum and Sphagnum tenellum . Current growing season shoots 
and leaves were collected from dominant shrub species (Salix 
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boothii and Betula alandulosa ) , as were roots up to approxi- 
mately 3 cm in diameter. Aboveground materials were placed in 
paper bags and below ground materials were placed in plastic 
bags. Initial sample preparation consisted of washing each 
sample with tap and distilled water and then subjecting it to 
ultrasound until no residues were observed in the water. Each 
sample was then oven-dried (50° C) , placed in paper bags and, 
together with National Bureau of Standards (NBS) quality control 
samples, were shipped to Energy Laboratories in Billings, 
Montana on September 28, 1987. A total of 67 samples (65 
natural and 3 NBS citrus leaf samples) were analyzed for Al, As, 
Cd, Ca, Cu, Fe, Pb, Mn, Ni and Zn. At the laboratory, each 
sample was digested with nitric acid and hydrogen peroxide 
according to method 3050 of the EPA Manual SW-846 (EPA 1982). 

Concurrent with vegetation sampling, a map was drawn 
delimiting the plant associations, areas of open water, transect 
locations, and man-made features such as roads. 

Canopy coverage and frequency data were summarized for the 
Swamp Gulch study area as a whole. Nomenclature follows Dorn 
(1984) for vascular plants and Crum et al. (1973) for mosses. 
The final plant association map was constructed with the aid of 
the field map, field notes, canopy cover data, and after a 
review of wetland classification schemes made by other authors. 



9.2 DATA QUALITY ASSURANCE/ QUALITY CONTROL 

The quality of vegetation chemical data were determined for 
both field and laboratory procedures. Field sampling included 
the collection of duplicate samples and the insertion of NBS 
Blind Field Standard (BFS) reference samples into the sample 
train. Laboratory QA/QC procedures consisted of splits of 
natural and field duplicate samples. Analytical results of 
these samples were used to calculate precision and accuracy 
statements following standard procedures developed for the EPA 
Contract Laboratory Program. Appendix D-1 presents the results 
of the data accuracy and precision statements. 

Accuracy is interpreted through the following example: for 
Al, we are 90% confident that the reported value is 44.5 ± 6.3% 
of the true value. Precision is interpreted as follows: for 
Al, we are 90% confident that the true value is within ±4.5% of 
the reported values. 
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9.3. RESULTS AND DISCUSSION 

9.3.1. Vegetation Patterns 

The Swamp Gulch wetland was a mosaic of plant 
associations, the presence of which was directly related to 
flooding frequency and duration. Land surface varied from areas 
of continually open water (about 4.6% of the site) to areas 
where the soil was periodically saturated (sometimes nearly year 
round) , to areas where the water table was rarely high enough to 
saturate the entire soil profile. The wettest areas were 
characterized by pure stands of Carex rostrata (Figure 9-1). 
This species was nearly ubiquitously present, being associated 
with Betula qlandulosa and Salix boothii, and with Pinus 
contorta. Carex rostrata was found growing as an emergent (in 
as much as 20 cm of water) and where the soil surface was only 
damp. A relatively dry site was occupied by Pinus contorta, 
Chimaphila umbellata and S phagnum tenellum . The driest areas 
were dominated by species indicative of the surrounding forest; 
these areas were characterized by Pinus contorta, Picea 
englemannii . and Arctostaphylos uva-ursi . A list of vascular 
plants and bryophytes identified within the Swamp Gulch wetland 
study area can be found in Appendix D-2. Beaver dams, extending 
across the study area near transect 6 (Figure 9-1) , contributed 
to the patterning of plant associations. The eastern two-thirds 
of the site had a higher water table than the western one-third, 
and was therefore dominated by Carex rostrata and associated 
communities. The western portion was dominated by Pinus 
contorta and forest grasses and forbs. 

9.3.2. Canopy Coverage and Production 

Carex rostrata was the dominant species throughout the 
study area. It was found at over 90 percent of the sampling 
stations and covered more than 55 percent of the study site 
(Table 9-1) . Betula glandulosa was the dominant shrub and Pinus 
contorta was the dominant tree. Salix boothii was also 
important; it covered 5.5 percent of the study area and was 
found at almost 20 percent of the sampling stations. Two 
mosses. Sphagnum tenellum and Isopterygium pulchellum were 
commonly encountered, especially in the wetter, eastern two- 
thirds of the site. Chimaphila umbellata was associated almost 
exclusively with Pinus contorta in a small area adjacent to 
Highway 200 (Figure 9-1) . Sphagnum tenellum grew profusely 
within this association. 

Many dead Pinus contorta can be observed, especially in 
the eastern portion of the wetland. These were from 20-30 feet 
tall and probably from 40-60+ years old. The most likely cause 
of death was a rise in the water table, which may have been the 
result of fluctuating climatic conditions. 
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Table 9-1. 



Mean percent canopy coverage and percent frequency 
for plant species. 





Life 


% Canopy 


% 


Species 


Forml 


Coveraae 
55.5 


Freauencv 


Carex rostrata 


GL 


90.4 


Betula glandulosa var. hallii 


S 


11.3 


31.7 


Pinus contorta 


T 


6.3 


8.7 


Salix boothii 


S 


5.5 


19.2 


Sphagnum tenellum 


M 


4.7 


12.4 


Chimaphila umbellata 


S 


3.0 


8.7 


Potomogeton sp. 


F 


2.8 


1.9 


Isopterygium pulchellum 


M 


2.8 


18.1 


Pyrola asari folia 


F 


1.2 


7.6 


Arctostaphylos uva-ursi 


HS 


1.1 


5.7 


Abies lasiocarpa 


T 


0.9 


1.0 


Juniperus communis 


S 


0.9 


1.9 


Picea engelmannii 


T 


0.9 


4.8 


Scirpus acutus 


GL 


0.9 


1.2 


Pseudotsuga menziesii 


T 


0.6 


1.0 


Equisetum sp. 


F 


0.5 


11.5 


Vaccinium scoparium 


HS 


0.3 


1.0 


Achillea millefolium 


F 


0.1 


1.9 


Elymus glaucus 


G 


<0.l2 


1.0 


Fragelia virginiana 


F 


<0.1 


4.8 


Festuca idahoensis 


G 


<0.1 


1.0 


Geum macrophyllum 


F 


<0.1 


1.0 


Phleum pratense 


G 


<0.1 


1.0 


Unknowns 




2.2 


1.0 




TOTAL 


101.8 





LF 



Life Form: GL = Grasslike, S = Shrub, T = Tree, M = Moss, 



F = Forb, HS = Half -Shrub and G = Grass 
2 Values of <0.1 were tabulated as 0.05. 



Herbaceous production on the Swamp Gulch wetland averaged 
3820 kg/ha (Table 9-2) . For comparison, vegetation composition 
and production were estimated at a pristine background wetland 
on Hardscrabble Creek (off the Alice Creek road) . The 
composition of the vegetation was very similar to the Swamp 
Gulch wetland. Carex rostrata dominated with about 60 percent 
canopy coverage and Salix boothii and Betula glandulosa were the 
dominant shrubs. Water temperature, pH and EC were 11° C, 7.4 
and 400 umbos/ cm, respectively. Herbaceous production at the 
background site was 3750 hg/ha, which was not significantly 
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Table 9-2. Peak standing crop (kg/ha) at the Swamp Gulch 
and Hardscrabble Creek (background) wetlands. 

Hardscrabble Creek Swamp Gulch Wetland 

(background) 
Mean^ STD Range of Values Mean STD Range of Values 

3750a ±532 3374-4127 3820a ±1550 1983 - 6776 

^ Means were not significantly different at the p<.05 level. 

different than production on the Swamp Gulch site (Table 9-2) . 

9.3.3 Element Concentrations in Plant Material 

This section presents data on the concentration of elements 
in the dominant plant species at the Swamp Gulch study site and 
at the Hardscrabble Creek (background) wetlands. These data are 
compared, contrasted and discussed with respect to water quality 
data. 

Statistical analysis was conducted to determine if 
elemental differences in plant tissues could be detected between 
the Swamp Gulch and the background sites. Large variations in 
the analytical data (and low sample numbers) resulted in large 
standard deviations. This situation restricted the usefulness of 
the statistical analyses by limiting the number of possible 
significant determinations. As a consequence, element data are 
discussed in a general manner without the use of statistics. 

A literature search revealed no data on the metal levels of 
the plant species occupying the study site. Considerable data 
exists however, for species of cattails (T ypha ) . Hutchinson 
(1975) reviewed the literature on metal levels in aquatic 
vegetation. EPA (1987a and b) conducted a literature review on 
metal concentrations in terrestrial macrophytes for the purpose 
of establishing plant and animal hazard criteria. These data 
sets will be compared, in a general way, with element data in 
the tissues of the macrophytes at the Swamp Gulch and 
Hardscrabble Creek wetlands. Metal levels in a variety of 
bryophytes, particularly Sphagnum spp., growing in metal- 
enriched areas have been reported. These data will also be 
compared with the bryophytes metal data from the Swamp Gulch and 
Hardscrabble sites. 

Water quality data are presented for three locations at the 
Swamp Gulch site. This is done to show the degree of metal 
enrichment in water as it flows past the mine site, and to help 
explain the level of metal attenuation by the wetland. 
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Background water chemistry levels were established from samples 
collected in Swamp Gulch above the Carbonate Mine. References 
will also be made to the water quality at station F-1, located 
directly down-gradient from the Swamp Gulch wetland. These data, 
therefore, should provide information as to the quantity of 
metals being attenuated by the wetland. 



Table 9-3. Concentration (ug/g) of elements in Carex rostrata . 



Background Wetland 
Element Matrix^ Mean(n=4) Range 



Swamp Gulch Wetland 
Mean(n=5V Range 



Al 



As 



Cd 



Ca 



Cu 



Fe 



Pb 



Mn 



Ni 



Zn 



AG 


61 


11-170 


232 


41-400 


BG 


408 


72-1085 


1720 


850-2670 


AG 


0.9 


0.352-2.1 


0.9 


0.6-1.4 


BG 


16 


1.8-35 


19 


1.5-41 


AG 


0.8 


0.35-1.6 


0.9 


0.5-3.5 


BG 


6.2 


0.6-21 


14 


6.5-17 


AG 


3285 


360-4530 


2456 


1110-3450 


BG 


6420 


2340-15000 


3500 


670-5150 


AG 


3.8 


0.352-8.7 


62 


0.35-160 


BG 


30 


5.1-78 


984 


88-2240 


AG 


1327 


350-3840 


6134 


2440-15300 


BG 


20800 


5070-58500 


37140 


24800-54000 


AG 


6.5 


1.0-22 


29 


5.2-70 


BG 


128 


2.0-490 


354 


190-800 


AG 


248 


81-460 


250 


160-430 


BG 


177 


88-380 


131 


43-190 


AG 


0.72 


0.7-0.7 


0.8 


0.7-1.2 


BG 


1.9 


0.7-5.5 


9.3 


0.7-17.0 


AG 


96 


14-320 


134 


69-290 


BG 


1096 


42-4210 


1296 


310-2700 



Ipiant material: AG=above-ground, BG=below-ground 
2 Values below the instrument detection limit were estimated by 
multiplying the instrument detection limit value by 0.7. 
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9.3.3.1 Effects of Carex rostrata in Remediating AMD 

Below-ground tissues had higher concentrations than above- 
ground tissues (Table 9-3). This was true for the background 
and the Swamp Gulch wetlands and for all the elements with the 
exception of Mn, which appeared to be more concentrated in the 
leaves than in the roots/rhizomes. These findings were consis- 
tent with previous research on metal uptake by species of Typ ha 
(Bayly and O'Neill 1972, Heil and Kerins 1988, Taylor and 

Crowder 1983a, 1983b) . In the present study, above-ground 
tissues concentration were different from below-ground 
concentrations by an order of magnitude for Al, As, Cd, Cu, Fe, 
Pb and Zn. 

Water quality samples collected down-gradient from the 
Swamp Gulch wetland indicated that some of the metals were being 
attenuated by Carex rostrata . Water at site F-1 had lower 
concentrations of Al, Cd, Cu, Fe, Pb, Ni and Zn than the 
influent water (Table 9-4) . Sulfate was also lower and pH was 
unchanged. 

Some elements (Al, Cu, Fe and Pb) were substantially 
elevated in plant tissues from the Swamp Gulch site, compared to 
background levels. These elevated levels were apparently not 

Table 9-4. Water quality^ in Swamp Gulch above and below the 
Carbonate mine site, and below the Swamp Gulch 
wetland at station F-1. 





Above Carbonate 


Below Carbonate 


Swamp Gulch 




Mine 


Mine 


Wetland 


Parameter 


(background) 


rAMD) 


rstation F-1) 


Al 


<0.1 


6.2 


0.4 


As 


<0.005 


<0.005 


<0.005 


Ca 


12.7 


22.0 


26.0 


Cd 


0.001 


0.038 


<0.001 


Cu 


<0.01 


0.61 


0.04 


Fe 


0.7 


25.8 


15.1 


Pb 


<0.01 


0.035 


<0.01 


Mn 


0.08 


3.73 


4.30 


Ni 


<0.03 


0.37 


<0.03 


Zn 


<0.01 


2.64 


0.12 


SO4 


10.7 


202.3 


142.0 


pH 


6.8 


3.6 


3.3 



^Values are means from three sampling dates: 4/29/87, 7/28/87 
and 1/11/88. 
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phytotoxic . Carex rostrata did not appear to be under abnormal 
stress at the Swamp Gulch site; production (3820 kg/ha) and 
cover (56 percent) were high, and individual plants appeared to 
be healthy. As previously discussed, production was 
statistically similar between the two sites and no difference in 
canopy cover was visually apparent. Phytotoxicity for Carex 
rostrata is discussed by element in the remainder of this 
section. 

Although Al was elevated in plant tissue (approximately four 
times background) at the Swamp Gulch site, the above-ground mean 
of 232 ug/g fell within the range for aquatic plants reported by 
Hutchinson 1975 (Table 9-5) . The below-ground tissue concentra- 
tion of 1720 ug/g , however , was substantially higher than data 
from non-enriched areas reported by other researchers (Table 9- 
5) . Heil and Kerins (1988) reported similar results for T ypha 
latifolia in a wetland constructed in Montana to control AMD; 
tissue concentrations of Al were higher from the area receiving 
AMD than from a control site. The uptake of Al by Carex 
rostrata may have contributed to the low quantity (0.4 mg/L) of 
Al measured in water flowing from the Swamp Gulch wetland (Table 
9-4) at station F-1. The most reasonable explanation for the 
higher levels of tissue Al in the Swamp Gulch wetland was that 
Carex rostrata has the ability to absorb this element if it is 
available. Water quality data showed that Al was substantially 
higher in the AMD entering the 



Table 9-5. 



Range and mean element concentrations (ug/g) in 
aquatic forbs and grasses (from Hutchinson 1975)'^ 



Metal 



Range 



Mean 



Al 
As 
Cd 
Cu 
Fe 
Pb 
Mn 
Ni 
Zn 



250-785 
2.8-202 
2.6-28 
2.5-243 

70-31, 500^ 
2.0-53 
100-23,000 
1.1-44 
26.5-1000 



366 

8.0 
48 
3170 
11 
2380 

143 



-"■Except where noted, these values are from non-enriched environ- 
ments. However, some of the plant species investigated may be 
metal lophytes . 

^From As-enriched environment. 

•^This range may be biased high due to analyses of vegetation that 
were not completely cleaned of Fe precipitates. 
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wetland compared to water above the Carbonate Mine site (Table 9- 
4) . This indicates, at least to some degree, that Al will be 
absorbed by Carex rostrata . 

The average concentration of Fe in the influent AMD (25.8 
mg/L) was 37 times higher than the water above the Carbonate Mine 
site (0.7 mg/L). Although these data are limited, the concentra- 
tion of Fe in water at station F-1 was 15.1 mg/L, which was less 
than the concentration in the influent water (25.8 mg/L). Tissue 
concentrations for Fe in the Swamp Gulch wetland were 4.6 (above- 
ground) and 1.8 (below-ground) times higher than the background. 
This suggests that a portion of the dissolved Fe was being 
absorbed by Carex rostrata . Iron precipitates adhering to 
organic matter also contribute to a reduction in dissolved Fe 
content as the water passes through the wetland. Leaf tissue 
levels of Fe at the Swamp Gulch wetland were within the range 
reported for aquatic plants (Table 9-5). Bayly and O'Neill 
(1972) reported that Tvpha latifolia had leaf tissue Fe levels 
from 200 to 1800 ug/g and rhizome Fe levels from 700 to 3600 
ug/g, depending on the season in which samples were collected. 
In a metal smeltering area, Fe concentrations in T ypha latifolia 
leaves and roots reached 905 and 57,138 ug/g, respectively 
(Taylor and Crowder 1983b) . These data were comparable to those 
measured in the Swamp Gulch wetland, as were values reported by 
Heil and Kerins (1988) . These researchers reported Fe concentra- 
tions in leaf tissue of 286 ug/g and root tissue levels of 40,884 
for Typha latifolia growing in an artificial wetland receiving 
AMD. 

Like Al and Fe, the concentration of Zn in water increased 
as it flowed past the Carbonate Mine site (Table 9-4) and was 
likewise attenuated as the water passed through the wetland. The 
influent concentration of Zn was 2.64 mg/L, while the concentra- 
tion at station F-1 was only 0.12 mg/L. This large decrease was 
not due to absorption by Carex rostrata since tissue concentra- 
tions were only slightly greater at the Swamp Gulch wetland 
compared to the background wetland (Table 9-3) . Discussion in 
the following sections will show that the shrubs accumulated 
substantial amounts of Zn. Tissue levels in Carex rostrata from 
both sites were similar to aquatic plant data from non-enriched 
areas (Table 9-5) . In a metal-enriched environment, root and 
leaf Zn concentrations in T ypha latifolia ranged from 13 to 101 
ug/g and from 24 to 572 ug/g, respectively (Taylor and Crowder 
1983b) . These researchers also reported that tissue levels of Zn 
were not correlated with soil/sediment levels. Phytotoxicity of 
Zn in terrestrial plants varies from 60 to more 800 ug/g 
depending on the plant species, the variety and the part of the 
plant analyzed (EPA 1987a) . These authors suggested a 50 ug/g 
tolerable Zn level and a 500 ug/g phytotoxic level for above- 
ground tissue. 
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Dissolved Cu increased in Swamp Gulch water as it flowed 
past the Carbonate Mine site. Copper concentrations in water 
increased from <0.01 mg/L above the mine to 0.61 ug/L below the 
mine. Apparently, Cu was attenuated in the wetland since the 
concentration in water at station F-1 was 0.04 ug/L. Vegetation 
data indicated substantially higher (16 and 33 times greater than 
background in above-ground and below-ground tissues, respective- 
ly) Cu concentrations in Carex rostrata tissue from the Swamp 
Gulch wetland compared to the background site, suggesting that 
this species was absorbing dissolved Cu from AMD. Leaf tissue 
concentrations in the Swamp Gulch wetland (0,35 -160 ug/g) were 
within the range reported for aquatic plants (Table 9-5) . 
However, the range of Cu values in the below-ground material (88- 
2240 ug/g) often exceeded the range for other aquatic plants. 
From a natural but metal-enriched wetland, Typha latifolia leaf 
concentrations reached 24 ug/g and root levels reached 265 ug/g 
(Taylor and Crowder 1983a) . The global diagnostic background 
level for Cu in above-ground tissue of terrestrial plants ranges 
up to 20 ug/g (EPA 1987b) . Phytotoxicity in terrestrial plants 
generally begins at levels >20 ug/g. Some Carex rostrata plants 
in the Swamp Gulch wetland had leaf tissue levels up to 160 ug/g, 
suggesting that some phytotoxicity may be occurring due to Cu. 
However, Carex rostrata may have an internal exclusionary 
mechanism that allows it to tolerate relatively high levels of Cu 
in its tissues. 

The concentration of Pb in the wetland influent AMD was 
0.035 mg/L, which was slightly greater than <0.010 mg/L in the 
water at station F-1 (Table 9-4) . These data suggest that some 
Pb (albeit, very little) was being attenuated within the wetland, 
a portion of which may have been absorbed by the sedge. Carex 
rostrata Pb tissue levels were 4.5 (above-ground tissue) and 2.8 
(below-ground tissue) times higher in samples from the Swamp 
Gulch wetland compared to the background wetland (Table 9-3). 
Above-ground tissue concentrations were within the range for 
aquatic plants reported by Hutchinson (1975) (Table 9-5). It has 
been difficult to establish phytotoxic levels for Pb in plants 
because of the great variability between species. Background 
levels in lettuce leaves, for instance, have been reported as 
high as 50 ug/g, while phytotoxic levels in alfalfa have ranged 
down to 10.8 ug/g (EPA 1987a). These authors felt that a leaf Pb 
level of 25 ug/g could, in general, be tolerated by grain crops 
in Montana. Since leaf tissue levels ranged from 5.2 to 70 ug/g 
in the Swamp Gulch wetland, some plants may be under stress due 
to high levels of Pb. However, it is also possible that Carex 
rostrata has a mechanism which allows it to tolerate higher Pb 
levels than crop species. 

The lower concentrations of the aforementioned metals (Al, 
Cu, Fe, Pb and Zn) in water at station F-1 compared to the 
influent AMD indicates that a portion of these elements were 
being attenuated within the wetland. Vegetation data suggests 
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that Carex rostrata absorbed a portion of these metals. 
Manganese, by contrast, was not accumulated by Carex rostrata 
even though the concentration of Mn was notably elevated (47 
times) in the influent AMD. Concentration in the influent AMD 
was 3.73 mg/L compared to the background water level of 0.08 
mg/L. Both above- and below-ground tissue concentrations were 
similar between the background and Swamp Gulch wetlands. 
Apparently, Carex rostrata has an external exclusionary mechanism 
that prevents the absorption of Mn. Not only was Mn not being 
absorbed by Carex rostrata, it was not being attenuated in any 
other manner. The concentration of Mn in water down-gradient 
from the wetland (at station F-1) was 4.30 mg/L, which was 
slightly greater than the influent concentration of 3.73 mg/L. 
In contrast, Typha latifolia can absorb substantial amounts of Mn 
from water, Heil and Kerins (1988). These investigators reported 
mean leaf tissue levels of 1053 and 1760 ug/g of Mn in T ypha 
latifolia from two artificial wetlands receiving AMD having an 
average of 1.35 mg/L of Mn in the influent water. Manganese 
tissue levels have ranged up to 23,000 ug/g in aquatic plants 
(Table 9-5). However, these data were from a study conducted in 
the Ukraine and it is not clear from Hutchinson's discussion 
whether the samples were from a metal-enriched area or from a 
metallophytic species. 

The concentration of As in the tissues of Carex rostrata 
were similar between the background and the Swamp Gulch wetlands 
(Table 9-3) . The concentration of As in the Swamp Gulch water 
did not increase as it flowed past the Carbonate Mine site (Table 
9-4) . The tissue concentrations observed in the Swamp Gulch 
wetland were within the range reported for aquatic plants (Table 
9-5) . Depending on the plant species and the particular tissue 
analyzed, phytotoxic symptoms in terrestrial plants have been 
reported from 5 to 20 ug/g As (EPA 1987a) . Most of these data 
were from crop species, which may be more sensitive to As than 
Carex rostrata . Based on the currently available aquatic and 
terrestrial plant data, as well as data from the Hardscrabble 
Creek (background site) , there is no reason to suspect that As is 
toxic to Carex rostrata at the levels observed. 

The concentration of Cd in water above the Carbonate Mine 
was 0.001 ug/g compared to 0.038 ug/g below the mine site (Table 
9-4) . Mean leaf tissue concentrations were similar between sites 
but the range of values suggests that Carex rostrata was 
absorbing and accumulating more Cd at the Swamp Gulch wetland 
than the background wetland. Leaf tissue levels at the back- 
ground site ranged from 0.35 to 1.60 ug/g, compared to 0.5 to 3.5 
ug/g at the Swamp Gulch wetland. Tissue levels from this study 
were similar to those of other aquatic plants (Table 9-5). A 
large body of data exists on Cd toxicity in crops. Background 
levels for Cd have been reported up to 3.1 ug/g in alfalfa tops 
and yield increases have been reported when Cd was present at 
much higher concentrations (EPA 1987a) . In a study of hazardous 

114 



levels in crops from a metal enriched area, 10 ug/g Cd in plants 
was considered tolerable while a tissue level of 50 ug/g would 
definitely cause a yield reduction. In light of these hazard 
levels and the data from the Swamp Gulch site, the levels of Cd 
in Carex rostrata appear to be normal and well below the 
phytotoxic level. It should be noted that while these levels may 
not be phytotoxic, they are at least an order of magnitude higher 
than concentrations found in the AMD water and represent notable 
bioaccumulation . 

Like most of the other elements, Ca was more concentrated in 
the roots/rhizomes than in the leaves (Table 9-3). As previously 
mentioned, Ca levels were higher in tissues from the background 
wetland compared to the Swamp Gulch wetland. This was most 
likely due to a greater availability of Ca at the background 
site, although the water and soil at that site were not analyzed. 
Tissue Ca levels from the Swamp Gulch wetland were within the 
ranges reported by Taylor and Crowder (1983b) for Typha . Bayly 
and O'Neill (1972) reported similar levels for rhizomes (2000- 
9000 ug/g) but higher values for leaves (5000-11,000 ug/g), 
compared to Swamp Gulch wetland values. 

The concentration of Ni in leaf and below ground material in 
the Swamp Gulch wetland were probably not elevated above 
background (Table 9-3). These levels fell within Hutchinson's 
(1975) range for aquatic plants (Table 9-5). Taylor and Crowder 
(1983b) reported leaf tip concentrations for Tvpha latifolia that 
ranged up to 91 ug/g and root levels ranging up to 388 ug/g in 
plants from a metal enriched area. Nickel leaf concentrations of 
467 ug/g have been obtained without toxic symptoms, indicating an 
internal exclusionary mechanism (Taylor and Crowder 1983a) . 
Below-ground tissue concentrations have been correlated with 
soil/sediment concentrations, but no correlation between leaf and 
soil levels have been found (Taylor and Crowder 1983b) . 

9.3.3.2 Effects of Salix boothii in Rem ediating AMD 

Like Carex rostrata the metal concentrations in Salix 
boothii were generally higher in below-ground compared to above- 
ground plant tissue (Table 9-6) . This pattern occurred at both 
the background and Swamp Gulch wetlands and for all elements, 
with the possible exception of Zn. Zinc appeared to be more 
concentrated in the leaves/stems than the roots/rhizomes, 
especially in the Swamp Gulch wetland. 

Aluminum was absorbed by Salix boothii and concentrated in 
the below-ground tissue as opposed to the leaves and stems. The 
level of Al in the below-ground tissue from the Swamp Gulch 
wetland was substantially higher (3.4 times) than the tissue 
concentration from the background site (Table 9-6) . However, 
above-ground tissue concentrations were similar between wetlands. 
The ability of Salix boothii to absorb and store Al may have 
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Table 9-6. Concentration (ug/g) of elements in Salix boothii . 



Background Wetland 
(n=2) 
Element Matrix^ Mean Range 



Al 



As 



Cd 



Ca 



Cu 



Fe 



Pb 



Mn 



Ni 



Zn 



Swamp Gulch Wetland 
(n=6) 
Mean Range 



AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 



51 
190 

0.352 
2.9 

0.6 
0.6 

6680 
8580 

0.352 
11 

1535 
3000 

0.72 
3.9 

110 
125 

0.72 
0.7 

130 
105 



34-68 


50 


110-270 


655 


0.35-0.35 


0.35 


2.3-3.4 


3.6 


0.5-0.6 


5.6 


0.5-0.7 


9.1 


6560-6810 


6817 


7590-9570 


4535 


0.35-0.35 


1.3 


8.7-13 


169 


890-2180 


132 


1830-4290 


18418 


0.7-0.7 


1.7 


1.7-6.0 


112 


99-120 


257 


120-130 


456 


0.7-0.7 


0.9 


0.7-0.7 


1.6 


110-150 


388 


100-110 


278 



35-74 
170-1640 

0.35-0.35 
0.35-10.0 

1.0-9.9 
2.3-17 

5040-7750 
3240-4860 

0.35-6.0 
0.35-390 

59-420 
4470-33000 

0.7-3.8 
9.2-350 

150-500 
63-1380 

0.7-1.6 
0.7-2.7 

210-500 
160-420 



^Plant material: AG=above-ground, BG=below-ground 

2values below the instrument detection limit were estimated by multiplying 
the instrument detection limit value by 0.7. 

contributed to the decrease in Al measured in water leaving the 
Swamp Gulch wetland (Table 9-4) . 

Arsenic tissue concentrations were similar between the 
background and the Swamp Gulch wetland. This was expected since 
As was not elevated in the influent AMD. In the above-ground 
tissue, As was well below the phytotoxic range (5-20 ug/g) 
reported for terrestrial plants (EPA 1987a) . Arsenic in the 
below-ground tissues ranged from 2.9 to 3.6 ug/g, indicating 
that As accumulates in these materials over time. 
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Salix booth ii tended to accumulate Cd, despite the fact 
that it was present at a low concentration in the influent water 
(Table 9-4) . The level of Cd in the AMD entering the wetland 
was 0.030 mg/L compared to 0.001 mg/L in the non-enriched water 
above the Carbonate Mine site. This increase in dissolved Cd 
was enough to increase tissue levels more than nine times 
background (Table 9-6). The absorbance of Cd by Salix boothii 
may have contributed to the reduced level of Cd measured in 
water leaving the Swamp gulch wetland at station F-1. Based on 
the tolerable level (10 ug/g) for terrestrial plants (EPA 
1987a) , it is reasonable to assume that Salix boothii is not 
experiencing toxic affects from Cd. 

Calcium behaved differently in Salix boothii than in Carex 
rostrata . For Carex rostrata . Ca levels were lower in plant 
tissues from the Swamp Gulch site compared to the background 
site. For Salix boothii . this pattern only occurred for the 
below-ground tissues. The concentration of Ca was 8580 ug/g in 
the roots at the background site compared with 4535 ug/g at the 
Swamp Gulch wetland. Calcium levels in the above-ground tissues 
were similar at the background (6680 ug/g) and the Swamp Gulch 
(6817 ug/g) wetlands, respectively. 

Copper tissue concentrations were higher in above- than 
below-ground material and were also higher at the Swamp Gulch 
wetland compared to the background wetland (Table 9-6) . The 
concentration of Cu was 3.6 (above-ground tissue) and 15.3 
(below-ground) times greater in tissues from the Swamp Gulch 
wetland than the background wetland. Salix boothii probably 
contributed to the removal of dissolved Cu from the AMD as it 
flowed through the wetland. Above-ground tissue levels for 
Salix boothii were well below the phytotoxic level of 20 ug/g 
for terrestrial plants (EPA 1987b) . 

Iron concentrations were generally lower in the tissues of 
Salix boothii than in Carex rostrata . In the above-ground 
tissues for instance, the Fe level in the Swamp Gulch wetland 
was 132 ug/g, which was notably lower than the 1532 ug/g in 
tissue from the background wetland on Hardscrabble Creek. Iron 
concentrations in below-ground tissues were higher in the Swamp 
Gulch wetland compared to background tissue levels. 

Lead tissue levels were lower in Salix boothii than in 
Carex rostrata. The Pb was more concentrated in below- than 
above-ground tissue and more concentrated in tissues from the 
Swamp Gulch wetland than the background site. Compared to 
phytotoxic values for terrestrial plants (>10 ug/g) , it is 
unlikely that Pb caused toxic effects on Salix boothii (EPA 
1987a) . 
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Manganese was also distinctly more concentrated in the 
below- than the above-ground tissue, and was more concentrated 
in tissue from the Swamp Gulch wetland compared to the back- 
ground site. These patterns differ from those of Carex rostrata 
where Mn levels were highest in the leaves and similar between 
sites. The data suggests that Salix booth ii may have a greater 
affinity of Mn than Carex rostrata . Despite the absorption of 
Mn by Salix boothii . the concentration of Mn in the water was 
not reduced. The influent water value was 3.74 mg/L and the 
value at station F-1 was 4.30 mg/L (Table 9-4). 

Salix boothii did not appear to accumulate Ni to a notable 
degree. Tissue levels at the background site were below the 
instrument detection limit of 1.0 ug/g, so values have been 
estimated (Table 9-6) . Tissue concentrations at the Swamp Gulch 
site were not phytotoxic. 

Salix boothii responded to higher Zn levels in AMD at the 
Swamp Gulch wetland by absorbing more. Tissue levels of Zn were 
higher at the Swamp Gulch wetland than the background site. 
Salix boothii tended to store this element in leaves and stems 
rather than in roots. At the background wetland, Zn levels were 
130 ug/g in the above-ground tissue and 105 ug/g in the below- 
ground tissue. Likewise, Zn levels at the Swamp Gulch wetland 
were 388 and 278 ug/g in the above- compared to below-ground 
tissues, respectively. The Zn levels in these wetlands were 
well below the 500 ug/g phytotoxic level for terrestrial plants 
recommended for the Helena Valley (EPA 1987a) . 

9'3.3.3 Effects of Betula alandulosa in Remediating AMD 

Betula qlandulosa tended to absorb and accumulate metals in 
much the same way as Salix boothii. In general, tissue con- 
centrations for the analyzed metals were higher in the below- 
ground than in the above-ground plant material. This was true 
for both wetlands studied and for all elements, with the 
exception of Mn and Zn (Table 9-7). These patterns were similar 
for the accumulation of Zn in Salix boothii and the accumulation 
of Mn in Carex rostrata . 

With the exception of Al and Ca, plant tissue metal 
concentrations were higher in the Swamp Gulch wetland than at 
the background site. Aluminum levels were similar between the 
wetlands, while Ca values were lowest in samples from the Swamp 
Gulch wetland. The relatively low tissue metal levels make 
phytotoxicity of Betula alandulosa unlikely. 

9.3.3.4 Effect s of Brvophvtes in Remediating AMD 

Metal levels were determined for two mosses, Isopteryqium 
pulchellum and Sphagnum tenellum . Isoptervaium pulchellum was 
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Table 9-7. Concentration (ug/g) of elements in Betula glandulosa. 



Element Matrix^ 



Al 



As 



Cd 



Ca 



Cu 



Fe 



Pb 



Mn 



Ni 



Zn 



Background 
Mean fn=2) 



AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 

AG 
BG 



35 
141 

0.352 
0.35 

0.352 
0.35 

6835 
2150 

0.352 
7.1 

73 
205 

0.7 
5.5 

81 
39 

0.72 
0.7 

185 
56 



Wetland 
Range 



Swamp Gulch 
Mean fn=6) 



Wetland 
Range 



29-41 


28 


12-270 


151 


0.35-0.35 


0.6 


0.35-0.35 


0.7 


0.35-0.35 


0.5 


0.35-0.35 


1.4 


5670-8000 


3470 


1600-2700 


1808 


0.35-0.35 


2.8 


0.35-14 


15 


57-89 


158 


120-290 


1417 


0.7-0.7 


1.4 


5.1-26 


19 


76-86 


452 


11-67 


186 


0.7-0.7 


1.1 


0.7-0.7 


0.7 


160-210 


229 


30-82 


147 



23-31 
67-250 

0.35-1.6 
0.35-1.6 

0.35-0.7 
0.5-2.0 

2730-4320 
1070-3060 

0.35-8.8 
0.35-48 

38-630 
3510-26200 

0.7-2.4 
4.2-49 

120-850 
110-350 

0.7-1.8 
0.7-0.7 

71-580 
84-320 



Ipiant material: AG=above-ground, BG=below-ground 

2values below the instrument detection limit were estimated by multiplying 
the instrument detection limit value by 0.7. 

present at both wetlands, while Sphagnum tenellum was found only 
at the Swamp Gulch site. 

Table 9-8 presents tissue concentration data for various 
bryophyte species in metal-enriched and non-enriched environ- 
ments. The range of values for a particular element varies 
widely due to the different species analyzed and because of the 
different environments in which they were found. For some of 
the elements (Fe, Pb, Mn, Ni and Zn) , the range of values for 
the non-enriched areas overlaps the range for the enriched Table 
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areas. Despite this apparent anomaly, these data indicate the 
general range of element concentrations for mosses. 

With the exception of Mn, the element levels in 
Isoptervaium pulchellum (Table 9-9) were within, or close to, 
the range for other mosses from enriched environments (Table 9- 
8) . Data from previous studies indicate a high level for Mn of 
412 ug/g in Sphagnum spp. from an artificial wetland that 
received AMD (Wieder et al. 1985). In contrast, Isoptervaium 
pulchellum accumulated 1423 ug/g of Mn at the Swamp Gulch 
wetland. This difference is likely due to the inherent 
abilities of Isoptervaium pulchellum to accumulate more Mn than 
the Sphaanum species studied by Wieder et al. (1985). However, 
this difference may also be due, at least in part, to the length 
of time that the two mosses were exposed to AMD. The Sphaanum 
moss was exposed over a 10 month period, whereas Isoptervaium 
pulchellum had been growing at the Swamp Gulch site for years. 
In any case, Isopteryaium pulchellum apparently can absorb and 
store relatively large amounts of Mn. 

With the exception of As and Ca, all elements were elevated 
in Isopterygium pulchellum at the Swamp Gulch wetland compared 
to background levels (Table 9-9) . This result was expected for 
As since it was not elevated in the influent water (Table 9-4) . 
The pattern for calcium was similar to that for the other plant 
species since the concentration was lower on the Swamp Gulch 
compared to the background wetland. 

Tissue concentrations of all the elements except Ca were 
lower in Sphaanum tenellum than in Isoptervaium pulchellum at 
the Swamp Gulch wetland. The concentration of elements in 
Sphagnum tenellum were within the range or were lower than the 
range for other mosses from enriched environments (Table 9-8) . 

9.3.4 Summary of Element Enrichment in Plants 

Table 9-10 lists those plant species analyzed at both 
wetlands and the elements that were enriched above background. 
Iron, Cu, Pb and Zn appeared to be enriched in all species. The 
moss, Isopteryaium pulchellum . accumulated all elements except 
Ca. Significant phytotoxicity of any of these elements is 
probably not occurring at the study site. Plant enrichment of 
Cd is of particular concern because the threshold level for safe 
intake of the element by animals is relatively low. This 
potential problem is discussed in subsequent sections. 
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Table 9-10, 



Plant species demonstrating element enrichment 
(i.e. greater than background levels) at the 
Swamp Gulch wetland site. 



Plant Species 



Enriched Elements 



Carex rostrata 
Salix booth ii 
Betula glandulosa 
Isopteryqium pulchellum 



Al, Cd, Cu, Fe, Pb, Zn 

Al, Cd, Cu, Fe, Pb, Mn, Zn 

Cd, Cu, Fe, Mn, Pb, Zn 

Al, Cd, Cu, Fe, Pb, Mn, Ni, Zn 



9.3.5 Vegetation Element Levels and Wildlife Hazards 

Concern has been expressed that the construction of 
artificial wetlands for AMD abatement may result in high levels 
of metals in the vegetation and that this would contaminate or 
poison the wildlife resource. The following sections address 
this concern by comparing animal hazard criteria levels with 
plant tissue data demonstrating elevated levels of the target 
elements. 

9.3.5.1 Dietary Hazard Levels for Animals 

There is a general lack of data on hazardous metal intake 
levels for wildlife species. However, there is considerable 
data for domestic animals. Table 9-11 lists the maximum 
tolerable levels of minerals in diets of cattle, horses and 
domestic rabbits. When these levels are fed for a limited 
period, they will not impair animal performance and should not 
produce unsafe levels in human food derived from these animals. 
These threshold values were generated from a variety of studies 
in which graded levels of the elements were ingested and the 
specific affects examined. Highly soluble and purified forms of 
these elements were generally used in the studies. This 
practice tended to bias the threshold levels low. In other 
words, a particular concentration of a metal in the diet may 
exceed the threshold level but not be dangerous. In addition, 
adverse physiological effects may not occur at a particular 
threshold level if the element has been complexed into plant 
matter during growth and is not digested by the animal. Animal 
tolerance to chemical elements will vary with age and 
physiological condition. For the purpose of discussion, these 
data (Table 9-11) will be used as a general guide for potential 
contamination of wildlife species utilizing the Swamp Gulch 
wetland. 
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Table 9-11. Maximum tolerable levels^ of dietary minerals for 

domestic animals (National Academy of Sciences 1980) 







Species 




Element 


Cattle 


Horse 
(200) 


Rabbit 


Al 


1000 


(200) 


As 


50 


(150) 


50 


Cd 


0.5 


(0.5) 


(0.5) 


Ca% 


2 


2 


2 


Cu 


100 


800 


200 


Fe 


1000 


(500) 


(500) 


Pb 


30 


30 


(30) 


Mn 


1000 


(400) 


(400) 


Ni 


50 


(50) 


(50) 


Zn 


500 


(500) 


(500) 


^ All values are in 


ug/g unless 


Otherwise spe 


cified. Values in 



parentheses were derived by extrapolation. 

9.3.5.2 Wildlife Observed and at Risk 

Wildlife observations were recorded during each field trip. 
Voles (Microtus spp. and possibly Clethrionomvs aapperi ^ were 
the animals most often observed. They were generally seen among 
Carex rostrata where the soil was saturated or nearly so. One 
cottontail rabbit f Svlvilacrus nuttalli) was seen in the center 
of the site. The only ungulates actually observed were mule 
deer (Odocoileus hemionus) . A total of six adult deer were seen 
during the study period while a fawn and doe were observed 
several times. Tracks of moose (Alces) and possibly elk f Cervus 
canadensis) were seen twice. Their use of the area was minimal. 
Evidence of beaver ( Castor canadensis ^ use of the site was low. 
Beaver tracks were observed only once and freshly chewed willow 
stems were seen on several trips. Beaver were probably moving 
through the area but they were definitely not residents of the 
site. Brook trout ( Salvelinus fontinalis ^ were often in the 
larger channels in the southeast portion of the wetland study 
area. An approximately eight inch brook trout was seen in a 
channel near station AAA-3. Only very small trout were observed 
m the Blackfoot River. Numerous frogs were observed throughout 
the Swamp Gulch wetland. Coyote ( Canis latrans) tracks were 
seen on only one occasion. Very few passerine birds were 
observed. A redtail hawk ( Buteo iamaicensis ^ and presumably a 
great horned owl (Bubo virainianus ) were observed once. 

Animal species most likely to be poisoned by heavy metals 
are those whose home ranges encompass a metalliferous environ- 
ment. In the case of the Swamp Gulch wetland, these would be 
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primarily herbivores such as voles, muskrats ( Ondatra zibethica ) 
and beaver, and carnivores such as weasel (Mustela frenata) and 
mink ( Mustela vison ) . Other possible residents of the study 
site that could be at risk are shrews ( Sorex spp.) and the skunk 
(Me phitis ) . Mammals and birds that feed over a large area are 
not likely to be poisoned by occasionally consuming vegetation 
or contaminated animals from the Swamp Gulch wetland. These 
include ungulates such as elk, moose and mule deer, and 
predators/ scavengers such as coyote, fox and various raptors. 
This discussion will focus on those species with small home 
ranges, meaning those more likely to consume large quantities of 
vegetation or other animals from the Swamp Gulch wetland. 

9.3.5.3 Hazards to Wildlife at the Swamp Gulch Wetland 

Cadmium appears to be the only element present in the 
macrophytes (Tables 9-3, 9-6 and 9-7) at levels that could be 
toxic to animals. The bryophytes contain metal levels that 
exceed maximum dietary intakes, however, these plants are not 
consumed in quantities that could cause metal poisoning. Of the 
macrophytes, Salix boothii had elevated Cd levels on the Swamp 
Gulch wetland that greatly exceeded the maximum dietary intake 
level. The Cd level in the above-ground tissues was 5.6 ug/g 
which was significantly greater than the dietary tolerable level 
of 0.5 ug/g (Table 9-11). Beavers consume large quantities of 
Salix boothii . so could be at risk for Cd poisoning. If beaver 
have been poisoned by consuming Salix boothii . it would help 
explain the absence of resident beavers at the site. Beaver are 
currently active in areas upstream and downstream of the Swamp 
Gulch wetland and were once very active at the site. 

In summary, Salix boothii on the Swamp Gulch site had a 
level of Cd that exceeded the maximum tolerable dietary intake 
for domestic animals. Beaver appear to be at particular risk 
from Cd poisoning because of their heavy utilization this shrub. 
Other animals, such as voles, may also be at risk. If these 
small rodents are accumulating high levels of Cd then the 
animals that prey on them (weasel, mink, owl) may also be 
experiencing toxic effects. It is less likely that infrequent 
visitors to the Swamp Gulch wetland such as elk, deer, moose, 
coyotes, foxes and birds would be poisoned. 

At this point, it is unknown whether true metallophytes 
exist at the Swamp Gulch study site. If these types of plant 
species are present, certain animals may be at special risk if 
they consume these plants. Future work at this site should 
include an examination of metal levels in tissues of the common 
herbivore (s) . These data would help determine the extent to 
which certain metals are present in the food chain. 
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10.0 MICROORGANISMS AND THE IMMOBILIZATION OF HEAVY METALS 
IN WETLANDS 

10.1 BACKGROUND 

This portion of the study characterized microbial com- 
munities and conditions associated with AMD in constructed and 
natural wetlands using a new method for measuring the popu- 
lations of microorganisms based on their unique biochemical 
components. MSI Detoxification, Inc. (MDI) uses this technology 
to characterize hazardous waste sites and proposed its 
implementation in this study. This new method allows the 
analysis of mixed microbial populations found in natural 
environments and does not require cultivation of organisms in 
the laboratory. Thus, it is a powerful tool which avoids 
isolation of microorganisms by traditional cultivation 
techniques, and therefore, may be more representative of the 
microorganisms present in situ. As a part of the study of the 
efficacy of natural wetlands in the remediation of acid mine 
drainage, it was determined that the role of microorganisms, 
specifically sulfate reducing bacteria, be investigated in the 
unique circumstances presented by the Swamp Gulch natural 
wetland in Lincoln, Montana. A reference for this wetland was 
provided by the Sand Coulee wetland, a constructed wetland near 
Tracy , Montana . 

The hypothesis for the perceived success of natural 
wetlands in remediating AMD considers that microorganisms have 
the ability to carry out the reduction of in-flowing sulfate. 
This, in turn, raises the pH of the system and produces sulfide. 
The sulfide may be involved in reductions of ferric iron (Fe-^"*") 
to ferrous iron (Fe2+) and subsequent immobilization, possibly 
as inorganic or organic sulfides. If the process (or the 
process initiator) is microbial rather than chemical, then 
sterilized soil samples should not cause the pH of AMD to rise 
after incubation under anaerobic conditions. 

To test this hypothesis in the Swamp Gulch environment, 
experiments were performed under which pH, sulfate, iron, and 
microbial biomass would be measured over time and related to AMD 
remediation. The experimental approach was divided into three 
Phases. Under Phase I, the microbial component of the Swamp 
Gulch AMD remediation activity was investigated through a 
comparison of remediation capacities in sterilized versus 
untreated core materials (soils, organic detritus, inorganic 
materials) . Phase II sought to verify the presumed microbial 
activities observed under Phase I by following pH, sulfate, and 
iron in wetland core materials in a controlled greenhouse 
environment. Under Phase II, materials from the Swamp Gulch 
natural wetland were compared with sediment core materials taken 
from the Sand Coulee constructed wetland. Phase III focused on 
the microbial component of the wetlands and compared biomass and 
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(inferred) microbial diversity in field samples extracted from a 
natural, constructed and control wetland (Hardscrabble Creek 
near Lincoln, Montana) . 

The characterization of the major types of microorganisms 
in artificial and natural wetlands and their correlation with 
effectiveness of AMD remediation are intended to allow us to 
develop rapid and effective methods to diagnose and renovate 
malfunctioning wetlands. 

10.2 PHASE I: SWAMP GULCH WETLAND BUFFERING CAPACITY AND THE 
EFFECT OF SOIL STERILIZATION ON WETLAND PH 

10.2.1 Introduction 

A preliminary experiment to compare the abilities of Swamp 
Gulch wetland soil samples to remediate AMD under various 
conditions was performed. The purpose of this experiment was 
two-fold: 1) to gain insight regarding the buffering capacity 
of Swamp Gulch soil on exposure to AMD, and 2) to determine if 
there was evidence that microorganisms beneficially raise the pH 
in Swamp Gulch wetland soil. 

10.2.2 Technical Approach 

Soil samples from the Swamp Gulch wetland were incubated 
with AMD under aerobic and anaerobic conditions. Initial pH 
values of the soil-AMD mixture and those of the same suspensions 
after incubation for nine days were compared. 

10.2.3 Materials and Methods 

10.2.3.1 Soil Sampling 

Soil samples (0.5 and 2.0 g) were obtained from cores taken 
from the Swamp Gulch wetland (Figure 10-1) . Each sample was 
placed in a plastic "Whirl-pak" bag and transported to the 
laboratory where it was transferred to a Waring Blender and 
homogenized for a total of 60 seconds. Soil samples (0.5 g and 
2.0 g) were weighed into tared 30 ml beakers (titration experi- 
ment) or Hungate roll-tubes (microbiological experiment) , 
respectively. 

10.2.3.2 Titration Experiment 

To determine soil buffering capacity, duplicate 0.5 g soil 
samples suspended in 10 ml of distilled water were titrated with 
AMD. Readings of pH were made for six minutes or until a stable 
value was obtained (< ten minutes) . Similarly, AMD was titrated 
into distilled water for comparative purposes. AMD titrated 
into distilled water was not replicated. All samples were 
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Figure 10-1. Location of wetland sites. 



continuously mixed using a magnetic stirring device during the 
titration. 

10.2.3.3 Microbiological Experiment 

In the microbiological study all sampling was carried out 
in triplicate and three treatments were considered: a sterilized 
control, an aerobic and an anaerobic treatment. All samples 
were incubated under ambient laboratory temperature conditions 
for a period of nine days. For the sterilized control, soil 
samples (2.0 g) suspended in 14 ml water or 14 ml AMD were 
sterilized on three consecutive days by autoclaving for 15 
minutes at 15 psi. Anaerobic treatment samples were gassed with 
Oo-free N2 and stoppered with Hungate stoppers. Aerobic 
treatment samples were incubated in 125 ml Erlenmeyer flasKs 
with loose fitting caps and incubated with vigorous shaking. 
Readings of pH were made initially and after the incubation 
period. 

After incubation and pH measurement, triplicate samples of 
each treatment were pooled to provide sufficient sample for iron 
analysis. The pooled samples were centrifuged at 10,000 x G for 
ten minutes and the supernatant fraction decanted. This 
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fraction was acidified with concentrated HCl and total dissolved 
iron determined by the phenanthroline color imetric method (APHA 
1985, Iron Determination, pp. 215-220). 

10.2.4 Results and Discussion 

10.2.4. 1 Titration Experiment 

Figure 10-2 depicts the titration of wetland soil 
homogenate and distilled water. After a volume of AMD was added 
to the stirred soil suspension, the pH fell to a minimum of 4.5 
during the six minute measurement period. This indicates that 
the soil has a capacity to adsorb acidity independent of 
microbial involvement, since the experiment was conducted on 
time scales too short for significant microbial activity. 

J-0-2-4.2 — Microbio logical Experiment 

The incubation under anaerobic conditions of AMD soil from 
the Swamp Gulch wetland raised the pH of the AMD by 0.72 units 
(Table 10-1) . When the soil and AMD were sterilized, no such 



Sample Sources: Swamp Gulch Soil & AMD Water 



PL, 




12345G70 9 

Volu me AMD Added, niL 

Figure 10-2. Titration of Swamp Gulch soil with AMD. 
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Table 10-1. pH and iron values^ for AMD incubated under various 
treatment conditions after a period of nine days. 



Pre Post 
Treatment Incubation pH Incubation pH 



Total 
Net Dissolved 
Difference Iron. mq/L 



Aerobic 
Anaerobic 

Sterile 
Control 



7. 08+0. 03a 
6. 28+0. 03b 



6. 35+0. 05b 



5. 17+0. 13c 
7. 00+0. 05a 



-1.92±0.10a 
+0.72+0. 06b 



162.9+2.5ab 
20. 7+0. O'^ 



6. 37+0. 10b +0.02+0.06c3 172. 8+5. la 



^ Means + standard deviation; n = 3 for pH data; n = 2 for total 

dissolved iron data. 
2 Treatment means (within columns) are significantly different 

(P=0.01) if followed by a different letter. 
^ The sterile control pH means from preincubation to postincu- 

bation were not significantly different from one another. 



increase in pH took place. Anaerobic conditions were necessary for 
this process. When AMD was incubated aerobically, the pH fell by 
1.92 units and the soil became bright orange, presumably due to the 
production of iron oxides. In all anaerobic incubations the soil 
maintained its dark brown/black color which would be expected if 
metal sulfides were being produced. Comparison of iron 
concentrations in the presence and absence of oxygen indicate that 
there was much less iron in the anaerobically treated sample. 
Presumably this anaerobic response is the consequence of microbial 
sulfate reduction to sulfide followed by the precipitation of iron 
sulfide, which is black in color. To verify that sulfides were 
present the soil was treated with 2N HCl and an H2S odor was 
detected. 

The Swamp Gulch wetland soil was able to raise the pH of the 
applied AMD from 2.5 to about 5.0 in the absence of microbiological 
contribution (titration experiment) . As more AMD was applied to 
the soil, this dropped to 4.5 (Figure 10-2). Because the pH scale 
is logarithmic, this represents a 50 to 100-fold increase in the 
ability of the soil to neutralize acidity (as hydrogen ions). 
Aerobic incubation resulted in a final pH (near 5) which was 
similar to that of the soil titration experiment. However, 
anaerobic incubation raised the pH to 7.0 (4.5 pH units) indicating 
that wetland soil combined with its indigenous anaerobic 
microorganisms was able to neutralize AMD acidity by a factor of 
greater than 10,000-fold. 
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These findings support the hypothesis that sul fate-reducing 
bacteria play a notable role in the remediation of AMD. Further 
evidence in support of the involvement of sul fate-reducing bacteria 
comes from the phospholipid fatty acid analyses of core samples 
taken from the same area of Swamp Gulch wetland (Phase III) . 

10.3 PHASE II: COMPARISON OF NATURAL AND CONSTRUCTED WETLAND 
SOILS IN THE REMOVAL OF ACIDITY, SULFATE AND IRON 

10.3.1 Introduction 

Biological oxidation of reduced sulfur due to the activities 
of sulfur oxidizing bacterial action on pyritic ores is the primary 
cause in the formation of AMD waste water problems. The resulting 
products in the waste are acidity, sulfate ion and iron. The 
sul fate-reducing bacteria remove sulfate by reducing it to sulfide, 
which can cause the co-precipitation of heavy metals as sulfides. 
The objective of these experiments was to compare acidity, sulfate 
and iron removal in a natural and an artificial wetland system. In 
addition these studies were carried using samples of sufficient 
size to minimize inhomogeneities which would be present in small 
samples. 

Iron, sulfate and pH measurements were included in the 
experimental protocol, as microbial sulfate reduction is the 
presumed basis for AMD remediation through immobilization of heavy 
metals. Iron was chosen as an indicator metal due to its ease of 
determination and its ubiquitous presence in the wetlands under 
study . 

10.3.2 Technical Approach 

The study was carried out in a greenhouse at the Montana State 
University Plant Growth Center. Three of the cores were obtained 
from the Swamp Gulch natural wetland near Lincoln, Montana (Figure 
10-1). These were designated SG-1, SG-2, and SG-3 for this study. 
The other three cores were obtained from the Sand Coulee 
constructed wetland large cell near Tracy, Montana (Figure 10-1) . 
These cores were designated SC-V, SC-VN, and SC-CH. 

The cores received AMD that had been collected at the Swamp 
Gulch site. AMD was eluted through all of the cores for the 
duration of the experimental study. The effluent (eluate) as well 
as the influent AMD (eluent) was analyzed for temperature and pH 
several times a week. Sulfate and iron concentrations were 
determined periodically over the course of the study. 

Direct microbiological parameter measurements were not made 
under this Phase of the study. Microbial diversity and biomass 
analyses comprised phase III of the investigation. 
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10.3.3 Materials and Methods 

10.3.3.1 Field sampling 

A site was selected at the Swamp Gulch wetland and was sampled 
to provide sediment cores for the greenhouse study. The selected 
site was between sites C3 and C2 (see Figure 4-1). It was assumed 
that the remediation capacity of the Swamp Gulch wetland was not 
saturated at this field location. Site vegetation was dominated by 
Carex with subsidiary Sphagnum moss. 

Three sediment samples (designated SG-1, SG-2 and SG-3) were 
collected on October 22, 1987. Each core vegetation mat was sawn 
along the sides and bottom with a keyhole saw, and removed from the 
wetland using a flat-bladed shovel. Each core was laid on its side 
on the ground surface and trimmed to fit tightly into the sediment 
container. The containers, five-gallon plastic buckets modified as 
shown in Figure 10-3, had been acid-washed and disinfected with 
chlorine bleach solution before they were taken to the field. 
After pressing the sediment firmly to remove air spaces along the 
container walls, native water from the sampling trench was poured 
over each sample so that standing water was visible just below the 
core top. A five-gallon container of the in situ water was 
collected so that it would be possible to replenish sample water in 
case of losses during transport. 

Sediment cores were transported to Bozeman the following 
afternoon. The next morning, samples were sprayed with MAVRIK 
Aquaflow, a broad-spectrum insecticide (active ingredient = 
(aRS, 2R) -f luvalinate [ (RS) -a-cyano-3-phenoxybenzyl- (R) -2-{2- 
chloro-4-trifluoromethyl) anilino}-3-methylbutanoate] ) . 
Fumigation, required of users of the MSU greenhouse facility, was 
carried out to prevent the spread of disease among experimental 
plants. It was believed that this treatment did not affect the 
microbial or chemical integrity of the samples due to its focused 
surface application. 

Samples from the Sand Coulee site near Tracy, Montana were 
collected on October 30, 1987, and placed into five-gallon 
containers in an identical manner as were those collected at the 
Swamp Gulch wetland (Figure 10-3). The three samples were 
collected about one meter from the middle baffle of the first cell 
of the wetland, in an area with sporadic Tvpha plants and pools of 
standing water on an irregular sediment surface. Such vegetation 
cover was typical for that cell (Figure 10-4). The absence of a 
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vegetative mat and the high water level in the artificial wetland 
precluded removal of a coherent core. Instead, the saturated 
sediment was shoveled into the container, preserving the depth 
sequence of the peat. This procedure was deemed acceptable for the 
Sand Coulee wetland samples because the sediment color and texture 
were homogeneous throughout the peat profile,- and because the 
sediment was fluid enough to fill the container without forming 
pockets of air. 
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Cores were extricated from the Sand Coulee constructed wetland 
large cell area (designated "SC" cores) for a variety of treat- 
ments. The core designated SC-V contained cattails which were 
cropped for this study. The core designated SC-VN represented a 
sample void of vegetation during sampling. The core designated 
SC-CH contained cattails and was treated with chloroform for a 
period of about two months prior to beginning this study at which 
time the chloroform addition was halted. The intent in adding 
chloroform was to chemically sterilize the core. The Sand Coulee 
constructed wetland cores were treated as replicates and studied in 
the same manner as were the Swamp Gulch wetland cores. 

10.3.3.2 Soil core flow characteristics 

The effluent flow rate from each of the cores varied, 
presumably due to flow characteristics unique to each wetland 
sample. Flow rates were easily controlled for most cores. 
However, the effluent flow rate of the SC-V core decreased after a 
period of time so that less than 200 ml per day would freely flow 
from the core unless extracted using a vacuum device. This 
impediment of flow in the SC-V core occurred 26 days into the study 
and the average overall flow rate dropped to 125 ml/day. 
Consequently, its behavior during the first 26 days was similar to 
the five other cores and subsequently was extremely slow (57 
ml/day) . The flow rate for all other cores ranged from 253 to 291 
ml/day. 

10.3.3.3 Sampling and preparation of AMD 

AMD was collected from the Swamp Gulch natural wetland site by 
siphoning AMD from a receiving pond, located adjacent to the 
wetland and on the north side of Montana Highway 200 (see Figure 3- 
1) . The collected AMD was stored in plastic 55-gallon drums that 
had previously been acid-washed with dilute HCl and rinsed several 
times with distilled water. The 55-gallon drums were stored closed 
in the greenhouse. 

As AMD water was needed for the experiment, it was siphoned 
from the 55-gallon drums and filtered through glass wool. The 
filtered AMD was stored in 5-gallon carboys under cover to avoid 
photooxidation resulting from exposure to sunlight in the 
greenhouse. This filtered AMD was used as the eluent for both the 
Swamp Gulch and Sand Coulee wetland cores . Temperature and pH 
determinations were made on the AMD eluent (influent) and eluate 
(effluent) collected from each core, every day effluent was 
sampled. 

10.3.3.4 Sulfate and iron determinations 

Sulfate and iron concentrations were determined several times 
during the course of the study period. Dissolved sulfate 
determinations were made on days 61, 82 and 88 of the study. 
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Dissolved ferrous and total iron measurements were made on days 42, 
61, 82 and 88. The schedule for these analyses was based on an 
estimate of the percent of one pore volume of water which had been 
eluted from each of the soil cores. The percent water contained in 
the soil was determined to be approximately 85% based on percent 
moisture determinations. Thus, the water contained in each of the 
core was estimated to be 16.9 L. The core effluent was sampled 
when approximately 13.5 and 20.3 L had passed through each core and 
again at the end of the study. Due to the innate low flow rate of 
the SC-V core it was not possible to obtain iron and sulfate data 
beyond 12.2 L of flow- through volume. 

Sulfate concentrations were determined by ion chromatography 
(subcontracted to Camas Laboratories, Inc., Missoula, Montana) 
using standard methods of analysis (APHA 1985, Sulfate Determina- 
tion, pp. 483-488) . Two-hundred mL aliquots of eluate for sulfate 
analysis were collected from each core and filtered through 
cellulose acetate membrane filters (nominal exclusion 0.45 um) . 
All of these samples were refrigerated (4°C) overnight, packed on 
ice in a cooler and transported to Missoula, Montana, for chemical 
analysis. All sulfate analyses were performed within one to three 
days of receipt of the samples. 

Iron analyses were performed using the phenanthroline 
colorimetric method (APHA 1985, Iron Determination, pp. 215-220). 
Eluate samples were collected for iron analysis at the same time 
and in a similar manner as were the eluate samples for sulfate 
analysis. For dissolved ferrous iron determinations the core 
eluates and AMD eluent were collected directly into HCl to acidify 
them immediately and minimize oxidation of the ferrous iron. All 
dissolved iron samples analyzed were determined on the same day as 
they were collected; the dissolved ferrous iron samples being 
determined within nine minutes of collection, and the total 
dissolved iron samples deteirmined within eight hours of collection. 
It should be noted that the collection of the samples from the SC-V 
core, which did not flow freely, was performed under a vacuum 
pressure . 

10.3.4 Results and Discussion 

10.3.4.1 Acidity determinations 

The temperatures of the influent AMD and the effluents from 
the cores were generally in the 18-2 1°C range. The influent AMD pH 
ranged from 2.17 to 2.70, with a mean and standard deviation 
of 2.44 + 0.09 (N=94). This indicated that AMD applied to the 
various cores remained stable with respect to pH over the study 
period (Figure 10-5) . 

The Swamp Gulch wetland core eluates exhibited pH values with 
nearly identical trends indicating that they were replicates 
(Figure 10-6) . The effluent pH values for the Swamp Gulch wetland 
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cores ranged from 5.67 to 7.58, with an overall mean and standard 
deviation of 6.43 + 0.37. 

The Sand Coulee wetland core eluates, while not true 
replicates, demonstrated some important trends. Eluate pH values 
from cores SC-V and SC-VN were similar (Figure 10-7). These pH 
values dropped at the start of the study, flattening out by 30 days 
into the experiment and ranged from 2.39 to 4.25, with an overall 
mean and standard deviation of 3.24 + 0.45. The Sand Coulee 
wetland core which had been treated with chloroform, SC-CH, showed 
an ascent in its pH for the initial 30 days of the study and a 
leveling of this pattern subsequently. The SC-CH eluate pH ranged 
from 4.25 to 5.73 with a mean and standard deviation of 5.12 ± 
0.39. 
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Figure 10-5. 



Influent (eluent) AMD pH measured over study 
period. 
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Figure 10-6. 



Swamp Gulch soil core's effluent (eluate) pH 
measured over study period. 
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Figure 10-7. 



Sand Coulee soil core's effluent (eluate) pH 
measured over study period. 
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For the SC-V and SC-VN cores for most of the study period, 
Sand Coulee constructed wetland soil raised pH of applied AMD from 
2.5 to about 3.0. It is questionable whether there was any acid 
neutralizing microbial component in these cores. Consequently, the 
0.5 increase in pH may reflect the innate acidity neutralizing 
capability of the Sand Coulee wetland sediment. 

It can be surmised that the application of chloroform to the 
SC-CH core influenced the rise in pH. Although the eluate pH of 
the SC-CH treatment never achieved the near-neutral pH (6.43 + 
0.37) typified by the Swamp Gulch wetland cores, it demonstrated 
remediation of acidity not evident in either of the other Sand 
Coulee wetland cores. 

Two hypotheses can be advanced which would explain this 
observation. Both of these are based on the principal of microbial 
selection. In the first, a population of microorganisms arises 
because the chloroform is lethal to most of the indigenous 
organisms. A population is selected due to exclusive effects 
imposed by the chloroform on most of the organisms. As this group 
of microorganism proliferates, its metabolic activities result in 
increased pH and amelioration of the acidity. In the second 
hypothesis, selection occurs through nutritional enrichment of a 
microbial population able to use the chloroform as a source of 
carbon and energy. Because they are able to use the chloroform, 
they come into dominance. Presumably they had the ability to 
remediate the AMD pH prior to the addition of the chloroform, but 
were unable to because they were starved for a carbon source. 

Because there was no apparent lethal effect of the chloroform 
in the SC-CH core, and because the cattails in this core did not 
die, the second hypothesis is more probable. Also, the latter 
hypothesis is supported by remediation of AMD pH in the Swamp Gulch 
wetland cores where a high level of detritus from decaying sedges 
and other plant material was present. Conversely, the Sand Coulee 
wetland cores lacked conspicuous amounts of detritus which would 
provide carbon and energy and allow indigenous microbes to 
proliferate. These findings indicate that the constructed wetland 
is ineffective in removing acidity, presumably because it lacks 
sufficient carbon to support indigenous microorganisms. Therefore, 
it should not be automatically assumed that poor performance in a 
constructed wetland is due solely to wetland size. This scenario 
is supported by observations on the phospholipid analysis of the 
Sand Coulee wetland in situ site (Phase III) . 

10.3.4.2 Sulfate determinations 

The influent AMD applied to the cores exhibited a sulfate 
concentration of 412 ± 56.1 mg/L, ranging from 336 to 482 mg/L 
(Figure 10-8) . The Swamp Gulch wetland cores are remediating the 
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Influent (eluent) AMD iron and sulfate 
concentrations measured over study period 
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added sulfate (Figure 10-9) . The eluate concentrations for the 
Swamp Gulch wetland core vary from to 158 mg/L with a mean and 
standard deviation of 59.5 ± 56.2 mg/L. In general, the ^^ „, 
remediation was better as the flow rates were slower (Table 10-2). 

Sulfate concentrations in the effluent increased as the flow 
rates were increased. This could be associated with the rate at 
which indigenous sul fate-reducing bacteria are able to facilitate 
sulfate reduction to sulfide, where greater retention times provide 
for greater sulfate reduction. This phenomenon would be 
particularly important where sul fate-reducing bacteria were either 
not present, present at low concentrations, or present, but not 
active (Phase III) . For example, in a constructed wetland if flow 
rates are high and/or channeling is evident, one would expect^poor 
sulfate-reduction performance to occur. Therefore channelization 
and high flow rates in constructed wetlands should be avoided by 
implementation of improved engineering design. 

Instead of remediating the sulfate concentrations the Sand 
Coulee wetland cores are releasing sulfate into the effluent (Table 
10-3 Figure 10-9) . The effluent concentrations for the Sand 
Coulee wetland cores range from 1385 to 2660 mg/L. Presumably this 
is because the in situ addition of sulfate at the Sand Coulee 
wetland site is much higher then that in the Swamp Gulch AMD input 
water, and it is beginning to be removed as the lower 
sul fate-containing Swamp Gulch AMD is being applied to the cores 
(Figure 10-9) . 

The flow rates for these cores (except SC-V) increased during 
the later part of the study, as they did for the Swamp Gulch 
wetland cores, and there was a subsequent decline m sulfate 
concentrations in the effluent, suggesting that leaching of excess 
sulfate had occurred in these cores. These high levels of sulfate 
suggest that sulfate reduction was not taking place in the Sand 
Coulee wetland soils. This could be due to the absence of 
sufficient carbon to allow the sul fate-reducing bacteria to 
proliferate (Lovley 1987). It could also be due to the acidity of 
the cores as sul fate-reducing bacteria do not grow well at pH 
levels <six (Widdel and Pfennig 1984). 

In the SC-CH core, chloroform may have killed the sulfate- 
reducing bacterial population that was present. Alternatively, 
chloroform treatment may have enriched other microorganisms, 
effectively displacing sul fate-reducing bacteria. Since the other 
Sand Coulee wetland cores did not reflect lower sulfate levels, it 
may have been that sulfate-reduction was not occurring as a 
consequence of unknown factors. 
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Table 10-2. Swamp Gulch sulfate concentrations vs. flow rate. 



Core Flow rate^ Mean sulfate^ 

Number fmL/dav) (jm/Id 

SG-1 219.3 35.3 ± 3.0 

401.7 12.7 + 3.6 

472.7 33.9 ± 2.5 

SG-2 213.5 0.0 ± 0.0 

255.7 49.0 + 6.2 

596.4 130.5 + 9.2 

SG-3 198.7 6.2 ± 5.7 

424.3 152.5 ± 7.8 

441.4 115.5 + 2.1 



1 n = 1 

2 n = 2 



Table 10-3. Sand Coulee sulfate concentrations vs. flow rate, 



Core Flow rate^ Mean sulfate2 

Number fmL/dav) (mg/L) _ 

SC-V 144.2 1399.0 ± 19.8 

82.4 2425.0 ± 173.9 

91.4 2540.0 ± 169.7 

SC-VN 243.9 1450.5 ± 26.2 

330.4 2353.5 ± 72.8 

584.3 1585.0 ± 35.4 

SC-CH 172.0 1588.0 + 25.5 

468.3 2245.0 ± 59.4 

340.0 2050.0 ± 169.7 



1 



n = 1 



2 n = 2 
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These findings indicate that sulfate remediation at the 
Sand Coulee constructed wetland was not occurring. In order to 
optimize wetland performance, construction specifications should 
favor the proliferation and activities of sul fate-reducing 
bacteria. This means that: 

° Anaerobic conditions must prevail within the sediment. 
Therefore slow flow rates and microbial population 
diversity enhancement is important. Varying wetland 
depth would also be beneficial. 

° Hetabolizcible carbon and other nutrients should be 
enhanced. This would include providing living plants 
and fertilizer. This would also increase the con- 
centrations of aerobic microorganisms (increased species 
diversity) whose activities would lead to more favorable 
conditions for the anaerobic sul fate-reducing bacteria 
(by removing oxygen and providing alternate carbon 
siibstrates) . 

° For the sul fate-reducing bacteria to proliferate, the pH 
of the constructed wetland should be near six or 
greater. Therefore it would be prudent to select wet- 
land construction materials that have an initial pH 
near six or above. The use of lime in a constructed 
wetland can help attain this goal. 

10.3.4.3 Iron determinations 

Iron analyses show a consistent total dissolved iron 
concentration in the influent AMD of 10.2 to 16.9 mg/L, 
averaging 13.3 mg/L (Figure 10-8). This iron was present in the 
ferric form, which was not surprising, since the source of the 
acid mine drainage was a shallow pond. 

The presence of ferrous iron in the effluent of the soil 
cores was not surprising, since the soil cores provide a 
reducing environment. It was unexpected that no apparent iron 
remediation occurred in any of the cores since all six of the 
wetland cores eluted more iron than was added to them by the 
eluent AMD. Conceivably the organic material acts as a cation 
exchange column. In the case of the Swamp Gulch wetland soils 
there appear to be hydrogen ions displacing metal ions at the 
exchange sites. These three cores remained fairly stable in 
their effluent iron concentrations over time. The eluate iron 
concentrations ranged from 77 to 196 mg/L for the Swamp Gulch 
wetland cores. These concentrations did not decrease over time 
or over % pore volume eluted (Figure 10-10) , nor did flow rate 
affect them. 

In the Sand Coulee wetland cores the iron concentrations 
averaged 203 + 136.6 mg/L. The iron concentrations of the SC-VN 
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(no vegetation) core effluent constantly declined (Figure 10- 
11). On day 42, the concentration of iron in this core effluent 
was 476 mg/L. By day 88, it had decreased to 35.5 mg/L. The 
SC-V (vegetation) and SC-CH (chloroform treated) cores were 
similar except that their iron concentrations were > 100 mg/L at 
the end of the study. The SC-CH core had high iron concentr- 
ations initially (450 mg/L) and began to decrease at a later 
date than the other cores (day 62) . By the final iron analysis, 
the Sand Coulee wetland iron levels had dropped by half. 

Thus all of the Swamp Gulch and Sand Coulee soil core 
eluate iron concentrations were higher then those of the 
influent (eluent) AMD. However, cores from each of the sites 
exhibit different trends. The Swamp Gulch wetland cores 
released about the same concentration of iron over time, 
irrespective of: 1) flow rate or, 2) percent pore volumes 
eluted. The Sand Coulee wetland cores eluted extremely high 
concentrations of iron early in the study, with a tendency 
toward decreasing iron concentrations with time. This may be 
related to flow rate or it may be a washing out of iron present 
prior to the beginning of the study. 

We suggest different hypotheses for these phenomena. In 
the Swamp Gulch wetland cores no apparent iron remediation took 
place. The iron seems to be released from the soil at a fairly 
steady rate suggesting that cation exchange with acidity 
(hydrogen ions) is occurring. 

In the case of the Sand Coulee wetland cores, the iron 
concentrations decreased over time. This could be indicative of 
either, 1) leaching of iron by the Swamp Gulch AMD, 2) or these 
cores may be remediating iron. Of these suppositions, it seems 
most likely that iron is being leached. When the total effluent 
volumes for the SC-V and SC-VN cores are compared with one 
another (Table 10-4) , the iron concentration in the effluent 
decreases as the flow through volume increases. In the SC-CH 
core, the chloroform may have altered the iron remediation 
properties . 

These findings indicate that iron removal at the natural 
wetland and at the constructed wetland sites was not occurring. 
In order to promote iron removal in constructed wetlands, 
materials should be selected that: 

° Provide as much cation exchange capacity as possible to 
remove acidity and metals. 

° Provide favoredsle conditions for the proliferation of 

sul fate-reducing bacteria. If these bacteria are able to 
reduce sulfate to sulfide, metals will be immobilized 
through to formation of metal sulfide precipitates. 
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Table 10-4. Sand Coulee iron concentrations vs. total flow 
volume eluted at end of study. 

Core Volume AMD^ Average Iron^ 

Number Eluted. mL ma/L 

SC-V 12,453 104.2 ± 1.1 

SC-VN 27,180 35.5 ± 0.0 

SC-CH 23,768 199.2 +3.5 



1 n = 1 

2 n = 2 



10.4 Phase III: MICROBIAL BIOMASS, COMMUNITY STRUCTURE AND 

PHYSIOLOGICAL STATUS ASSESSMENT OF NATURAL AND CONSTRUCTED 
WETLANDS 

10.4.1 Introduction 

Quantitative extraction of cellular components of sedimentary 
microbiota permits an assessment of the viable microbial community 
biomass, composition and, to a degree, physiological status, 
without the problems associated with direct enumeration or cultural 
methods (White 1986) . In particular, phospholipid ester-linked 
fatty acids (PLFA) have proven to be reproducible, capable 
indicators of a wide suite of microbiota encountered in the 
environment (e.g., Guckert et al. 1985). The phosphate of the 
microbial phospholipids has been shown to have a rapid turnover in 
both living and killed cells (White et al. 1979), indicating that 
an accurate analysis selective for this lipid class would permit a 
quantitative estimation of the viable microbiota. The objective of 
these studies was to use PLFA analysis to assist in the assessment 
of the role of microorganisms in the amelioration of AMD in natural 
and constructed wetlands by demonstrating the range of biomasses 
and the differences in the structure and physiological status of 
the viable microbial communities within these wetlands. 

10.4.2 Technical Approach 

A total of twenty wetland sediment samples were collected in 
the field and subjected to PLFA analysis in conjunction with the 
wetlands study. 
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10.4.3 Materials and Methods 

10.4.3.1 Field sampling 

Sediment samples were collected from several wetland sites: 
Hardscrabble Creek, Swamp Gulch and Sand Coulee (Figure 10-1) . 
Hardscrabble Creek is a natural wetland site. The Swamp Gulch 
wetland site receives inputs of AMD (pH =2.5) from the 
Carbonate mine. Hardscrabble Creek was chosen as a "control 
site" that did not receive AMD, but was in close proximity to 
Swamp Gulch. Sediment samples were collected from the 20-21 cm 
depth horizon for both sites in individual replicate cores. The 
Sand Coulee wetlands are constructed wetland cells which began 
operation within the past two years. In both Sand Coulee 
wetland cells studied, the surface sediment (approximately 1-2 
cm) and 20-21 cm depth horizon were sampled in independent cores 
collected along a transect across the wetlands. 

The principal differences in the wetland sediments observed 
included the dense vegetation mat and complex root systems of 
the natural wetlands which held the wetland sediment together 
into a cohesive unit, whereas the man-made wetland sediment 
lacked this vegetation and cohesiveness. All wetland sediment 
samples were frozen and then lyophilized (freeze-dried) . After 
drying, the samples were inventoried, randomly relabeled, and 
total dry weights were determined. 

10.4.3.2 Statistical analyses 

All statistical analyses utilized SPSSX statistics programs 
available on the Montana State University Honeywell 66/DPS 
mainframe computer with CP-6 operating system. Data transforma- 
tions were used to meet the homogeneity of error variance 
assumptions for the analysis of variance model used (Winer 
1971). In general, a log^g transformation was performed, except 
when the variable was a percentage, in which case an arcsine 
square root transformation was made (Winer 1971) . Following the 
analysis of variance (ANOVA) , Tukey's Honestly Significant 
Difference test (HSD, SPSSX programs) was calculated for a 
multiple comparison of means. Significant difference maps were 
produced keeping the within-experiment, family-wise error rate 
set at a = . 05 . 

10.4.3.3 Lipid nomenclature 

Fatty acids are named as the total number of carbon atoms, 
followed by the number of double bonds closest to the aliphatic 
(w) end with the two numbers being delimited with a colon (i.e., 
carbon atom number: double bond number) . Furthermore, the 
geometry ('c' for cis and 't' for trans ) of the double bond is 
also denoted with the double bond number (e.g., 18:lw7c). The 
prefixes 'i', 'a', and 'br' refer to iso, anteiso and methyl- 
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branching of the unconfirmed position, respectively. 
Cyclopropyl fatty acids are designated as 'cy' with the ring 
position in parenthesis relative to the aliphatic end (e.g., cy 
19:0). 

10.4.4 Results 

10.4.4.1 Total Biomass 

The total PLFA quantified can be used as an estimate of the 
biomass of the viable microbiota in the wetland sediments. The 
total nmoles of PLFA per gram dry weight of wetland sediment is 
shown for each site and depth horizon in Figure 10-12. His- 
tograms which are labeled with symbols for the same group (e.g., 
Bl and Bl+2, Figure 10-12) are included in the same homogeneous 
subset as defined by the Tukey HSD test. Histograms with no 
common symbols (e.g., Nl and N2, Figure 10-13) are significantly 
different by the Tukey test. The letters for these symbols 
indicate the parameter measured (B=biomass, N=nmoles/gdw, M=mole 
percent, T= trans/cis . C=cy c 1 opr opy 1 / c i s ) . The numbers differen- 
tiate the groups, with 1 being the group with the highest mean 
value of the parameter measured. 

At the 20 cm depth horizon, therefore, the Hardscrabble 
Creek and Swamp Gulch sites have a significantly greater micro- 
bial biomass than the Sand Coulee sites (Figure 10-12) . The 20 
cm horizons for the Lincoln sites are as dense in microbial 
biomass as the surface sediments of the Sand Coulee sites. 

10.4.4.2 Community structure 

In addition to estimating total microbial biomass, the PLFA 
profiles can be used to compare the microbial community 
structure for the sampled wetlands. Table 10-5 includes the 
PLFA profiles for all sites and depth horizons expressed as the 
mole percent of the entire PLFA profile. The fatty acids are 
listed in their elution order off the chromatography column. 
The data are expressed as average ± one sample standard 
deviation. Statistically significant differences within 
individual PLFA are marked with an '*• (ANOVA, a = 0.05). 

While the presentation of a complete data set such as shown 
in Table 10-5 is important, microbial community structure 
interpretations are difficult using individual PLFA. Table 10-6 
describes some PLFA groups assembled to aid interpretation. The 
eukaryotes are characterized by polyunsaturated fatty acids 
(Shaw 1966). The bacterial group is a combination of the 
terminally-branched fatty acids common to gram-positive-like 
organisms (Kaneda 1977) ; the PLFA indicative of sulfate-reducing 
bacteria (Bowling et al. 1986); the cyclopropyl fatty acids and 
their immediate precursor, the major product of the anaerobic 
desaturase fatty acid biosynthetic pathway, 18:lw7c (Fulco 
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Figure 10-12. Total microbial biomass as measured by analysis of 
phospholipid ester-linked fatty acids. 
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Figure 10-13. Sulfa te-reducing bacterial phospholipid fatty acid 
biomarker expressed as density per gram of soil and 
proportion of total microbial biomass. 
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1983) . The anaerobic desaturase group includes the w7c products 
and the corresponding cyclopropyl and mid-chain branched isomers 
resulting from transmethylation, as well as the trans -isomers 
which are hypothesized to be formed by an isomerase of the w7c 
products of anaerobic desaturase (Guckert et al. 1986). The 
sul fate-reducing bacteria can also be further defined based on 
the genera-specific PLFA biomarkers (Dowling et al. 1986; Edlund 
et al. 1985) . 

When the PLFA groups defined in Table 10-6 are used to test 
for changes in microbial community structure, the significant 
difference map shown in Table 10-7 is produced. The map is set 
up such that treatment means for each PLFA group increase from 
left to right, and those connected by a common line segment are 
not significantly different by the Tukey test (a = 0.05). 

The results presented in Table 10-7 will be discussed in 
more detail later. Briefly, the Sand Coulee wetland surface 
depth horizons (site codes = 3 and 5) were characterized as 
having the highest eukaryotic and lowest bacterial proportions. 
There were no significant differences in the proportions of 
gram-positive- like organisms or the total anaerobic desaturase 

Table 10-6. Description of phospholipid ester-linked fatty acids 
(PLFA) used to define microbial groups for data 
interpretation. 



I 
I 



MICROBIAL GROUP 



PLFA USED TO DEFINE GROUP 



Eukaryotes 



20:4w6 + 20:5w3 + 18 : 3w3 + 18 : 2w6 + 18 : 3w6 + 
16:lwl3t 



Bacteria 

Gr am-pos it ive 

Anaerobic 
Desaturase 

18 : lw7c 

Sulfate Reducing 
Bacteria 

Desulfobacter spp. 

Desul f ovibr io spp. 



il5:0 + al5:0 + il7:lw7c + I0mel6:0 + il7:0 
+ al7:0/17:lw9c + cyl7:0 + 18:lw7c + cyl9:0 

il5:0 + al5:0 + il7:0 + al7 :0/17 : lw9c 

16:lw7c + 16:lw7t + 10mel6:0 + cyl7:0 + 
18:lw7c + 18:lw7t + cyl9:0 

18 : lw7c 



il7:lw7c + 10mel6:0 + cyl7:0 
10mel6:0 + cyl7:0 
il7:lw7c 
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products. However, the principal unmodified product of anaero- 
bic desaturase, 18:lw7c, was in the greatest proportion in the 
Hardscrabble Creek site (site code = 1) . All sulfate-re- ducmg 
bacterial groups were in the greatest proportion of the total 
microbial biomass in the Sand Coulee 20 cm depth horizons (site 
codes 4, 6) with the lowest proportions being in the Swamp Gulch 
20 cm (site code 2) and Sand Coulee surface (site codes 3, 5). 

Since sul fate-reducing bacterial biomass and activity had 
been assumed to have a possible influence on AMD amelioration, a 
further illustration of sul fate-reducing bacterial distribution 
between the sites was done (Figure 10-13). The mole percent 
histograms graphically display the relationship shown in the 
significant difference map (Table 10-7) . The total sul fate- 
reducing bacterial biomass estimates in nmole per gram dry 
weight wetland sediment are also shown indicating that although 
the Swamp Gulch site has the lowest proportion of biomass as 
sul fate-reducing bacteria (group M2, Figure 10-13), the density 
of sul fate-reducing bacteria (per gram dry weight) for this site 
along with Hardscarbble Creek is included in the highest group 
(group Nl, Figure 10-13). In general, the Sand Coulee sites are 
lower in sul fate-reducing bacterial biomass density with the 
higher values occurring in the surface sediments. Greater total 
biomass and species diversity was evident at the Hardscrabble 
Creek and Swamp Gulch wetland sites. 

Even though sul fate-reducing bacteria are present in the 
Sand Coulee constructed wetland they are ineffective in removing 
sulfate. Where these organisms are present (at the surface) 
environmental conditions (e.g. , ©2 levels, nutrient 
availability) would preclude their activity. Therefore, 
specifications for wetland construction should address carbon 
and other nutrient additions as well as reduced oxygen levels 
necessary to enhance the densities and activities of these 
organisms . 

10.4.4.3 Physiological stress 

In addition to microbial biomass and community structure 
estimates the PLFA can be used to describe aspects of the com- 
munity physiological status. Specifically, the trans/cis ratio 
of monoenoic PLFA and the cyclopropyl/cis PLFA ratios have been 
suggested as possible "starvation or stress" lipid indices 
(Guckert et al. 1986). These ratios for each wetland site and 
depth horizon are shown in Figure 10-14 . The Sand Coulee sites 
can all be characterized by having higher cyclopropyl/sis ratios 
(group CI, Figure 10-14), whereas the Swamp Gulch site is char- 
acterized by an extremely high trans/cis ratio (group Tl, Figure 
10-14) . The control Hardscrabble Creek site is low for both of 
these "stress" lipid ratios (groups T2 and C2, Figure 10-14). 
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Table 10-7. 



Significant difference map generated from 
Tukey's HSD test^. 



Site 
Codes: 



1 
2 
3 
4 
5 
6 



Lincoln, Hardscrabble Creek (20 cm) 
Lincoln, Swamp Gulch (20 cm) 
Sand Coulee Large Bog (surface) 
Sand Coulee Large Bog (20 cm) 
Sand Coulee Small Bog (surface) 
Sand Coulee Small Bog (20 cm) 



MICROBIAL GROUP ' 

Eukaryotes 

Bacteria 

Gram Positive 

fbranched chains) 
(Anaerobic 

desaturase) 
18:lw7c 

Sulfate Reducing 
Bacteria 
Desulfohacter spp. 

Desulfovihrio spp. 



Low 

6 4 



High 
3 5 



1 
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5 
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6 


2 


3 


5 


1 


4 


6 


2 


3 


5 


1 


6 


4 



Within experiment familywise error rate set at a=0.05 
for 2*arcsine square root transformed mole % data of 
Table 10.5 grouped according to designations in Table 
10.6. 

2 

Treatment means for each group of fatty acids increase 
from left to right, and those connected by a common 
line segment are not significantly different for this 
test. 
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Figure 10-14. Microbial community trans/cis and cyclopropyl/cis 
phospholipid ester-linked fatty acids ratios. 
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10.4.5 Discussion 

The results of the PLFA data presented here need to be 
discussed in the context of results of the other projects 
conducted within the wetlands study. The Swarop Gulch wetland 
site has had in situ AMD inputs. The natural wetland sediments 
(Phase III) , as well as sediment cores brought into the labora- 
tory and continually flushed with AMD (Phase II) , are able to 
effectively buffer the AMD water. The cores analyzed in the 
laboratory have been shown to have a concomitant removal of 
sulfate from the AMD during the buffering process. When 
sediment cores from the Sand Coulee large wetland cell were 
exposed to AMD, no buffering occurred except when a core was 
treated with chloroform. No apparent sulfate utilization was 
measured from the Sand Coulee wetland cores. 

All of these wetland sediments are bacterial ly dominated 
communities (Table 10-5, note proportions of bacteria PLFA as 
described in Table 10-6) as opposed to a prevalence of 
eukaryotic (i.e., higher life forms such as green algae, fungi, 
etc.) organisms. The Sand Coulee wetland surface horizons show 
a higher eukaryotic input (Table 10-7) . The natural wetlands, 
which contain more vegetation and root material, support a much 
higher microbial biomass density than the constructed wetlands 
of the Sand Coulee cells (Figure 10-12). 

The apparent correlation of sulfate loss with AMD buffering 
suggests an important role for sul fate-reducing bacteria in the 
AMD amelioration process. Sul fate-reducing bacteria data 
displayed in Figure 10-13 indicate that, as was seen for total 
microbial biomass, the natural wetlands (Hardscrabble Creek, 
Swamp Gulch) support a higher sul fate-reducing bacterial biomass 
density than the artificial wetlands. Interestingly, the 
sul fate-reducing bacterial proportion of the total microbial 
biomass shows the opposite trend with the non-buffering Sand 
Coulee site having a higher proportion of the total microbial 
community as sul fate-reducing bacteria. The data in Figure 10- 
13 suggest several important aspects of AMD buffering by wetland 
systems: 1) If total sul fate-reducing bacterial biomass density 
is an important factor in the buffering process, we might 
predict that the Hardscrabble Creek control site would also be 
able to adequately buffer AMD; 2) The higher proportion of 
sul fate-reducing bacteria in the less diverse Sand Coulee 
microbial coxununities with the lower sul fate-reducing bacterial 
densities suggests the importance of other microbial populations 
to the sul fate-reducing bacteria. Sul fate-reducing bacteria are 
obligate anaerobes which are considered to be terminal carbon 
users when sufficient sulfate levels are present. Organisms at 
this nutritional level rely upon other aerobic, facultatively 
anaerobic and obligately anaerobic primary degraders to generate 
the carbon and energy substrates required for sul fate-reducing 
bacterial growth from the more complex molecular inputs into the 
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wetland (Dowling 1987) . Therefore, in the constructed wetland 
the availability of additional simple carbons sources (e.g., 
viable plants) would be beneficial to the sulfate-reducmg 
bacterial populations. 

This discussion suggests that if the primary degraders 
could be stimulated by an input of carbon and energy sources the 
sul fate-reducing bacterial populations might be better "fed," 
thereby increasing their density. If the density of sulfate— 
reducing bacterial populations is actually associated with AMD 
buffering, this artificial wetland system might then provide 
acidity and sulfate amelioration. One obvious difference in the 
wetlands that might account for the differences in microbial 
population proportions is the lack of vegetation in the Sand 
Coulee wetlands. Root systems provide nutritional inputs. 
Also, bioturbation of burrowing micro-eukaryotes associated with 
this vegetation can have an effect on the wetland sediment. _^ 
Note that after the Sand Coulee surface horizons the Swamp Gulch 
and Hardscrabble Creek sites also have higher proportions of 
eukaryotic biomass than the corresponding 20 cm depth horizons 
at the sand Coulee cells (Table 10-7) . 

It should not be assumed that wetland size alone governs 
performance. These findings indicate that the amount of 
microbial biomass is deficient in constructed wetlands when 
compared to natural wetlands. Enhancement of the microbial 
components would be expected to markedly affect the amelioration 
of sulfate and acidity due to AMD. 

The PLFA group made up of the products of the anaerobic 
desaturase fatty acid biosynthetic pathway shows no significant 
differences in their proportion to the overall microbial 
community for any of the sites (Table 10-7) . However, propor- 
tions of the immediate biosynthetic products (18:lw7c and 
16:lw7c) and their modification products (cyclopropyl and trans ) 
differ between sites. Table 10-7 indicates that 18:lw7c is 
significantly higher in proportion for the Hardscrabble Creek 
profiles, followed by Swamp Gulch and the 20 cm depth horizons 
of Sand Coulee wetlands. 

Figure 10-14 indicates that the cyclopropyl/ci^ ratios and 
trans/cis ratios also differ between sites, suggesting a 
difference in the microbial community physiological status. A 
full discussion of the lipid biochemistry and microbial physiol- 
ogy behind these proposed microbial "stress" indices is given in 
Guckert et al. (1986). The interpretation of the higher 
cyclopropyl/cis ratios at the Sand Coulee site must take into 
consideration the input of cyclopropyl PLFA from 
sul fate-reducing bacteria, where physiological response of 
cyclopropyl PLFA is unknown. It is evident that the Sand Coulee 
microbial communities are stressed in a way that might be 
analagous to anaerobic 
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microcosms which have been shown to have dramatic loss of cis - 
monoenoic PLFA with the corresponding production of cyclopropyl 
PLFA (Guckert et al. 1985). 

The other dramatic physiological shift shown in Figure 10- 
14 is the extremely high trans/cis ratio in the Swamp Gulch 
microbial community. This ratio generally ranges from 0.01 to 
0.08 for most bacterial cultures and 0.01 to 0.09 for marine, 
mangrove and estuarine sediment surfaces (Guckert et al. 1986). 
When Vibrio cholerae was starved for 7 days, the trans/cis ratio 
was 0.85 (Guckert et al. 1986). Significantly higher trans/cis 
ratios (0.15) were found in a highly consolidated sediment 
surrounding the burrow of a mud shrimp in estuarine sediment 
(Dobbs and Guckert 1988) . Four hours after an estuarine 
sedimentary microbial community was disturbed by sieving, the 
trans/cis ratio increased to 0.95 (Findlay et al. 1988). It is, 
therefore, unusual to find an in situ trans/cis ratio above 0.50 
as was observed in the Swamp Gulch sediments. This finding 
might be a response to low pH, metal toxicity, or another 
unknown stress. The differences between the trans/cis and 
cyclopropyl/cis stress of the Swamp Gulch and Sand Coulee 
wetlands are unknown. The control site, Hardscrabble Creek, is 
low for both "stress" lipid indices. It was previously 
suggested that this community might buffer AMD as the Swamp 
Gulch wetland has. It is unknown if the trans/cis ratio would 
follow the same pattern following AMD input. 

The microbial populations at both Swcimp Gulch and Sand 
Coulee wetlands are physiologically stressed. These physiolog- 
ical indices represent a means by which wetland performance can 
be measured. For example, if amendments applied to constructed 
wetlands resulted in improved AMD eumelioration and the stress 
indices dropped in magnitude, this would indicate that the 
overall "health" of the microbial community improved. This 
would translate to enhanced acidity removal, sulfate control and 
metal immobilization. 
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Appendix A-2. Procedures for calciilating accuracy and precision 



Procedures fog calculating accuracy. 

1. From the summarized information on the accuracy form 
the percent recovery of the blind field standard and the 
recovery of each lab spiked sample pair are calculated 
separately as: 

%Recovery of BFS = VA x 100; or %Recovery of Spike = SSR-SR x 100 

VK ' SA 

where: 

VA = analytical value of BFS 

VK = known (or certified) value of BFS 

SSR = spiked sample results 

SR = sample results 

SA = spike added 

Perfect accuracy would be 100 percent recovery. 

2. Calculate the standard deviation of all pairs. 

1/2 



SD = 



(Recovery^ - Recoveryp,vr, . ) 



n-1 



where recoveryj is the individual recoveries, recoverygyg , 
is the average recovery, and n is the number of values. 

To validate recovery data, the individual recoveries are 
compared with the average recovery value to identify individ- 
ual values that lie outside the range or reasonableness. 
Chauvenet's criterion is used to identify individual recovery 
values that lie outside this range. 

To use Chauvenet's criterion, the screening variable is 
computed for recovery values that are suspected of lying 
outside the range of reasonableness. 

Screening Variable = (Recovery^ - Recoverygvg . ) /SD 

The calculated screening variable is then compared to the 
maximum allowable value (Table B-1) for the appropriate 
number of recovery determinations. The suspect recovery 
value is set aside (set aside values are called "outliers") 
if the calculated screening variable equals or exceeds the 
maximum allowable value. 

If outliers are identified using Chauvenet's criterion, a new 
average recovery and a new standard deviation are recalcu- 
lated using the remaining "good" values, and Chauvenet's 
criterion is reapplied. This procedure is repeated until all 
surviving recovery values pass Chauvenet's criterion. 
(Usually one application and one recalculation are enough.) 
The final average recovery and final standard deviation are 
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Appendix A-2. Continued. 



calculated from the "surviving" recovery values. The final 
average recovery value is used to eliminate any bias from 
the laboratory data. 

4. The range of uncertainty (R) in the recovery is then calcu- 
lated. 



+ R = + tSD/(n)l/2 



where; 



R is the range of uncertainty expressed as a 
percent 

t is the value of the t distribution for the 

selected confidence level (90 percent) and (n-1) 
degrees of freedom 

n is the number of samples 

SD is the standard deviation. 

The range of uncertainty, is used in conjunction with the 
average recovery to determine if bias adjustments are 
required. 

Together, the final average recovery value for BFS and lab 
spike and the corresponding range of uncertainties constitute 
the QA statements of accuracy for a particular sampling 
program. 

The completeness of accuracy data is that percentage of the 
total number of samples that remained after outliers are 
identified and set aside with Chauvenet's criterion. 
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Appendix A-2. Continued. 



Table B-1. Chauvenet's Criterion, 



CHAUVENET'S CRITERION 
FOR REJECTING A SUSPECTED VALUE^ 



Number of 


Samples 


n 




3 




4 




5 




6 




7 




8 




9 




10 




11 




12 




13 




14 




15 




16 




17 




18 




19 




20 




21 




22 




23 




24 




25 




26 




27 




28 




29 




30 




40 





Maximum Allowable Values for 
(Recoveryi - Recoveryayg) /sD 



1.901 
1.983 
2.015 
2.111 
2.164 
2.195 
2.214 
2.228 
2.279 
2.318 
2.348 
2.373 
2.393 
2.409 
2,424 
2.435 
2.445 
2.454 
2.462 
2.469 
2.475 
2.480 
2.485 
2.502 
2.517 
2.530 
2.543 
2.555 
2.634 



^Based on "t" distribution rather than the traditional "normal" 
distribution. 



"Individual Recovery = Recovery^ 



Average Recovery = Recoveryg 



vg 
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Appendix A-2. Continued. 



Procedures for calculating precision. 

1. From the summarized information on the precision form the 
relative percent difference (RPD) of each replicate pair 
(blind field replicates and laboratory duplicates) are 
calculated separately as: 

RPD = |S-D| X 100 
(S+D)/2 

where: 

RPD = Relative Percent Difference 
S = First sample value (original) 
D = Second sample value (duplicate) 

Perfect precision would result in 0% RPD. 

2. Any RPD value exceeding the control limit of +20% is evalu- 
ated as an outlier using the Dixon's Q method. This i^ethod 
compares the difference between the suspected outlier and the 
value nearest to it in size with the difference between the 
highest and lowest values. The ratio of these differences 
(without regard to sign) is the Dixon's Q. 

Q = [ suspect RPD value - nearest RPD value| 
(largest RPD value - smallest RPD valuej 

The critical values of Q for P=0.05 are given in Table B-2. 
If the calculated value of Q exceeds the critical value the 
suspect RPD value is rejected, and not used in the following 
precision QA statement calculations. 



3. 



Calculate the relative standard deviation (RSD) of each 
duplicate pair (field duplicates and lab duplicates are 
treated separately) . 

RSDeach pair = SD/Mean 



Calculate the RSD for all the pairs 
duplicates are treated separately) . 



(field duplicates and lab 



RSDoverall = 



n 

z 

i=l n-1 



pairs 



1/2 
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Appendix A-2. Continued. 



5. Calculate precision as a percent. 

. . ^n-lpairs ^ RSDoverall 
Precision (%) = x 100 

(n-lpairs) ^/^ 
where: t is value from the 90% probability level 



Calculate the maximum uncertainty for any individual test 
result as: 



X + xtSDoverall/(n-lpairs)^/^ 



where: x is the reported test result 

t is the value of the t distribution at the 90% 
probability level. 



7. The completeness of precision data is that percentage of the 
total number of RPD values that remained after outliers are 
identified and set aside with Dixon's Q ratio. 
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Appendix A- 2. Continued, 



Table B-2 . Critical values of Q* 


(P=0.05) 




Sample size 




Critical value 


4 




0.831 


5 




0.717 


6 




0.621 


7 




0.570 


8 




0.524 


9 




0.492 


10 




0.464 






4 



E.P. King. 1958. J. Am. Statisti. Assoc. Vol. 48, 531. 
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APPENDIX B 

Influent Flow Rates, Precipitation and Water Chemistry 
Data from The Swamp Gulch Wetland 
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26.54 


38,06 


27.05 


36 ,84 16.11 


1 4,75 


1 6 , 1 1 


12.23* 



Total 1141.32 
Ma;c 171.62 
Min 36.84 



Note: * Record maximum and minimum for the pe^riod 
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Appendix B-2. Precipitation at the Swamp Gulch wetland study. 



Date 



cm 



m 



Date 



cm 



in 



6/09 


0.25 


0.10 


6/10 


0.25 


0.10 


6/16 


1.02 


0.40 


6/17 


0.30 


0.12 


6/18 


1.27 


0.50 


6/21 


0.13 


0.05 


TOTAL 


3.22 


1.27 


7/01 


0.13 


0.05 


7/05 


0.13 


0.05 


7/09 


0.13 


0.05 


7/10 


2.34 


0.92 


7/16 


0.71 


0.28 


7/17 


5.72 


2.25 


7/18 


0.89 


0.35 


7/20 


0.08 


0.03 


7/21 


0.10 


0.04 


7/22 


1.42 


0.56 


7/30 


0.18 


0.07 


TOTAL 


11.83 
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0.08 
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0.13 


0.05 
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0.25 


0.10 
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0.10 


0.04 
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0.05 


0.02 


8/25 
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0.42 


8/26 


0.03 


0.01 


8/27 


0.03 


0.01 


8/28 


0.05 


0.02 


TOTAL 


2.62 


1.03 


9/27 


0.08 


0.03 


9/30 
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0.04 


TOTAL 


0.18 


0.07 
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0.08 


0.03 


10/16 


0.38 


0.15 


TOTAL 


0.46 


0.18 
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0.13 


0.05 


11/03 


0.20 


0.08 


11/14 


0.25 


0.10 


11/19 


1.07 


0.42 


TOTAL 


1.65 


0.65 



12/02 


0.36 


0.14 


12/04 


0.20 


0.08 


12/08 


0.18 


0.07 


12/09 


0.25 


0.10 


12/11 


1.52 


0.60 


12/12/87-01/11/88 




4.24 


1.67 + 


1/13 


0.18 


0.07 


1/14 


0.38 


0.15 


1/15 


0.10 


0.04 


1/20 


0.31 


0.12 


1/22 


0.23 


0.09 


1/23 


0.13 


0.05 


1/30 


0.20 


0.08 


1/31 


0.15 


0.06 


TOTAL 


12/87 - 


1/88 




5.92 


2.33 


2/01 


0.05 


0.02 


2/02 


0.18 


0.07 


2/03 


0.48 


0.19 


2/04 


0.20 


0.08 


2/06 


0.05 


0.02 


2/07 


0.64 


0.25 


2/08 


1.35 


0.53 


2/09 


0.28 


0.11 


2/10 


0.08 


0.03 


2/11 


0.10 


0.04 


2/12 


0.30 


0.12 


2/13 


0.05 


0.02 


2/14 


0.41 


0.16 


2/15 


0.66 


0.26 


2/17 


0.05 


0.02 


2/21 


0.76 


0.30 


2/22 


0.13 


0.05 


TOTAL 


5.77 


2.27 
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APPENDIX D 

Precision and Accuracy of Vegetation Data, and Plant Species 
Observed at the Swamp Gulch Wetland 
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Appendix D-1. Accuracy and precision (at 907, con-fidence level) 
of vegetation data. 



El ement Accuracy l PrecisionZ 

Al 44.5+ 6.37. 

As 91.4+36.27. 
Cd NC 

Ca 85.8+ 4. 27. 

Cu 80.8+ 5.87. 

Fe 60.0+ 2.97. 

Pb 90.2+12.67. 

Mn 75.4+ 4,27. 
Ni NC 

Zn 64.2+12.17. 



+ ^. 


Tj/. 


± 5. 


67. 


±13. 


57. 


+ 2. 


37. 


+ 3. 


57. 


+ 5. 


57. 


+ 9. 


47. 


+ 1. 


47. 


+65. 


67. 


+ 7. 


97. 



1 Based on the analysis o-f MBS reference samples. 

2 Based on field replicate samples. 

NC = Not calculated because the NBS value and the reference 

samples concentration were below the instrument detection 
1 imi t . 
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Appendix D-2. Vascular plants and bryophytes of the Swamp Gulch 
Wetland study site, 1987. 



VASCULAR PLANTS 
ASTERACEAE 



Greene 
B. & H. 



Greene 



Achillea millefolium L. 
Agoseris aurantiaca (Hook.) 
Anaphalis margaritacea (L.) 
Aster junciforrais Rydb. 
Cirsium arvense (L.) Scop. 
Erigeron peregrinus (Pursh) 
Hieracium scouleri Hook. 
Senecio serra Hook. 
Tragopogon dubius Scop. 

BERBERIDACEAE 

Berberis repens (Lindl.) G. Don 

BETULACEAE 

Alnus incana (L.) Moench 

Betula glandulosa Michx. var Hallii (Howell) Hitchc, 

BRASSICACEAE 

Vicia americana Muhl . 
Thapsis arvense L. 

CARYOPHYLLACEAE 

Cerastium berringianum Cham. & Schlecht. 

CUPRESSACEAE 



Juniper us communis L. 



Carex aquatilis Wahl. 
Carex illota Bailey 
Carex nebraskensis Dewey 
Carex rostrata Stokes 
Scirpus acutus Muhl. 



CYPERACEAE 



ERICACEAE 



Arctostaphylos uva-ursi (L.) Spreng. 
Chimaphila umbellata (L.) Bart. 
Ledum glandulosum Nutt. 
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Appendix D-2. Continued. 



Pyrola asarifolia Michx. 
Vaccinium sp. 



EQUITACEAE 



Equisetum sp. 

Melilotus officinalis (L.) Pallas 

JONCACEAE 

Juncus ensifolius Wikst. 

ONAGRACEAE 

Epilobiuin glaberrimum Barbey 

PINACEAE 

Abies grandis (Dongl.) Forbes 

Abies lasiocarpa (Hook.) Nutt. 

Picea engelmannii Parry ex Engelm. 

Pinus contorta Dongl. ex Lond. 

Pinus ponderosa Dongl. ex Laws. & Laws. 

Pseudotsuga menziesii (Mirbel) Franco 

POACEAE 

Agrostis sp. 

Bromus carinatus H & A 

Bromus ciliatus L. 

Danthonia unispicata (Thumb.) Munro ex Macoun 

Elymus glaucus Buckl. 

Elymus trachycaulus (Link) Gould ex Skinners. 

Festuca idahoensis Elmer 

Muhlenbergia andina (Nutt.) Hitchc. 

Phleum pratense L. 

Scolochloa festucacea (Willd.) Link. 



ROSACEAE 



Fragaria virginiana Miller 
Geum macrophyllum Willd. 
Rosa woodsii Lindl. 



Populus tremuloides Michx. 
Salix boothii Dorn 



SALICACEAE 
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Appendix D-2. Continued. 



SPARGANIACEAE 

Sparganium minimum Fries 

NON-VASCULAR PLANTS 

BRYOPHYTES 

Brachythecium sp. 

Cratoneuron filicinum (Hedw.) Spruce 

Isopterygium pulchellum 

Kopterygium pulchellum (Hedw.) Jaeq. 

Leptobeyum pyriforme (Hedw.) Wils. 

Polytrichum juniperinum Hedw. 

Sphagnum tenellum Ehrh. ex Hoffm. 
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